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Preface 

T he SeaWiFS Project continually strives to improve the quality of the SeaWiFS data products and to expand 
the product suite. This effort requires the SeaWiFS staff to be knowledgeable of recent developments 
in marine optics, atmospheric correction approaches, and vicarious calibration techniques. This capability is 
largely maintained through communication with the user community, particularly the SIMBIOS Project staff 
and science team which provides most of the in situ bio-optical and atmospheric validation data. As a result, 
a number of significant improvements in the sensor calibration, the atmospheric correction scheme, and certain 
masks and flags have been developed since the third reprocessing in mid-2000. 

One particularly important improvement in data quality resulted from the recalibration of the Marine Optical 
Buoy (MOBY) to account for spectrometer stray light. MOBY data is used for the vicarious calibration of 
SeaWiFS. The recalibration required development of a portable version of the National Institute of Standards 
and Technology (NIST) Spectral Irradiance and Radiance Responsivity Calibrations using Uniform Sources 
(SIRCUS) which was deployed to the MOBY facility in Honolulu, Hawaii by NIST staff (Carol Johnson and 
Steve Brown). The cost of the portable system was shared by the SeaWiFS Project, the SIMBIOS Project, 
and the MODIS science team. The recalibration significantly increased the total radiances in bands 1 and 2, 
thereby increasing the retrieved water-leaving radiances and decreasing the frequency of negative water-leaving 
radiances in coastal regions. Another significant improvement is the implementation of a revised near-infrared 
reflectance model in the atmospheric correction scheme. The improved algorithm is designed to work in turbid 
Case-2 water with no impact on corrections in clear waters. 

The research community has always stated a requirement for daily mean PAR, which is needed for computing 
primary production and the surface heat budget. With this reprocessing, daily mean PAR is now included in the 
archive product suite. This PAR product was under development and was tested in collaboration with members 
of the SIMBIOS science team (Robert Frouin and Menghua Wang) for over two years and the algorithm is 
described in Chapt. 8. 

The overall improvements in the derived products, especially in coastal waters is impressive. I would like to con- 
gratulate the SeaWiFS Calibration and Validation Team and others who have helped realize these improvements 
and who have provided validation data. These collaborations underscore the fact that the SeaWiFS Project and 
the science community are partners in the overall SeaWiFS program and the Project looks forward to future 
cooperation. 


Greenbelt, Maryland 
November 2002 


— C. R. McClain 
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Abstract 

The efforts to improve the data quality for the Sea- viewing Wide Field-of-view Sensor (SeaWiFS) data products 
have continued, following the third reprocessing of the global data set in May 2000. Analyses have been ongoing 
to address all aspects of the processing algorithms, particularly the calibration methodologies, atmospheric 
correction, and data flagging and masking. All proposed changes were subjected to rigorous testing, evaluation 
and validation. The results of these activities culminated in the fourth reprocessing, which was completed 
in July 2002. The algorithm changes, which were implemented for this reprocessing, are described in the 
chapters of this volume. Chapter 1 presents an overview of the activities leading up to the fourth reprocessing, 
and summarizes the effects of the changes. Chapter 2 describes the modifications to the on-orbit calibration, 
specifically the focal plane temperature correction and the temporal dependence. Chapter 3 describes the 
changes to the vicarious calibration, including the stray light correction to the Marine Optical Buoy (MOBY) 
data and improved data screening procedures. Chapter 4 describes improvements to the near-infrared (NIR) 
band correction algorithm. Chapter 5 describes changes to the atmospheric correction and the oceanic property 
retrieval algorithms, including out-of-band corrections, NIR noise reduction, and handling of unusual conditions. 
Chapter 6 describes various changes to the flags and masks, to increase the number of valid retrievals, improve 
the detection of the flag conditions, and add new flags. Chapter 7 describes modifications to the level-la and 
level -3 algorithms, to improve the navigation accuracy, correct certain types of spacecraft time anomalies, and 
correct a binning logic error. Chapter 8 describes the algorithm used to generate the SeaWiFS photosynthetically 
available radiation (PAR) product. Chapter 9 describes a coupled ocean-atmosphere model, which is used in 
one of the changes described in Chapter 4. Finally, Chapter 10 describes a comparison of results from the third 
and fourth reprocessings along the U.S. Northeast coast. 


Prologue 

The Sea-viewing Wide Field-of-view Sensor (SeaWiFS) 
Project has continued efforts to improve the overall quality 
and accuracy of the SeaWiFS data products. The algo- 
rithm improvements used through the third reprocessing 
(May 2000) were described in Volumes 9, 10, and 11 of the 
SeaWiFS Postlaunch Technical Report Series , (McClain 
et al. 2000a, 2000b, and O’Reilly et al. 2000, respectively). 
Since that time, work has been ongoing to identify further 
improvements and resolve issues with the data products. 
Specific objectives of the improvements included: reducing 
the number of pixels with negative water-leaving radiance 
( Lwn ) values, particularly in bands 1 and 2; increasing the 
number of valid retrievals and the geographic coverage, by 
identifying overly restrictive mask conditions or logic er- 
rors which exclude otherwise valid data; and improving 
the consistency of the atmospheric correction. Numerous 
other improvements or corrections were identified, for ex- 
ample, to the calibration procedures, bio-optical models, 
and navigation processing. In addition, an entirely new 
product— photosynthetically available radiation (PAR) — 
was added to the product line. 

These efforts culminated in the fourth global data re- 
processing, which was completed in July 2002. Major im- 
provements in the products resulting from the algorithm 
changes are as follows: the fraction of band 1 retrievals 
with negative Lwn values was reduced by half; the total 
number of valid data points was increased by 25%, as was 


the daily geographic coverage; and the monthly geographic 
coverage was increased by 6%. In addition, the agreement 
between the satellite and in situ data was significantly im- 
proved, particularly in the band 1 and 2 matchups and the 
chlorophyll matchups, and the long term consistency of the 
products was also improved. 

A short synopsis of each chapter in this volume is given 
below. 

1. Introduction to the Fourth SeaWiFS Reprocessing 

The fourth reprocessing of the SeaWiFS global data set 
in July 2002 was preceded by over two years of data analy- 
sis, algorithm development, testing, and evaluation by the 
SeaWiFS Project. Numerous issues that were known at 
the time of the third reprocessing were addressed, and so- 
lutions were developed. Additional algorithm refinements 
were proposed during this time. A number of proposed 
changes were deferred, because they either did not show a 
clear improvement in data quality, or were not sufficiently 
well developed to be accepted for operational use. All of 
the accepted changes were subjected to a rigorous, step- 
by-step evaluation, and detailed results were made avail- 
able for review and comment by the scientific community, 
well in advance of the actual reprocessing. This chapter 
describes the activities that culminated in the fourth Sea- 
WiFS data reprocessing. 

2. Changes to the On-Orbit Calibration of SeaWiFS 

Monthly lunar calibrations and daily solar calibrations 

are used to track the on-orbit stability of the radiometric 
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calibration of the SeaWiFS instrument. Analyses of these 
data by the SeaWiFS Calibration and Validation Team 
(CVT) since the third reprocessing showed residual errors 
in the detector temperature corrections, the lunar calibra- 
tion time corrections, and the mirror side reflectance cor- 
rections. These analyses, and the resulting changes to the 
on-orbit calibration methodologies, are presented in this 
chapter. Revised detector temperature corrections allowed 
the lunar calibration time series to be characterized by a 
single, time-dependent exponential function. A time cor- 
rection was applied to the normalization of the lunar cali- 
bration time series; time-dependent mirror side reflectance 
corrections were also implemented. These changes to the 
calibration methodologies enabled the production of a more 
robust on-orbit calibration of SeaWiFS for the fourth re- 
processing of the SeaWiFS mission data set. 

3. Changes to the Vicarious Calibration of SeaWiFS 

The SeaWiFS Project vicariously calibrates the visible 
bands of the instrument against in-water measurements 
from the Marine Optical Buoy (MOBY) in order to achieve 
the necessary calibration accuracy for the SeaWiFS re- 
trievals of water-leaving radiances. The SeaWiFS CVT 
implemented several changes in the vicarious calibration 
procedure for the fourth reprocessing of the SeaWiFS mis- 
sion data set. The most significant change is the use of 
stray light corrected MOBY data in the vicarious calibra- 
tion. More stringent data quality screening criteria were 
also implemented for both MOBY and SeaWiFS data to 
determine the vicarious calibration matchups. Finally, an 
inverse vicarious calibration procedure was implemented in 
which the vicarious gains are computed at the top of the 
atmosphere. The SeaWiFS top-of-the-atmosphere (TOA) 
radiances are compared with MOBY water-leaving radi- 
ances that have been propagated to the top of the atmos- 
phere using the SeaWiFS-retrieved atmospheric correction 
parameters. These changes to the vicarious calibration 
procedure enabled the production of a more robust set of 
vicarious gains for the fourth reprocessing of the SeaWiFS 
mission data set. 

4. Modifications to the SeaWiFS NIR Correction 

This chapter describes several changes made to the 
near-infrared (NIR) correction implemented during the 
fourth reprocessing of SeaWiFS data. The changes made 
were both to correct flaws in the actual implementation 
of the NIR correction and to improve the estimation of 
NIR reflectance. The changes include revised water ab- 
sorption coefficients, the addition of a scaling factor to al- 
low for the gradual introduction of the NIR correction as 
the derived chlorophyll concentration increases, improve- 
ments to the iteration control, and an alternative model 
for the backscatter estimate. The revised NIR correction 
reduces discontinuities in the aerosol model selection at 
the boundaries where the correction is introduced. It also 
reduces the attenuation effect that the original correction 
had on high chlorophyll concentrations. 


5. Changes to the Atmospheric Correction Algorithm 

and Retrieval of Oceanic Optical Properties 

In preparation for the fourth SeaWiFS reprocessing, a 
series of algorithm changes were implemented to enhance 
the performance of the atmospheric correction process and 
to improve the quality and consistency of oceanic opti- 
cal property retrievals. Included in these changes was 
the introduction of a filtering process to reduce the rel- 
ative noise between the two NIR channels. In addition, 
several modifications were made to improve the handling 
of the SeaWiFS out-of-band response, and the normaliza- 
tion of water-leaving radiances was extended to account 
for Fresnel transmittance effects through the air-sea in- 
terface. These and other algorithm updates are described 
within this chapter. 

6. Masks and Flags Updates 

The flags and masks used for the SeaWiFS level-2 and 
level -3 processing were updated for the recent fourth re- 
processing. This chapter discusses the changes and why 
they were made. In many cases, underlying algorithms 
were changed. Some flags changed their states to either 
flagging (noting a condition), or masking (denoting data 
excluded from the product) to allow more data to be kept 
or to improve its quality. New flags were introduced either 
as a part of new algorithms or to denote the status of the 
data more clearly. The flag and mask changes significantly 
contributed to the improvement in the data quality and 
increased the amount of data retrieved. 

7. Level- la and Level -3 Processing Changes 

Improvements were made in the level -la and level-3 
data processing for the fourth reprocessing. Two of these — 
handling of the spacecraft time tag anomalies, and modifi- 
cation of the level -3 space binning program — address spe- 
cific problems in the data and processing logic, which ex- 
cluded otherwise valid data from processing. The third 
improvement, an update to the navigation algorithms, im- 
proves the overall data quality by reducing the maximum 
navigation errors. The net effect of these changes is an in- 
crease both in coverage and quality of the SeaWiFS data 
products. 

8 . The SeaWiFS PAR Product 

A new SeaWiFS derived product was developed which 
provides an estimate for the amount of PAR reaching the 
ocean surface over a 24 hr period. A description of the 
algorithm is provided in this chapter, followed by com- 
parisons with in situ observations. The in situ observa- 
tions include several years of data covering a wide range 
of solar illumination conditions. The results indicate good 
algorithm performance, with root mean squared (RMS) 
differences between satellite-retrieved and observed daily 
average PAR within a few einsteins per square meters per 
day. 
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9. A Partially Coupled Ocean-Atmosphere 
Model for Retrieval of Water-Leaving 

Radiance from SeaWiFS in Coastal Waters 

The global atmospheric correction algorithm for Sea- 
WiFS tends to over correct for the atmosphere in coastal 
waters because of water-leaving radiance (Lw) in the NIR 
part of the spectrum, A r . This Lw(K) phenomenon oc- 
curs particularly in water with high inorganic particulate 
levels. An iterative solution is used to solve this prob- 
lem. A bio-optical model is used to determine the NIR 
backscatter from the backscatter at 670 nm, and specifi- 
cally addresses inorganic particulates. This solution re- 
quires compensation for absorption by chlorophyll, detri- 
tal pigments, and gelbstoff (colored dissolved organic mat- 
ter). The L\v{ A r ) is found and removed from the total 
radiance so that the standard atmospheric models can be 
applied. Chlorophyll a concentrations, C a , in coastal and 
Case- 2 waters are reduced to appropriate levels. The algo- 
rithm cannot yet correct areas where negative Lw occurs 
at 670 nm. 

10. A Comparison of SeaWiFS LAC 
Products from the Third and Fourth 

Reprocessing: Northeast US Ecosystem 

A comprehensive set of 2,524 local area coverage (LAC) 
scenes of the northeast US continental shelf and adjacent 


waters scanned by SeaWiFS from 4 September 1997 to 
11 November 2002 was processed using NASA standard 
algorithms and methods employed in the fourth reprocess- 
ing. Estimates of chlorophyll a concentration (C a ) and nor- 
malized water-leaving radiance ( Lwn ) from the fourth re- 
processing are statistically compared with results obtained 
previously from the third reprocessing. Chlorophyll a con- 
centration from the fourth reprocessing is lower than those 
from the third reprocessing in Gulf Stream and Sargasso 
Sea waters, along the outer continental shelf and slope wa- 
ter, and over the deep basins in the Gulf of Maine. In 
nearshore shelf waters approximately less than 50 m, C a 
from the third and fourth reprocessings are comparable, 
except over the shoal water on Georges Bank, Nantucket 
Shoals, and the northern nearshore Gulf of Maine, where 
C a from the fourth reprocessing is approximately 1.1-1. 2 
times greater than C a from the third reprocessing. The 
median Lwn( 412) and median Lwat(443) values from the 
fourth reprocessing are substantially greater than values 
from the third reprocessing, and the frequency of negative 
water-leaving radiances for the 412 nm and 443 nm bands 
is significantly lower with the fourth reprocessing. Statis- 
tical match-up comparisons between SeaWiFS C a and in 
situ C a indicate that the fourth reprocessing improved the 
accuracy of C a estimates for this region. 
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Chapter 1 


Introduction to the Fourth SeaWiFS Reprocessing 

Gene C. Feldman 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland 

Frederick S. Patt 

Science Applications International Corporation 
Beltsville, Maryland 

Abstract 

The fourth reprocessing of the SeaWiFS global data set in July 2002 was preceded by over two years of data 
analysis, algorithm development, testing, and evaluation by the SeaWiFS Project. Numerous issues that were 
known at the time of the third reprocessing were addressed, and solutions were developed. Additional algorithm 
refinements were proposed during this time. A number of proposed changes were deferred, because they either 
did not show a clear improvement in data quality, or were not sufficiently well developed to be accepted for 
operational use. All of the accepted changes were subjected to a rigorous, step-by-step evaluation, and detailed 
results were made available for review and comment by the scientific community, well in advance of the actual 
reprocessing. This chapter describes the activities that culminated in the fourth SeaWiFS data reprocessing. 


1.1 INTRODUCTION 

Throughout the life of the SeaWiFS mission, the NASA 
SeaWiFS Project, working in close collaboration with 
members of the ocean color research community, has con- 
ducted continual evaluations of the calibration, algorithms, 
and operational procedures used to process SeaWiFS data. 
Significant progress has been made in a large number of ar- 
eas which have not only continuously improved the scien- 
tific quality of the SeaWiFS data products since launch, 
but have also expanded the product suite to include a 
number of valuable additional products. Some of these 
products are the Normalized Difference Vegetation Index 
(NDVT), PAR, and Land Surface Reflectance. Before im- 
plementing any changes, however, to the operational prod- 
ucts that are archived and distributed by the Goddard 
Space Flight Center (GSFC) Distributed Active Archive 
Center (DAAC), the SeaWiFS Project conducts a very 
comprehensive set of evaluations for each of the proposed 
changes, and solicits feedback from the scientific commu- 
nity on the validity of the changes and their effects on the 
derived products. 

In July 2002, the SeaWiFS Project completed the most 
recent round of evaluations and concluded that the im- 
provements justified another reprocessing. To facilitate 
the exchange of information with as broad a representation 
of the international ocean color community as possible, a 
comprehensive Web site was developed that contained all 


the details of the analyses performed, complete statistical 
analyses of the results, and a change-by-change accounting 
for the resulting differences in the derived products. The 
Web sitef contained a summary of the results for each of 
the detailed incremental tests that were run to assess the 
effect of each change proposed for the fourth reprocess- 
ing of the SeaWiFS global data. These tests were all run 
for the month of May 1999. The summary for each test 
included a statement of the changes added for that test, 
and a brief discussion of the changes in the products that 
resulted from the test. The discussions included links to 
text, images, and plots on the Web site which described the 
changes in detail and illustrated the results. In addition 
to these detailed tests, additional analyses were performed 
on the complete SeaWiFS mission period; these results are 
presented on the Web site as well. 

While science is a continuous process, and there will al- 
ways be additional improvements that can be made as un- 
derstanding improves, the SeaWiFS Project believes that 
the results of the most recent evaluations clearly demon- 
strate that better products are produced today than during 
the third SeaWiFS reprocessing. This belief was strongly 
supported by the feedback from the ocean color scientists, 
who reviewed the posted analyses and performed addi- 
tional evaluations of their own. 


f http : //seavif s . gsf c . nasa . gov/SEAWIFS/RECAL/Repro4 
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In parallel with the changes to the processing algo- 
rithms, substantial advances were made in the compu- 
tational capabilities, including both the system hardware 
and the data processing system software. As a result, the 
SeaWiFS Project was able to reprocess all of the stan- 
dard archive products for the entire mission in less than 
one week — achieving a reprocessing rate of approximately 
400:1 throughput capacity. This made possible multiple, 
full mission tests well in advance of the final reprocessing 
and delivery of the archive products to the DA AC, and 
provided both the Project and the scientific community 
with unprecedented knowledge of the full effects of the al- 
gorithm changes. 

1.2 REPROCESSING MOTIVATION 

At the time of the third reprocessing (May 2000), a 
number of issues remained to be resolved with the Sea- 
WiFS products and algorithms. The issues spanned nearly 
all aspects of the level -2 algorithms, including the sen- 
sor calibration, atmospheric correction, in-water retrievals, 
and the flags and masks. In addition, improvements were 
being developed in areas such as navigation, level -3 bin- 
ning, and new products, which could be applied to the 
mission data set only through reprocessing. This section 
summarizes the issues that led to the fourth SeaWiFS re- 
processing. 

1.2.1 Calibration 

The problem of low and negative water-leaving radi- 
ances ( Lwn ) * n SeaWiFS bands 1 and 2 (412 and 443 nm, 
respectively) received a great deal of attention during the 
preparation for the third reprocessing. Indeed, every 
level -2 algorithm change was evaluated with regard to its 
effect on the number of pixels with negative Lwn val- 
ues. Nonetheless, this problem remained largely unre- 
solved, particularly in coastal areas. In addition, the 
match-up comparisons consistently showed that SeaWiFS 
band 1 radiances were low relative to the in situ data, by 
about 15% on average (Bailey 2000). 

The calibration temporal corrections for bands 7 and 
8 were represented as a series of quadratic and linear seg- 
ments, to accommodate periodic fluctuations in the radio- 
metric response. These provided a reasonable representa- 
tion of the instrument response, but required periodic up- 
dating and did not extrapolate well for future processing. 
In addition, the temporal response for all of the bands, de- 
termined from the lunar calibration data, was based on an 
assumption that bands 3 and 4 had no significant temporal 
change. 

The radiometric correction to account for the differ- 
ences of the two sides of the half-angle mirror is referred 
to as the mirror side correction . These corrections also 
showed changes over time, which at that time were cor- 
rected in stepwise fashion, by adjusting the correction con- 
stants for each calibration time segment. This had the 


disadvantage of introducing discontinuities into the mirror 
side corrections. 

The vicarious calibration of the visible bands using the 
MOBY data was performed at the surface, by comparing 
MOBY and SeaWiFS radiances. This required an iterative 
approach, because there was no direct conversion of surface 
radiance differences to TO A calibration corrections. In ad- 
dition, the data selection criteria and processing methods 
were developed independently for the vicarious calibration 
and validation activities, and were not consistent. 

1.2.2 Atmospheric Correction 

The third reprocessing included an iterative method 
of correcting the NIR bands for surface radiance prior to 
their use in atmospheric correction (Siegel et al. 2000). Al- 
though this approach was selected from several candidates 
as providing the best overall results at that time, it suffered 
from several disadvantages: 

a) Small, but measurable effects at very low chloro- 
phyll values, which theoretically should have had 
no effect on these bands; 

b) Underestimation of high chlorophyll values; 

c) Inapplicability to Case-2 waters; 

d) Incorrect absorption coefficients for the SeaWiFS 
bands passes; and 

e) Occasional convergence problems. 

Very clear atmospheric conditions, which resulted in 
low or negative band 8 aerosol reflectance, caused large 
uncertainties in the aerosol model selection or outright fail- 
ure of the atmospheric correction. This frequently resulted 
in chlorophyll retrieval failures under what, in fact, were 
nearly ideal viewing conditions. 

The radiometric resolution of bands 7 and 8 also con- 
tributed to model selection uncertainty (Hu et al. 2000a). 
At low aerosol reflectance levels, one-bit changes in the 
band 7 or 8 signal resulted in significant changes in single- 
scattering aerosol reflectance ratio, e values, which some- 
times resulted in excessive model switching. 

1.2.3 In- Water Retrievals 

The effects of out-of-band contributions to the Sea- 
WiFS radiances were already handled in the third repro- 
cessing by applying a correction to the band-averaged nor- 
malized water-leaving radiances (Lwn)- There were two 
remaining issues: 

1. The solar irradiance used for the computation of 
remote sensing reflectance from the corrected Lwn 
values was still band averaged. 

2. The corrections were determined from a chlorophyll- 
based water-leaving radiance model (Gordon et al. 
1988) and used a log-linear interpolation of 13 dis- 
crete wavelengths to the SeaWiFS band passes. 
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1.2.4 Flags and Masks 

Most of the flag algorithms used for SeaWiFS level -2 
processing were developed prior to the SeaWiFS launch; a 
number had not been reviewed since, or had been superfi- 
cially modified (i.e., refinement of thresholds). The results 
of flagging and masking, specifically the number of pixels 
and bins that were affected, were not well understood. 

1.2.5 Anomalous Chlorophyll 

A small, but significant, number of pixels produced 
anomalous high chlorophyll values in areas of moderate 
or low chlorophyll. These values were averaged into the 
level -3 time binned data products, and continued to be 
visible as speckles , i.e., isolated high values, in the prod- 
ucts even over long compositing periods. 

1.2.6 Other Processing Levels 

The navigation processing performed during level -la 
conversion has undergone continual evaluation and peri- 
odic improvement. These changes are normally incorpo- 
rated into operational processing when they are accepted, 
but can only be applied to the mission data set through 
reprocessing. In addition, the detection and handling of 
time anomalies in the SeaWiFS data stream was entirely a 
manual process, and a problem in the first level -3 binning 
step caused otherwise valid level -2 scenes to not be binned. 

1.2.7 New Products 

An ongoing effort has been to identify additional geo- 
physical parameters that could be produced from SeaWiFS 
data and eventually delivered as archive products. Numer- 
ous candidate products were under consideration as of the 
third reprocessing. The most promising of these were the 
PAR reaching the ocean surface, and the NDVI over land. 

1.3 SOLUTIONS 

The algorithmic changes implemented in response to 
the above stated issues arose from a variety of sources over 
a period of about two years starting in May 2000. Some 
of these were already in progress at that time. Others 
resulted from analysis of specific problems, or were de- 
veloped, based on new information that became available 
during the time period. This section lists the algorithm 
changes (which are described in the other chapters of this 
volume) and the time frame in which they were developed. 

1.3.1 On-orbit Calibration Corrections 

By the time of the third reprocessing, a correlation 
had already been identified between the fluctuations in the 
band 7 and 8 time response and the focal plane tempera- 
ture. This allowed for a much simpler temporal correction; 


however, this modification was not fully developed for op- 
erational use at the time. During the next several months, 
the analysis of the temperature correlation was completed, 
and a new set of temperature corrections was computed to 
replace those based on the prelaunch testing. The tem- 
poral correction was then revised in the form of a (single) 
decaying exponential function covering the entire mission. 
The same functional form was also applied to bands 1-6. 
The implementation of this in a revised calibration table 
was completed by December 2000. The exponential func- 
tion was recomputed periodically using additional lunar 
calibration data after that time; the calibration was final- 
ized just prior to the fourth reprocessing in July 2002. 

The form of the mirror side correction was also revised 
in this time frame. A time-varying correction, in the form 
of a series of linear segments, was implemented in Novem- 
ber 2000. 

In early 2002, analysis of the water-leaving radiances 
from a mission-long test showed that a residual time drift 
remained in the clear- water values. At around the same 
time, analysis of the lunar calibration data indicated a 
measurable temporal response change in bands 3 and 4 
(Kieffer et al. 2002), which were used to normalize the lu- 
nar time series for the other bands. A time correction was 
incorporated into the lunar data normalization, and re- 
vised temporal correction functions were computed for all 
bands. This was completed in March 2002. 

1.3.2 MOBY Stray Light Correction 

One of the most significant improvements to the Sea- 
WiFS calibration arose from the discovery in early 2001 of 
a stray light effect in the MOBY data calibration (Clark 
et al. 2001). It was quickly determined that this effect had 
caused the measurements at the blue end of the spectrum 
(bands 1 and 2) to be biased low. Because all of the vis- 
ible bands are vicariously calibrated using MOBY as the 
reference, this stray light effect was largely responsible for 
the observed radiance bias in the in situ match-up com- 
parisons, and also contributed substantially to the number 
of negative Lwn values. 

The first set of corrected MOBY measurements were 
received in February 2002. Several additional corrected 
measurement sets were received prior to the fourth repro- 
cessing, and were incorporated into the final vicarious cal- 
ibration of the visible bands. 

1.3.3 Vicarious Calibration Procedures 

In addition to the use of the stray light corrected MOBY 
data, significant changes were made to the vicarious cali- 
bration procedures for the fourth reprocessing. The new 
level -2 processing software, developed for the third repro- 
cessing, made possible the comparison of radiances at the 
top of the atmosphere (the inverse calibration ) instead of 
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the ocean surface (forward) as before. This allowed the cal- 
ibration to be performed in a single step, instead of itera- 
tively, significantly reducing the processing time required. 
The inverse calibration procedure was already under de- 
velopment as of the third reprocessing, and was shown to 
produce consistent results with the forward calibration. 

In August 2001, a detailed comparison was performed 
of the selection and processing procedures for the vicari- 
ous calibration and in situ match-up comparisons, and a 
consistent set of procedures was developed. Additional se- 
lection criteria for the MOBY data were developed in April 
and May of 2002, and were used for the final vicarious cal- 
ibration. 

1.3.4 Flags and Masks 

A detailed evaluation was undertaken of the flags used 
as masks for the third reprocessing, beginning in Septem- 
ber 2001. Specifically, the objective was to determine the 
number of observations that were affected by each flag, 
and whether the flag algorithm could be modified to re- 
duce the number of affected pixels without affecting data 
quality. This analysis, which continued through December 

2001, resulted in the changes made to the stray light flag 
and the sensor zenith angle limit. 

The cloud and ice flag was subjected to multiple chan- 
ges. In conjunction with the temporal calibration change of 
band 8, the threshold was adjusted to compensate for the 
change in average radiance. A correction of the algorithm 
for air mass was proposed in September 2000. In January 

2002, however, a substantial revision was proposed to the 
cloud and ice algorithm: use of the band 8 surface and 
aerosol reflectance in place of TOA reflectance. This re- 
moved the effect of increasing Rayleigh radiance at high 
solar zenith angles, which caused the flag to be set for 
valid, open ocean pixels at high latitudes. 

The absorbing aerosol flag was analyzed starting in July 
2000, and this revealed some implementation issues that 
were causing inconsistent performance of this flag. In ad- 
dition, the thresholds were adjusted to compensate for the 
vicarious calibration changes in bands 1 and 2. 

During the final evaluation of the algorithm changes, it 
was shown that the speckles were still present in the data, 
and the revised cloud flag had increased the incidence of 
the speckles at high latitudes. Analysis of these observa- 
tions showed that many of them exhibited distinctly non- 
physical spectral behavior. As a result, a spectral test was 
added to the algorithm for the chlorophyll failure flag. 

Additional flag changes were made as follows: 

• The turbid water flag was revised over the period 
from November 2001-January 2002, to be based on 
a threshold applied to the band 6 remote sensing 
reflectance. 

• The high total radiance (high L t ) flag, which was 
set if any band had a value above the knee in the 
bilinear gain, was modified in January 2002 to apply 
this test only to bands 7 and 8. 


• The high glint flag was removed as a mask; this 
change had already been suggested, but not evalu- 
ated, by the time of the third reprocessing. 

1.3.5 NIR Correction 

The effect of the NIR surface radiance correction at low 
chlorophyll levels was known as of the third reprocessing. 
A modification was already being designed to phase in the 
correction around a chlorophyll level of lmgm“ 3 . This 
modification was completed by the end of 2001. 

A change to the NIR absorption coefficients, to better 
correspond to the SeaWiFS band passes, was introduced in 
November 2001. In December 2001, the NIR backscatter 
model was changed from chlorophyll- to reflectance-based, 
at a reference wavelength of 670 nm. Improved control 
of the NIR correction iteration process was completed in 
February 2002. 

1.3.6 NIR Noise Reduction 

A method to reduce the effect of NIR band digitiza- 
tion error on the atmospheric model selection was already 
developed as of the third reprocessing, but had not been 
accepted for operational use at that time. The general ap- 
proach of filtering radiances was developed starting in Oc- 
tober 1998. It was substantially enhanced in June 2000, 
and was also applied to the Ocean Color and Temperature 
Scanner (OCTS) global data processing by the NASA Sen- 
sor Intercomparison and Merger for Biological and Inter- 
disciplinary Oceanic Studies (SIMBIOS) Project in 2001. 
After extensive analysis of the results, it was decided in 
June 2002 to apply this algorithm only to the 1 km LAC 
and high resolution picture transmission (HRPT) data, 
and not the global area coverage (GAC) data. 

1.3.7 Miscellaneous Changes 

A number of other changes were also made: 

1. A modification was developed in April 2001, to use 
different solar irradiance values for calculations be- 
fore and after the out-of-band correction was ap- 
plied to the Lwn calculation. 

2. The atmospheric model selection algorithm was 
modified to revert to a fixed model for very clear 
atmospheres (i.e., band 8 aerosol reflectance below 
a threshold). 

3. In June 2002, a correction to the Lwn calculation 
was implemented for the Fresnel reflectance at the 
air-water interface. 

4. Under certain geometric conditions, certain aerosol 
models cross over in e space, i.e., when plotting e 
as a function of scattering angle, causing discon- 
tinuities when the aerosol path radiances are ex- 
trapolated into the visible. A correction was de- 
veloped, also in June 2002, which identifies these 


7 



Algorithm Updates for the Fourth SeaWiFS Data Reprocessing 


model crossing conditions and revises the model se- 
lection result accordingly. 

5. Storing chlorophyll in the level -2 product as a scaled 
integer introduced artifacts at low chlorophyll val- 
ues. The solution was to store the level -2 chloro- 
phyll value as a floating point variable instead of an 
integer. 

1.4 TESTING AND EVALUATION 

A nearly final set of algorithm changes was implemented 
by December 2001. At that time, work began on a strategy 
for systematic testing and evaluation of the results, with 
many opportunities for review and feedback by the scien- 
tific community. This section presents the overall testing 
approach, the specific test cases run, and a summary of 
the results. 

1.4.1 Testing Approach 

The month of May 1999 was selected for the case-by- 
case evaluation of the algorithm changes. For each case, a 
number of items were generated: 

■ Images of daily, 8-day, and monthly binned prod- 
ucts; 

■ Difference and ratio images; 

■ Difference and ratio histograms; 

■ Two-dimensional histograms; 

■ Sample and bin gain-or-loss statistics; and 

■ A summary of the most significant effects. 

These results were posted on the SeaWiFS Reprocessing 
Web site, and the scientific community was invited to re- 
view and comment on the results. 

In addition, multiple, mission long tests were run, and 
were used to generate comparisons of in situ matchups, 
clear water and deep water analyses, and mapped images. 
All of these results were posted on the Web site as well. 

1.4.2 Test Cases 

The following describes the test cases that were run. 
Initially, a baseline run was generated using the configura- 
tion for the third reprocessing. Each case was then com- 
pared with the previous one to generate the statistics and 
difference images. The vicarious gains were recomputed 
for several cases as required. 

Case 1 , MOBY Stray Light Correction: The new visible 
band vicarious gains were derived from the MOBY data 
corrected for stray light effects, along with refinements to 
the vicarious calibration procedures. The configuration 
was otherwise identical to that for the third reprocessing. 

Case 2, Updated Calibration Table: The new calibra- 
tion table implements changes in the calibration time de- 
pendence (from piecewise quadratic to exponential), tem- 
perature corrections and mirror side corrections for all 


bands. The change also included a re-evaluation of the 
vicarious gains in bands 1-7, and an adjustment to the 
cloud mask threshold to compensate for a reduction in the 
band 8 calibrated radiance. 

Case 3 , Revised Out-of-Band Correction : The revised 
out-of-band correction is applied prior to use in down- 
stream computations (e.g., chlorophyll retrieval). A re- 
vised set of correction factors was generated using the re- 
cently published clear-water reflectance model (Morel and 
Maritorena 2001). 

Case 4 , Miscellaneous Changes and Corrections: The 
high L t flag, which was previously set when any band ex- 
ceeded the bilinear calibration knee value, was changed to 
only test the NIR bands, where the increased quantization 
of the radiance can significantly affect aerosol retrievals. 
The high L t flag was applied as a mask at level -3. In ad- 
dition, the threshold for the high sensor zenith angle flag 
was increased from 56° to 60°. A new turbid water flag 
was introduced, but this flag is only informational. This 
stage in the test sequence also introduced the improved 
handling of low aerosol conditions, a modified NIR iter- 
ation control, and the phase-in of the NIR correction. In 
addition, the transmittance in the cloud albedo calculation 
was modified to improve behavior at high solar and view- 
ing angles, and the NIR water-leaving radiance model was 
updated with new absorption and backscatter coefficients. 
The cloud albedo threshold was adjusted to compensate 
for the algorithm change. The cloud albedo changes were 
superseded by subsequent modifications (Case 7). 

Case 5, Change in NIR Correction Method: The NIR 
correction changed from chlorophyll-based to reflectance- 
based. This correction, as implemented in the third re- 
processing, tended to depress the 510-865 nm Angstrom 
values, which often caused the selection of oceanic aerosol 
models and resulted in reduced chlorophyll retrievals. 

Case 6 , Change to Glint Mask Use: This case incorpo- 
rated the change to eliminate masking based on the high 
glint flag. At the third reprocessing, a glint correction was 
introduced for pixels where glint radiance was less than the 
threshold set to trigger the high glint flag. During prepara- 
tion for the fourth reprocessing, it was demonstrated that 
this correction was reasonable up to the level where the 
high L t flag is triggered. 

Case 7, Improved Cloud Flagging Method: The new 
cloud flagging algorithm was based on band 8 surface- and 
aerosol reflectance, which better compensates for the in- 
creasing Rayleigh path radiance with solar and viewing an- 
gles. The third reprocessing algorithm, based on the TOA 
reflectance, was found to be overly restrictive at higher 
solar and viewing angles. 

Case 8 , Reduce Stray Light Mask: The stray light 
masking algorithm was changed to reduce the number of 
stray light masked pixels from three to two before and af- 
ter a bright target. The original masking correction scheme 
for stray light (GAC resolution), determined before launch, 
was conservative with the intention that the masking and 
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Table 1. The following table shows a summary of the median percent change for each incremental change to the 
SeaWiFS processing code relative to the previous change. The column headers “M” and “D” denote “Monthly” and 
“Daily,” respectively. The “Total” results axe versus the third reprocessing. 


Case 

Filled Bins 
M D 

Samples 

M D 

c a 

M 

D 

Lwn( 412) 
M D 

Lwn( 490) 
M D 

£ww(555) 

M D 



r a (865) 
M D 

i 

0.6 

1.9 

2.4 

2.8 

-6.0 

-5.2 

17.0 

16.5 

6.2 

6.2 

4.7 

4.8 

0.0 

0.0 

2 

-0.1 

-0.4 

-0.5 

-1.0 

1.1 

1.1 

-1.8 

-1.7 

-1.4 

-1.3 

-1.0 

-0.9 

-3.3 

-3.4 

3 

0.0 

-0.2 

-0.3 

-0.3 

-3.7 

-3.8 

2.3 

2.3 

-0.2 

-0.2 

-1.3 

-1.4 

0.0 

0.0 

4 

0.8 

3.4 

3.0 

3.3 

-1.4 

-1.5 

-0.1 

-0.1 

0.6 

0.5 

-0.3 

-0.5 

0.8 

0.9 

5 

0.0 

-0.1 

-0.1 

-0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

6 

0.2 

0.9 

0.8 

1 0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

7 

1.4 

1.5 

1.5 

1.3 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

8 

3.0 

16.0 

16.0 

16.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

9 

0.1 

0.4 

0.4 

0.5 

3.2 

3.2 

-0.6 

-0.6 

-0.6 

-0.6 

1.0 

1.1 

-1.2 

-1.1 

Total 

6.0 

1 

24.7 

24.6 

24.3 

-7.1 

-7.6 

17.3 

16.7 

5.1 

4.7 

3.6 

2.9 

-1.3 

-3.4 


correction schemes be revisited after launch. The change 
is only implemented for GAC resolution data. Stray light 
corrections were applied to the additional unmasked pix- 
els. 

Case 9: Add Fresnel Transmittance Correction , Chloro- 
phyll Failure Condition Modification , and Final Gains: The 
normalization of water-leaving radiance was modified to 
include a correction for Fresnel transmittance through the 
water-atmosphere interface. In addition, a few other mi- 
nor changes were introduced at this stage. 

The number of match-up points used in the vicarious 
calibration was further reduced by the exclusion of six 
points with moderate glint contamination. The effect of 
removing these points was to very slightly raise the gains 
(all things being equal). Additionally, four more lunar cal- 
ibration points were added to the data set used to compute 
the time dependence correction for the calibration table. 
This case also included the fix for the model crossing prob- 
lem. 

1.4.3 Cumulative Effects of Changes 

Table 1 shows a summary of the median percent change 
for each incremental change to the SeaWiFS processing 
code relative to the previous change. The statistics pre- 
sented include the change in the number of filled bins 
and the total number of valid samples; the median change 
in chlorophyll; the Lwn values at 412, 490, and 555 nm 
(bands 1, 3, and 5, respectively); and the aerosol optical 
thickness (AOT, r G ) at 865 nm (band 8). The statistics 
were computed for monthly and daily binned files, using 
the May 1999 test runs. 

The statistics show that the most significant changes 
were an increase of nearly 25% in the number of valid data 
samples (with an equivalent change in the number of filled 
bins in the daily product), and an increase of about 17% 
in the 412 nm (band 1) Lwn values, mainly due to the im- 
provement in the vicarious calibration from the MOBY 


stray light correction. Also significant is the fact that 
Cases 6, 7, and 8, which were intended only to increase 
the number of valid retrievals, had no measurable effect on 
radiometry. 

1.5 OTHER OPTIONS NOT ACCEPTED 

During the period between the third and fourth Sea- 
WiFS data reprocessings, a number of algorithm changes 
were considered, but were either deferred for further anal- 
ysis or rejected. These are summarized below. 

Revised At-launch Gains: Prior to the SeaWiFS launch, 
two separate calibrations were performed, one by Santa 
Barbara Research Systems (SBRS) in 1993 and the other 
by the National Institute of Standards and Technology 
(NIST) in 1997. The NIST 1997 gains have been used 
for all processing to date, modified by the vicarious gains 
determined using the MOBY data (bands 1-6) and atmos- 
pheric model assumptions near MOBY (band 7). Recent 
analysis has shown that the vicarious gains are more con- 
sistent with the SBRS 1993 calibration. For ocean data, 
the vicarious calibration largely cancels out any effect on 
the bands 1-6 Lwn values, but the reflectances would be 
affected, and vicarious gains are not used for land data 
processing. The effect of this change was not evaluated in 
time for the fourth reprocessing. 

NIR Band Alternative Calibration Schemes: The cali- 
bration for the NIR bands has been based on the prelaunch 
gain for band 8 since launch, and the vicarious calibration 
based on atmospheric model assumptions for band 7 since 
shortly after launch. Several alternative approaches were 
investigated, including use of a single model, or selection of 
models based on in situ AOT measurements, for the band 
7 vicarious calibration; vicarious calibration of both bands 
7 and 8 using the AOT measurements; and incorporation 
of local wind speed measurements at the MOBY location. 
None of these were found to produce improvements over 
the current approach, although the AOT-based vicarious 
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calibration improved the AOT matchups performed by the 
SIMBIOS Project. 

Alternate Cloud Shadow Algorithm: The SeaWiFS low 
water-leaving radiance flag is used to detect cloud shadows 
based on a threshold on the band 5 Lwnv alues. An algo- 
rithm based on multiple bands might be more robust, but 
this was not investigated for the fourth reprocessing. 

Revised Solar Irradiance Model: The current solar irra- 
diance model used in the algorithms is based on the work of 
Neckel and Labs (1984). Other projects such as the Moder- 
ate Resolution Imaging Spectroradiometer (MODIS), have 
switched to a more recent model developed by Thuillier 
et al. (1998). Thuillier et al. have also been working on 
additional refinements to the model, and a new version 
has been developed (Thuillier et al. 2003). A decision was 
made to postpone implementation of the so-called Thuil- 
lier model until the latest spectrum is officially published 
and the effects on SeaWiFS products are fully quantified. 

Alternative Aerosol Model Selection Algorithm: The 
current aerosol model selection algorithm is designed to 
choose 2 models from the full set (currently 12) and inter- 
polate linearly between them [the ratio of band i to band 
j is hereafter denoted as K/j], based on the ratio of band 
7 to band 8 (A 7 / 8 ) single-scattering aerosol reflectance (e) 
value. This approach has problems under several com- 
monly encountered conditions: 

a. Two models have nearly identical NIR band ratios 
for certain geometric conditions, though they differ 
in the visible bands. 

b. The aerosol reflectance values for the visible bands 
change discontinuously when the model pairs 
change. 

c. The uncertainty in the e value is large for small NIR 
band values, causing unnecessary model switching. 

Alternative model selection algorithms were discussed 
to reduce the effect of these problems, by using more than 
two models. One scheme proposed was a gaussian weigh- 
ing of multiple models; tests of this did not work well, and 
the suspected cause was a nonuniform distribution of the 
models in e space. Another approach has been proposed 
which would involve a weighted linear fit to multiple mod- 
els in the vicinity of the computed e value. Investigation 
of this approach was in progress at the time of the fourth 
reprocessing. A simple fix for a specific model selection 
problem was described in Sect. 1.3.7, item 4. 

NO 2 Correction: At the fifth SIMBIOS Science Team 
Meeting (January 2002), it was suggested that NO 2 ab- 
sorption may be significant. The current atmospheric cor- 
rection algorithm does not consider NO2 absorption, how- 
ever, the topic bears future study. 

Pressure Dependence of the Band 7 Oxygen Absorption 
Correction: The current oxygen absorption correction for 
band 7 is independent of atmospheric pressure. The effect 
of pressure on this correction is small and likely within 
the uncertainty of the correction. The effort involved was 
determined to be unwarranted for the fourth reprocessing. 


f/Q Correction: The SeaWiFS Project has been work- 
ing with D. Antoine and A. Morel on the issue of the bidi- 
rectional reflectance at the water’s surface, f/Q, since Oc- 
tober 1999. The approach has been to evaluate the effect 
of the correction on reducing viewing angle dependencies 
in the Lwn values. Thus far, the algorithm has not been 
found to improve the quality of the SeaWiFS products. 
The results of the f/Q correction have been provided to 
Antoine and Morel to continue the investigation. 

Revert to At-Launch Kd( 490) Algorithm : The 1^(490) 
(diffuse attenuation coefficient at 490 nm) algorithm was 
changed with the third reprocessing, from using a ratio of 
A 2 /5 to a ratio of A3/5. This change was made because of 
uncertainties in the band 2 radiances. With the calibration 
now based on the stray light corrected MOBY data, the un- 
certainties in the band 2 radiances have been reduced. The 
at-launch algorithm for LQ(490) is more sensitive, partic- 
ularly for the lower values found in open ocean conditions. 
For these two reasons, it was proposed that the at-launch 
algorithm be re-implemented. Preliminary analysis showed 
there was general agreement between the current algorithm 
and the previous version. In turbid waters, however, band 
2 radiances become negative more often than band 3. This 
causes the A 2 /s algorithm to produce unrealistically high 
values. For this reason, the decision was made to continue 
with the A3/5 algorithm for this reprocessing. 

Log-Scaled Chlorophyll : A proposal to replace the chlo- 
rophyll values with their logarithm in the level -2 product 
was considered, to account for the log-normal distribution 
of chlorophyll. This was initially a response to the arti- 
facts resulting from the use of scaled integers (Sect. 1.3.7). 
It also reduced the effect of anomalous chlorophyll values 
(Sect. 1.2.5). The change was determined to be unaccept- 
able because of the effect on global chlorophyll statistics, 
and the floating point storage solution was implemented 
instead. 

1.6 CONCLUSIONS 

The SeaWiFS Project undertook the preparation for 
the fourth reprocessing with a number of issues to resolve. 
More importantly, perhaps, the Project decided from the 
outset that each change would be accompanied by a techni- 
cal justification; the effects of each change would be clearly 
demonstrated by testing, analysis, and evaluation; and the 
scientific community would be engaged to the maximum 
extent possible throughout. Some changes were already 
in progress as of the third reprocessing, while others were 
proposed when the effort was nearly completed. Most were 
developed from analyses conducted by Project members, 
but a few were developed from work done for other sen- 
sors, and one (e.g., the MOBY stray light correction) came 
about as a result of external events that could not have 
been foreseen in advance. Several promising ideas were 
ultimately shown to produce no tangible benefit. 
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The results posted from the tests and evaluations show 
substantial improvement in overall data quality in a num- 
ber of areas compared to the third reprocessing. Nonethe- 
less, the number of potential changes deferred for further 
analysis is a clear indication that further improvements are 
possible. In addition, some issues (e.g., negative Lwn val- 
ues and absorbing aerosols) were only partially addressed 
by the improvements described in this volume. Future 
changes will be implemented following the same approach 
to data analysis, algorithm development, testing, and eval- 
uation as summarized in this chapter. 
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Chapter 2 


Changes to the On-Orbit Calibration of SeaWiFS 

Robert E. Eplee, Jr., Robert A. Barnes, 
and Frederick S. Patt 
Science Applications International Corporation 
Beltsville, Maryland 

Abstract 

Monthly lunar calibrations and daily solar calibrations are used to track the on-orbit stability of the radiometric 
calibration of the SeaWiFS instrument. Analyses of these data by the SeaWiFS CVT since the third repro- 
cessing showed residual errors in the detector temperature corrections, the lunar calibration time corrections, 
and the mirror side reflectance corrections. These analyses, and the resulting changes to the on-orbit calibra- 
tion methodologies, are presented in this chapter. Revised detector temperature corrections allowed the lunar 
calibration time series to be characterized by a single, time-dependent exponential function. A time correction 
was applied to the normalization of the lunar calibration time series; time-dependent mirror side reflectance 
corrections were also implemented. These changes to the calibration methodologies enabled the production of 
a more robust on-orbit calibration of SeaWiFS for the fourth reprocessing of the SeaWiFS mission data set. 


2.1 INTRODUCTION 

The SeaWiFS CVT undertook extensive analyses of 
the SeaWiFS on-orbit calibration data and ocean data in 
preparation for the fourth reprocessing of the mission data 
set. These analyses showed three sources of residual errors 
in the on-orbit calibration data: a) an annual periodicity 
in the lunar calibration time series for some bands, b) a 
residual time drift in the lunar calibration data, and c) a 
change in the half-angle mirror side reflectances with time. 
A number of changes were implemented to the on-orbit cal- 
ibration methodologies to address these residuals. These 
changes, which have resulted in a more robust on-orbit cal- 
ibration of SeaWiFS, are detailed in the following sections 
of this chapter. 

2.2 TEMPERATURE CORRECTION 

The ambient temperatures at which the SeaWiFS de- 
tectors operate are affected by the solar insolation on the 
spacecraft. The SeaWiFS calibration incorporates correc- 
tions for variations in the radiometric response of the in- 
strument with changing detector temperatures (Barnes et 
al. 1994). The prelaunch corrections do not fully account 
for the variation in radiometric response with changing 
detector temperatures observed over the course of a year. 
Examination of the SeaWiFS lunar calibration time se- 
ries (Fig. 1) shows an annual periodicity in some bands 
which corresponds to the variation in detector temper- 
atures (Fig. 2) as the Earth-sun distance changes over 


a year. Through the third reprocessing, the CVT dealt 
with these annual periodicities with a series of piecewise 
quadratic and linear segments in the time correction ap- 
plied to the ocean data (Barnes et al. 2001 and Eplee and 
Barnes 2000). This approach required regular updates to 
the operational calibration table. 

For the fourth reprocessing, the periodicities in the lu- 
nar time series were used to compute revised detector tem- 
perature correction factors for the SeaWiFS calibration. 
The SeaWiFS level-lb calibration equation includes a cor- 
rection for variations in the radiometric response of the de- 
tectors as a function of temperature, fp (Eplee and Barnes 
2000). This correction has the form 

MX) = 1 + K 2 (X)(T-T te{ ), (1) 

where A is the wavelength, K 2 is the temperature correc- 
tion factor, T is the detector temperature from the sensor 
output data, and T re f is the detector reference temperature 
(20°C). 

During the prelaunch calibration of SeaWiFS, the tem- 
perature correction factors were derived for each band over 
a limited range of detector temperatures. These values are 
shown in Table 2. To investigate the periodicities in the 
lunar calibration time series, each band in the time series 
was fit to a decaying exponential function of time. The 
residuals of these fits showed linear dependencies on de- 
tector temperature. The resulting slopes of the fit residual 
versus detector temperature trends were used to revise the 
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Fig. 1. The initial lunar calibration time series. The lines are fits to the data. 



Fig. 2. The band 7 and 8 detector temperatures during lunar calibrations. 


temperature correction factors. As shown in Table 2, the 
changes in the correction factors varied with band, with 
the largest changes affecting bands 6-8. 


Table 2. The detector temperature correction fac- 
tors are shown. 


Band 

Prelaunch K 2 

[(°c) -1 ] 

Revised K 2 
[ ( °C)-1] 

1 

9.01 x 10~ 4 

7.664 x 10~ 4 

2 

5.85 x lO" 4 

5.540 x 10- 4 

3 

4.20 x 10" 4 

3.392 x 10" 4 

4 

3.90 x lO" 4 

3.057 x 10" 4 

5 

3.91 x 10" 4 

3.045 x 10" 4 

6 

1.51 x 10" 4 

-3.443 x 10~ 5 

7 

1.06 x 10“ 4 

-4.495 x 10" 5 

8 

7.8 x 10" 5 

-1.485 x 10“ 3 


These revised temperature correction factors were used 
to reprocess all of the lunar calibration data. The resulting 
time series (Fig. 3) no longer showed the annual periodici- 
ties. The revised lunar calibration time series was fit with 
a single decaying exponential function of time, /(t), for 
each band. This function has the form 

f(X,t) = ao(A) - oi(A) [l - (2) 

where t is the time of the calibration, to is the reference 
time for the time series, a 0 is the initial value of the func- 
tion, ai is the difference between a 0 and the asymptotic 
value of the function, and a 2 is the time constant of the 
function. 

As time increases, the function approaches the asymp- 
totic value a 0 - a\. These time correction functions were 
incorporated into the revised calibration table used in the 
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fourth reprocessing. Unless a significant change in the be- 
havior of the radiometric response of the detectors occurs, 
piecewise updates to the calibration table should not be 
necessary for operational processing of the ocean data. 

2.2.1 Effects Over One Orbit 

The effects of the revised detector temperature correc- 
tions can be observed over a single orbit. As part of the 
the fourth reprocessing analysis effort, these effects were 
examined for one orbit on 1 May 1999 for bands 7 and 8. 
Bands 7 and 8 reside on a single focal plane, for which the 
instrument telemetry provides one focal plane temperature 
measurement. 

The detector temperature for bands 7 and 8 for a typ- 
ical orbit on 1 May 1999 is shown in Fig. 4 as a function 
of the latitude of the center pixel of the scan line. The 
digitization of the temperature is 0.2°C. The behavior of 
the other bands is similar. During the part of the orbit 
when science data are not being collected, the so-called 
back orbit, the detectors are turned off and the focal plane 
heater is turned on. As the spacecraft crosses over the 
North Pole, the focal plane heater is turned off and the 
detectors are turned on. Typically, the focal plane heater 
keeps the detectors at a higher temperature than that of 
their normal operation. As shown in Fig. 4, the detector 
temperature falls from the initial temperature of approx- 
imately 16.7°C to a more typical temperature of approx- 
imately 14.1°C within about the first 5° of latitude. The 
temperature continues to fall to approximately 13.7° C over 
the next 15° of latitude. This temperature, which is the 
nominal operating temperature of the detectors for this 
orbit, is stable through 55° of latitude. As the spacecraft 
approaches the equator, solar heating raises the tempera- 
ture to approximately 13.9°C, and then to approximately 
14.1°C. For most of the orbit, from 75°N to 60°S, the de- 
tector temperature is stable to within two counts (0.4° C). 

The detector temperature corrections for bands 7 and 8 
are shown over the orbit as functions of latitude in Fig. 5. 
The prelaunch corrections are plotted as dashed lines and 
the revised corrections are plotted as solid lines. The vari- 
ation in the prelaunch corrections over the orbit is essen- 
tially unity. Over the temperature range of 14-17°C, the 
revised temperature correction for band 8 is approximately 
0.05% per count and the correction for band 7 is approx- 
imately 0.02% per count. The change in the temperature 
correction for band 8 is approximately 0.1% over the course 
of the orbit. Because the band 8 radiance is used in the 
atmospheric correction algorithm to determine the aerosol 
abundance (Eplee et al. 2001, and Robinson and Wang 
2000), the decrease in radiance in the Southern Hemisphere 
may affect the retrieved aerosol abundance between the 
third and fourth reprocessings for a single pixel. 

The combined temperature and time corrections affect 
the total radiances. In this analysis, single time corrections 
for bands 7 and 8 were computed for 1 May 1999 and were 


used to scale the detector temperature corrections over the 
orbit. The A 7 / 8 ratio, normalized by the vicarious gain 
of band 7, was used to determine the aerosol type in the 
atmospheric correction algorithm (Eplee et al. 2001 and 
Robinson and Wang 2000). [The gain of a given band is 
hereafter denoted as G* , e.g., G 7 is the vicarious gain for 
band 7.] Figure 6 shows the changes in the normalized 
radiance ratio (G 7 x A 7 / 8 ) between the third and fourth 
reprocessings as a function of latitude. The plot shows 
that for the Southern Hemisphere the changes in the radi- 
ance ratio between the third and fourth reprocessings may 
cause changes in the retrieved aerosol type. The changes 
in aerosol retrievals between the third and fourth repro- 
cessings, shown in Figs. 5 and 6, will propagate through 
to changes in the retrieved normalized water-leaving radi- 
ances and chlorophyll concentrations. 

2.2.2 Effects Over Time 

The effects of the revised detector temperature cor- 
rections and resulting exponential time correction are ob- 
served on a monthly basis in the lunar calibration data. 
As part of the fourth reprocessing analysis effort, these ef- 
fects were examined daily for bands 7 and 8, where the 
temperature correction revisions are the most significant. 

Figure 7 shows the mean daily detector temperature 
for bands 7 and 8 over the course of the mission. The tem- 
peratures range from 14-20°C on an annual cycle as the 
Earth-sun distance changes. The combined temperature 
and time correction to the radiances for band 8 over the 
course of the mission are shown in Fig. 8. The prelaunch 
temperature correction and quadratic time correction (the 
third reprocessing) are shown as dashed lines and the re- 
vised temperature correction and exponential time correc- 
tion (the fourth reprocessing) are shown in solid lines. Fig- 
ure 9 shows the corrections for band 7. For both bands, 
the effect of the revised temperature correction on the com- 
bined radiance correction is apparent. 

Bands 7 and 8 are used for atmospheric correction of 
the ocean data. The band 8 radiance is used to deter- 
mine the aerosol abundance, and the L(A 7 / 8 ) value (nor- 
malized by G7), is used to determine the aerosol type. 
Figure 10 shows the changes in the normalized radiance 
ratio [G 7 x L( A 7 / 8 )] between the third and fourth repro- 
cessings as a function of time. At any point during the 
mission, changes in the band 8 radiance and in the radi- 
ance ratio between the third and fourth reprocessings may 
result in changes in the retrieved aerosol abundance and 
type for a given pixel. These changes in the aerosol re- 
trievals will propagate through to changes in the retrieved 
normalized water-leaving radiances and chlorophyll con- 
centrations. 

2.3 LUNAR NORMALIZATION 

For one analysis of a fourth reprocessing mission test, 
mean global Lwn values in bands 1-6 were computed from 
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Fig. 3. The revised lunar calibration time series. The lines axe fits to the data. 



Fig. 4. Band 7 and 8 detector temperatures are plotted over one orbit from 1 May 1999. 



Fig. 5. The band 7 and 8 detector temperature corrections are plotted over the orbit. The two lower lines 
show the initial corrections. The two upper lines show the revised corrections. 
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Fig. 6. Relative changes in the L( X 7 /s) values versus latitude resulting from the new temperature correction 
in the fourth reprocessing. 



Day of Year 

Fig. 7. The band 7 and 8 detector temperatures are plotted over the mission. 



Day of Year 

Fig. 8. The band 8 combined time and temperature corrections are plotted over the mission. The dashed 
line shows the initial corrections. The solid line shows the revised corrections. 
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Day of Year 

Fig. 9. The band 7 combined time and temperature corrections are plotted over the mission. The dashed 
line shows the initial corrections. The solid line shows the revised corrections. 



Day of Year 

Fig. 10. Relative changes in the L( A 7 / 8 ) values versus time resulting from the new temperature correction 
in the fourth reprocessing. 


the 8-day composite files for clear water (water with a 
depth of at least 1 km and with a chlorophyll concentra- 
tion of less than 0.15 mgm~ 3 ). The mean radiances in clear 
water should be constant with time if the radiometric cal- 
ibration of the instrument is stable (Eplee et al. 2001, and 
Eplee and McClain 2000). This analysis showed decreases 
in the radiances for bands 1-6 of several percent over the 
course of the mission. Figures 11 and 12 show the clear 
water time series for bands 3 and 5 as dashed lines. Band 
3 exhibited a change of —2% over 1,000 days and —3.4% 
over the mission (a rate of change of —1.9 x 10~ 5 per day). 
Band 5 exhibited a change of -2.6% over 1,000 days and 
—4.4% over the course of the mission (a rate of change of 
—2.5 x 10 -5 per day). These results imply that there was a 
small residual decrease in the radiometric response of Sea- 


WiFS with time that was not corrected by the previous 
lunar calibration-based time corrections of the SeaWiFS 
data. With 4.5 years of data available, small residual drifts 
can become significant . 

The CVT, in analyzing the lunar data to derive the 
corrections for the instrument, references the lunar cal- 
ibration time series to the mean of the observations for 
bands 3 and 4 (Barnes et al. 2001, and Eplee and Barnes 
2000). This normalization is done to reduce the systematic 
noise in the lunar observations which arises from the lack 
of understanding of how lunar libration affects the long 
term measurements of the moon. The initial assumption 
was that the mean of bands 3 and 4 did not change with 
time. The clear water analysis shows that this assumption 
is incorrect. 
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Fig. 11. The band 3 global clear water radiances are plotted over the mission. The radiances are mean values 
computed for all pixels where the water depth is greater than 1 km and the chlorophyll concentration is less 
than 0.15 mg m~ 3 . 



Fig. 12. The band 5 global clear water radiances are plotted over the mission. The radiances are mean values 
computed for all pixels where the water depth is greater than 1 km and the chlorophyll concentration is less 
than 0.15 mg m~ 3 . 


Recently, a comparison of the first 50 SeaWiFS lunar 
measurements with the United States Geological Survey 
(USGS) Robotic Lunar Observatory (ROLO) lunar irradi- 
ance model (Kieffer et al. 2002) verified that bands 3 and 
4 have small drifts with time. The comparisons show that, 
on average, the radiometric responses of bands 3 and 5 
have been decreasing at a rate of 0.35% per 1,000 days, or 
0.60% over the course of the mission. The effect of this 
time drift on the Lwn values would be an order of magni- 
tude larger. To correct the residual drift in the SeaWiFS 
time correction, a time correction of —3.5 x 10 -6 per day 
was applied to the mean of bands 3 and 4 before the mean 
was used to normalize the lunar time series. The overall 
time corrections for bands 1-8 were recomputed and the 


vicarious calibration of SeaWiFS were run again. 

The fourth reprocessing was run with the updated time 
correction and the 8-day composite clear water analysis 
was rerun. This analysis showed no significant changes in 
bands 1-6 over the course of the mission. Figures 11 and 
12 show the updated clear water time series for bands 3 
and 5 as solid lines. Band 3 exhibited a change of +0.16% 
over 1,000 days and +0.27% over the course of the mission 
(a rate of change of +1.5 x 10 -6 per day). Band 5 exhibits 
a change of -0.38% over 1,000 days and -0.64% over the 
course of the mission (a rate of change of —3.5 x 10 -6 per 
day). Any residual time drift in the instrument calibration 
was reduced by a factor of 8-10. 
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Fig. 13. The mirror side A reflectances for band 1 are plotted over the mission. The lines are the three 
piecewise linear fits to the data. 


2.4 MIRROR SIDE CORRECTIONS 

The SeaWiFS level -lb calibration equation includes a 
correction for the different reflectances of the two sides 
of the instrument’s half-angle mirror (Eplee and Barnes 
2000). Over the course of the mission, the SeaWiFS mir- 
ror side reflectances have changed at the one-count level, 
giving rise to artifacts in the ocean data which are visible 
on close examination. Through the third reprocessing, the 
problem was addressed by updating the operational cali- 
bration table with step-function changes to the mirror side 
correction factors. For the fourth reprocessing, the mir- 
ror side reflectance time series was fit with a set of three 
piecewise linear corrections. An example reflectance time 
series is shown for band 1, mirror side A in Fig. 13. The 
piecewise linear corrections were incorporated into the re- 
vised calibration table used in the fourth reprocessing, thus 
minimizing the future updates required for the operational 
calibration table. 


2.5 SUMMARY 

Several changes were implemented in the SeaWiFS on- 
orbit calibration methodology. Periodicities in the lunar 
calibration time series were used to revise the detector 
temperature corrections, thus allowing the revised time 
series to be fit by a single decaying exponential function 
of time. Comparisons were made of the SeaWiFS lunar 
measurements to the USGS ROLO lunar irradiance model 
to develop a time correction for the lunar data normal- 
ization. Time-dependent mirror side reflectance correc- 
tions were also implemented. These changes in calibration 
methodology resulted in a more robust on-orbit calibra- 
tion of SeaWiFS for the fourth reprocessing of the mission 
data set. An additional result of these changes is that 
fewer updates to the operational calibration table will be 
required to maintain a stable calibration of the SeaWiFS 
ocean data. 
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Chapter 3 


Changes to the Vicarious Calibration of SeaWiFS 

Robert E. Eplee, Jr., and Robert A. Barnes 
Science Applications International Corporation 
Beltsville, Maryland 

Sean W. Bailey 

Futuretech Corporation 
Greenbelt , Maryland 

P. Jeremy Werdell 

Science Systems and Applications , Incorporated 
Lanham, Maryland 

Abstract 

The SeaWiFS Project vicariously calibrates the visible bands of the instrument against in-water measurements 
from MOBY in order to achieve the necessary calibration accuracy for the SeaWiFS retrievals of water-leaving 
radiances. The SeaWiFS CVT implemented several changes in the vicarious calibration procedure for the fourth 
reprocessing of the SeaWiFS mission data set. The most significant change is the use of stray light corrected 
MOBY data in the vicarious calibration. More stringent data quality screening criteria were also implemented 
for both MOBY and SeaWiFS data to determine the vicarious calibration matchups. Finally, an inverse vicarious 
calibration procedure was implemented in which the vicarious gains are computed at the top of the atmosphere. 
The SeaWiFS TOA radiances are compared with MOBY water-leaving radiances that have been propagated to 
the top of the atmosphere using the SeaWiFS-retrieved atmospheric correction parameters. These changes to 
the vicarious calibration procedure enabled the production of a more robust set of vicarious gains for the fourth 
reprocessing of the SeaWiFS mission data set. 


3.1 INTRODUCTION 

The SeaWiFS CVT implemented several changes in the 
SeaWiFS vicarious calibration procedure for the fourth re- 
processing. The in situ measurements are provided by 
stray light corrected MOBY data. The SeaWiFS data is 
composed of 101x101 pixel extracted scenes (subscenes) 
centered on MOBY. More stringent data quality screening 
criteria were imposed on both the MOBY and SeaWiFS 
data sets. Vicarious gains were generated using the in- 
verse calibration procedure. These changes are detailed in 
the following sections of this chapter. 

3.2 STRAY LIGHT CORRECTION 

In vicariously calibrating SeaWiFS, the assumption is 
made that the uncertainties in the MOBY radiances are 
small (Eplee et al. 2001). Analyses of the MOBY data 
performed by the MOBY Project has shown that stray 
light within the instrument is biasing the MOBY measure- 
ments in the blue end of the visible spectrum (Clark et al. 
2001). As part of an ongoing intra-agency collaboration, 


the MOBY Project is working with NIST to characterize 
the stray light within the MOBY spectrometers and to 
develop stray light corrections for the MOBY data. The 
MOBY Project plans to reprocess the entire MOBY time 
series, starting with the more recent data. The full set of 
stray light corrected data over the time range of Novem- 
ber 1999 to March 2002 were implemented, along with 
targeted match-up data from September 1997-November 
1999, into the vicarious calibration for the fourth repro- 
cessing. The revised vicarious gains provide water-leaving 
radiances (Lw) that are higher for bands 1 through 5 
(555 nm), compared with the gains for the third reprocess- 
ing. 

3.3 DATA SCREENING PROCEDURES 

Before matchups between MOBY and SeaWiFS can 
be used in vicarious calibration, both data sets must be 
screened for data quality. The data screening procedures 
for both data sets have been enhanced for the fourth re- 
processing. 
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MOBY-derived values of Lw are computed by prop- 
agating the upwelling radiance (L u ) measurements made 
at depth to the surface using the diffuse attenuation co- 
efficient (Clark et al. 2001). Measurements of Lw com- 
puted from the L u measurement at the top arm of MOBY 
(1.5m depth) were compared to those computed from the 
L u measurement at the middle arm (5 m depth). The ma- 
jority of these measurement pairs are in good agreement; 
however, divergence between the two water-leaving radi- 
ances on the order of 5-30% can occur. These are likely 
the result of unfavorable sky and or sea conditions. Con- 
sequently, the following exclusion criteria for MOBY data 
were implemented: 

1) Data from both the top and middle arms must ex- 
ist;, and 

2) The combined uncertainty in the Lw measurements 
for all bands from both arms must be less than 10%, 
as determined by the square root of the sum of the 
squared absolute percent differences for bands 1-5, 
corrected for the cosine of the solar zenith angle. 

The currently available set of stray light corrected 
MOBY data has yielded 163 points where the skies are 
clear over MOBY at the time of the SeaWiFS overpass 
(Fig. 14). The exclusion criteria imposed on the MOBY 
data have reduced the number of matchups for the vicar- 
ious calibration to 45. Further data will be incorporated 
into future vicarious calibrations of SeaWiFS as additional 
stray light corrected MOBY data become available and as 
future refinements of the MOBY stray light correction are 
implemented. In all probability, these additions to the 
match-up data set will not significantly alter the results 
obtained for the fourth reprocessing. 

The SeaWiFS match-up data consist of 101x101 pixel 
subscenes centered on the MOBY location. The size of 
these extracts ensures that stray light within SeaWiFS will 
be handled properly by the level -2 conversion software in 
the vicinity of MOBY. The analysis for each matchup is 
performed on 5x5 pixel subscenes centered on the MOBY 
location. Each subscene is checked for data quality, screen- 
ing out matchups for clouds and cloud shadows, stray light, 
sun glint, high satellite and solar zenith angles, high L t in 
bands 7 and 8, and aerosol optical depths greater than 
0.1. All 25 pixels in a subscene are required to be clear 
for the matchup to be used in the vicarious calibration. 
These quality constraints imposed on the SeaWiFS data 
have reduced the number of matchups to 23. The vicari- 
ous calibration used in the fourth reprocessing is based on 
23 SeaWiFS-MOBY matchups selected from the original 
pool of 163 potential matchups. 

3.4 CALIBRATION METHODOLOGY 

In order to reduce the uncertainty in the vicarious gains, 
the vicarious calibration procedures for the fourth repro- 
cessing were modified. The major change was switching 
from the forward vicarious calibration procedure, which 


compares the MOBY and SeaWiFS data at the sea sur- 
face, to an inverse calibration procedure, which makes the 
comparison at the top of the atmosphere. The calibration 
methodologies discussed here apply to the SeaWiFS visi- 
ble bands, bands 1-6. The vicarious calibration of band 7, 
was performed for the fourth reprocessing to optimize the 
retrieval of the aerosol optical properties, as described in 
Eplee et al. (2001). 

The forward vicarious calibration procedure compares 
normalized water-leaving radiances ( Lwn ) between Sea- 
WiFS and MOBY to determine the vicarious gains for Sea- 
WiFS bands 1-6. For a given matchup, the Lwn values for 
SeaWiFS, L(\s), is averaged over the 5x5 pixel subscene 
and this average value L(Xs) is divided by the MOBY ra- 
diance L(Xm) for that matchup. The mean of this ratio, 
L(X S /m)> is computed over all of the individual matchups. 
The vicarious gain for each band is adjusted iteratively 
over multiple runs of level -2 conversion software on the 
match-up data set until the mean L(X S /m) value converges 
to unity. The forward vicarious calibration procedure al- 
lows the effects of iterative NIR correction algorithms to be 
incorporated into the vicarious gains. The vicarious gains 
were derived for the third and previous reprocessings us- 
ing the forward vicarious calibration procedure (Eplee et 
al. 2001). 

The inverse vicarious calibration procedure compares 
TOA radiances between SeaWiFS and MOBY to deter- 
mine the vicarious gains for SeaWiFS bands 1-6 (Fig. 15). 
For a given matchup, the atmospheric correction param- 
eters retrieved for the SeaWiFS data are used to normal- 
ize the water-leaving radiances measured by MOBY and 
to propagate the resulting normalized water-leaving radi- 
ances to the top of the atmosphere. For the matchup, the 
vicarious gains are computed from the Lt(Xs/M) values in 
each pixel, averaged over the 5x5 subscene. The overall 
vicarious gains are the means of the gains computed for 
the individual matchups. 

The inverse vicarious calibration procedure generates 
gains which agree with those produced by the forward cal- 
ibration to within 0.04% for the same calibration condi- 
tions. At the same time, the inverse calibration offers a 
number of advantages over the forward calibration. The 
inverse calibration runs more quickly than the forward cal- 
ibration, because the gains are computed directly from a 
single run of the level -2 conversion software. The inverse 
calibration allows the effect of an instrumental offset to 
be investigated, although no evidence of an offset has been 
observed to date. Additionally, the inverse calibration pre- 
serves spectral shape information in the matchups by com- 
puting a set of gains for each matchup. The forward cali- 
bration computes the gains from the mean of the match-up 
ratios, thus diluting the spectral shape information. The 
only drawback to the inverse calibration procedure is that 
it cannot incorporate the effects of iterative NIR correction 
algorithms. Consequently, the inverse calibration proce- 
dure can only be used when the NIR correction algorithms 
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Fig. 14. Current SeaWiFS-MOBY matchups are shown. The distribution of the matchups in time is primarily 
due to the availability of stray light corrected MOBY data. 


SeaWiFS Data: 


MOBY Data: 



Fig. 15. The diagram shows how SeaWiFS and MOBY data are used to compute vicarious gains for SeaWiFS 
from ratios of TOA radiances. 
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Fig. 16. The vicarious gains for band 1 are shown. The full match-up data set is represented by pluses, while 
the matchups used in the vicarious calibration are represented by filled circles. The horizontal line shows the 
vicarious gain for the band. 



Fig. 17. The vicarious gains for band 2. The symbols used are the same as in Fig. 16. 



Fig. 18. The vicarious gains for band 3. The symbols used are the same as in Fig. 16. 
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Fig. 19. The vicarious gains for band 4. The symbols used are the same as in Fig. 16. 



Fig. 20. The vicarious gains for band 5. The symbols used are the same as in Fig. 16. 



Fig. 21. The vicarious gains for band 6. The symbols used are the same as in Fig. 16. 
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Table 3. These vicarious gains were used for the fourth reprocessing. 


Band 

1 

2 

3 

4 

5 

6 

7 

8 

Gain 

1.013007 

0.996384 

0.962951 

0.982130 

0.991338 

0.956581 

0.9380 

1.0000 


have been disabled or have been enabled with a transi- 
tion (Chapt. 4 in this volume) at low chlorophyll concen- 
trations. The vicarious gains were derived for the fourth 
reprocessing using the inverse calibration procedure. 

3-5 VICARIOUS GAINS 

Throughout the analyses performed in preparation for 
the fourth reprocessing, the vicarious calibration was up- 
dated as specific changes were implemented in the data 
processing methodology that would effect the derivation 
of the vicarious gain. Such changes included the MOBY 
stray light correction, the switch to the revised calibration 
table, the implementation of the ramp-in of the NIR cor- 
rection algorithm at low chlorophyll levels, and the Fresnel 
transmittance correction at the air-sea interface. 

The vicarious calibration data set used for the fourth 
reprocessing is composed of 23 SeaWiFS-MOBY matchups 
selected from an original pool of 163 potential matchups. 
Time series of the TOA gains for bands 1-6 are shown 
in Figs. 16-21. The 163 potential matchups in each band 
are shown by pluses and the 23 matchups used in comput- 
ing the overall vicarious gains are shown by filled circles. 
The data quality screening applied to both the MOBY 
and SeaWiFS data has significantly reduced the scatter in 
the individual gains used to compute the overall vicarious 
gains. The vicarious gains used for the fourth reprocessing 
are shown in Table 3 and are plotted as horizontal lines in 
Figs. 16-21. It should be noted that the computed gains 
pass through the clusters of matchups that were excluded 
from the vicarious calibration. 

As was discussed above, the gain for band 7 used in the 
fourth reprocessing was computed to optimize the retrieval 
of the aerosol optical properties, as described in Eplee et 


al. (2001). 

As a validation of the vicarious gains, the Lwn values 
retrieved by SeaWiFS for the 23 matchups were compared 
with the water-leaving radiances measured by MOBY and 
were then normalized using the SeaWiFS-derived atmos- 
pheric correction parameters. The RMS differences be- 
tween the SeaWiFS retrievals and the MOBY measure- 
ments are shown as percentages for bands 1-5 in Table 4. 
The radiances in band 6 are close to zero, so the RMS dif- 
ference for band 6 is not a meaningful number. For bands 
1-5, the differences increase with decreasing wavelength, 
which is to be expected if the major source of the uncer- 
tainties in the vicarious gains is the extrapolation of the 
atmospheric correction from the NIR to the visible. 


Table 4. SeaWiFS and MOBY Lwn agree to bet- 
ter than 10%. 


Band 

i 

2 

3 

4 

5 

RMS Differences [%] 

9.74 

7.86 

4.46 

3.69 

2.41 


3.6 SUMMARY 

The CVT has improved the quality of both the MOBY 
and SeaWiFS data sets which are used in the vicarious cal- 
ibration for the fourth reprocessing. The stray light correc- 
tion of the MOBY data has improved the spectral shape 
of the vicarious gains. The stricter data quality screen- 
ing imposed on both the MOBY and SeaWiFS data has 
reduced the uncertainty in the gains. In addition, an in- 
verse calibration procedure was implemented to determine 
the vicarious gains. The result of these changes in the vi- 
carious calibration methodology is a more robust vicarious 
calibration of SeaWiFS for the fourth reprocessing. 
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Modifications to the SeaWiFS NIR Correction 

Sean W. Bailey 
Futuretech Corporation 
Greenbelt , Maryland 

Wayne D. Robinson and Bryan A. Franz 
Science Applications International Corporation 
Beltsville, Maryland 

Abstract 

This chapter describes several changes made to the NIR correction implemented during the fourth reprocessing 
of SeaWiFS data. The changes made were both to correct flaws in the actual implementation of the NIR 
correction and to improve the estimation of NIR reflectance. The changes include revised water absorption 
coefficients, the addition of a scaling factor to allow for the gradual introduction of the NIR correction as the 
derived chlorophyll concentration increases, improvements to the iteration control, and an alternative model 
for the backscatter estimate. The revised NIR correction reduces discontinuities in the aerosol model selection 
at the boundaries where the correction is introduced. It also reduces the attenuation effect that the original 
correction had on high chlorophyll concentrations. 


4.1 INTRODUCTION 

One of the primary assumptions used in the SeaWiFS 
atmospheric correction process is what is known as the 
black pixel assumption. This assumption states that the 
contribution to the total radiance by water-leaving radi- 
ance in the NIR region is negligible, and therefore, the 
TOA radiance in the NIR bands is solely due to atmos- 
pheric path radiance. For most conditions, this assump- 
tion is valid, however, in turbid or highly productive wa- 
ters, this assumption breaks down. Various methods have 
been proposed to account for NIR water-leaving radiance 
for these situations. The SeaWiFS Project implemented 
an NIR correction (Siegel et al. 2000) with the third re- 
processing in May 2000. Since then, refinements to the 
implementation of the NIR correction algorithm, as well as 
significant changes to the algorithm have been developed. 
While the basic concept for estimating NIR water-leaving 
radiance proposed by Siegel et al. (2000) is maintained, 
the methods have changed. This chapter describes these 
modifications. 

4.2 BACKSCATTER MODEL 

The NIR correction implemented with the third repro- 
cessing used the particulate backscatter (£>b P ) estimate de- 
scribed in Loisel and Morel (1998) to estimate backscatter- 
ing (bb) at 555 nm and the spectral backscattering function 


described in Morel (1988) to extrapolate this value into the 
NIR. In practical application, the NIR correction is applied 
most often in highly productive waters (>3mgm” 3 C a ) 
and turbid or Case-2 waters. The Loisel and Morel model 
is explicitly Case-1, developed with a data set of chloro- 
phyll values ranging from 0.01-4.53 mg m -3 , with the ma- 
jority of the data having values of less than 1.5 mg m -3 . 
The Loisel and Morel model is not applicable to Case-2 
waters. Gould et al. (1999) developed a spectral depen- 
dence model for the scattering coefficient in Case-1 and 
Case-2 waters. This model was implemented in place of 
the Morel (1988) model as being more appropriate to the 
practical application of the NIR correction in SeaWiFS. 

A reflectance-based estimate of backscatter (Sydor and 
Arnone 1997) at a reference wavelength of 670 nm has 
been adopted in place of the chlorophyll-based estimate 
described in Loisel and Morel (1998). This reflectance- 
based estimate relies on the assumption that water ab- 
sorption is dominate at wavelengths greater than 650 nm. 
This assumption allows for the estimation of fy>(670) from 
R rs (670) by (3): 

6 6 (670) = *rs(670)^®, (3) 

where a w is the absorption coefficient for water. 

There can, however, be a measureable amount of par- 
ticulate absorption (a p ) at 670 nm. Two additional absorp- 
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Fig. 22. Water absorption spectra with SeaWiFS spectral response overlayed. The open circles indicate 
values for water absorption cited in Siegel et al. (2000), and the solid diamonds indicate revised values. 


tion terms have, therefore, been added, particulate absorp- 
tion (Bricaud et al. 1998), 

a p (670) = 0.019890 C^* 818 , (4) 


and detrital-gelbstoff absorption (a d g) (Chapt. 9 in this 
volume), 


a dg (670) = 0.15 - 0.19 


i? rs (555) - i? rs (670) 
i? rs (555) 


The total absorption (a to t) term becomes, 


(5) 


a tot (670) = o„,(670) + a p (670) + a dg (670), (6) 


that the retrieval of aerosol properties can be performed 
correctly. Because these aerosol properties should be inde- 
pendent of in- water constituents, it is assumed that the 
aerosol products (e.g., Angstrom) should not be corre- 
lated with chlorophyll retrievals or water-leaving radiance 
(or reflectance). Visual inspection of LAC resolution im- 
agery processed with the chlorophyll-based NIR correction 
showed patterns in the Angstrom product that mirrored 
the chlorophyll product, with distinct boundaries at the 
transition regions of the image where the NIR correction 
was applied. These same scenes were processed with the 
reflectance-based NIR corrected; the Angstrom product 
showed a correlation with the in-water constituents only 
under extremely turbid conditions. 


resulting in a new 6^(670) equation, 

& 6 (670) = fir 8 (670) at °^ 0) . 


(7) 


Because there are no reliable in situ measurements of 
NIR water-leaving radiance, validating the NIR correction 
is difficult. In order to evaluate whether or not the changes 
to the NIR correction algorithm constituted an improve- 
ment over the correction implemented with the third re- 
processing, the Angstrom exponent product was used as 
a proxy. The goal of the NIR correction is to account 
for water-leaving radiance in the NIR bands under condi- 
tions for which the black-pixel assumption is not valid, so 


4.3 ABSORPTION COEFFICIENTS 

The NIR correction algorithm uses the absorption coef- 
ficient of water in its formulation. The values used with the 
Siegel et al. (2000) algorithm, implemented with the third 
reprocessing, were derived from published values [Hale and 
Query (1973) and Smith and Baker (1981)] for the nom- 
inal center wavelengths of SeaWiFS NIR bands. Given 
the strong absorption of water in the NIR and the 40 nm 
bandpasses for the SeaWiFS NIR bands (Fig. 22), these 
numbers are underestimated. A revised set of bandpass- 
averaged numbers were generated (Table 5), based on the 
newer published values for water absorption than those 


27 



Algorithm Updates for the Fourth SeaWiFS Data Reprocessing 


cited in Siegel et al. (2000), Kou et al. (1993), and Pope 
and Fry (1997). 


Table 5. Revised absorption coefficient for water 
based on full SeaWiFS bandpasses. 



Band 6 
670 nm 

Band 7 
765 nm 

Band 8 
865 nm 

Q "U)\ 

0.4346 

2.5500 

4.2860 

Revised a w 

0.4458 

2.9530 

4.8680 

Percent difference 

2.5770 

15.8039 

13.5791 


f Siegel et al. 2000. 


4.4 SCALING FACTOR 

A change was made in the NIR correction to introduce 
a scaling factor to the NIR radiances as a function of the 
amount of chlorophyll. The NIR radiances are modified 
by a factor that is zero for chlorophyll from 0-0.7 mgm” 3 , 
1.0 for chlorophyll of 1.3mgm” 3 and above, and linearly 
increasing from 0-1.0 in the 0.7-1. 3 mgm -3 chlorophyll 
range. The change was made for several reasons. 

Determining the NIR radiances is a computationally 
intensive iterative process that only makes a significant 
correction at chlorophyll values above 2 mgm -3 (Siegel et 
al. 2000). Below this value, the assumption of no NIR 
radiances is reasonable. Performing the iterations above 
chlorophyll concentrations of 0.7 mgm” 3 increases the pro- 
cessing speed and still addresses regions where the non-zero 
NIR is important. 

Another reason for not having the NIR correction ap- 
plied at low chlorophyll concentrations is that the vicarious 
calibration process can be performed much more quickly 
if the iteration for the NIR does not have to be done. As 
the chlorophyll at the vicarious calibration site rarely gets 
higher than 0.1 mgm” 3 , the threshold of 0.7 mgm” 3 allows 
for the vicarious calibration to be performed without the 
NIR iteration. In addition, users of the chlorophyll data 
had concerns during the third reprocessing that the NIR 
algorithm was changing the low chlorophyll values. This 
occurred because the initial first guess NIR was computed 
for a chlorophyll value of 0.3 mgm” 3 , which made a small, 
but noticable difference in the chlorophyll retrievals made 
with and without the NIR algorithm enabled. The current 
use of zero NIR removes this problem. 

The scaling factor for the transition from no NIR at 
a chlorophyll concentration of 0.7 mgm” 3 , to the full NIR 
at 1.3 mgm” 3 linearly increases from 0-1.0 to remove the 
occurrence of any artifacts in the chlorophyll distribution 
that would be caused by the abrupt change of the NIR ra- 


diances at one chlorophyll value. The transition is instead 
distributed over a 0.6 mgm” 3 range, where the NIR correc- 
tion has a small effect on the chlorophyll. The transition 
starts well above the typical conditions used for vicarious 
calibration and well below where the effect of the NIR cor- 
rection is significant. 

4.5 ITERATION CONTROL 

In order to improve the estimate of the NIR reflectance, 
the correction algorithm was implemented in an iterative 
manner. The original iteration control was quite simplis- 
tic. The algorithm would iterate until retrieved chlorophyll 
concentration changed by less than 20% from the previous 
iteration, for a maximum of 10 iterations. If the chloro- 
phyll retrieval failed during the iteration process, the it- 
eration was reinitialized once with a chlorophyll value of 
5 mgm” 3 . The iteration control was modified to improve 
the behavior. The following changes were implemented: 

1. Dampening between iterations: The NIR reflectance 
returned is averaged with the previous iteration re- 
sult to minimize the possibility of large fluctuations 
between iterations. 

2. Iterations stop when the average NIR water-leaving 
reflectance changes by less than 2%, or Rj. s (670) < 
0 . 0 . 

3. The initial C a value is changed from 0.3 mgm” 3 , to 
0.0 mgm” 3 . The initial i2r S (670) value is set to 0.0. 

4. The iteration is reinitialized on any iteration where 
the chlorophyll retrieval fails, up to the maximum 
allowed iterations (currently 10), using: 

a) C a = i5.0, where i is the number of iterations, 
and 

b) J? rs (670) = 5.0(0.00032 + 0.00021C o ). 

4.6 DISCUSSION 

Changes were made to the methods used by the Sea- 
WiFS project to estimate water-leaving radiance in the 
NIR bands to allow for atmospherically correcting data 
where the black pixel assumption is not valid. These chan- 
ges included a scaling factor, a revision of the iteration con- 
trol for the algorithm, an update to the absorption coeffi- 
cients for water used in the algorithm, and the switch to a 
reflectance-based estimate of backscatter with a new func- 
tion to estimate NIR backscatter from a reference wave- 
length in the visible spectrum. These changes produce a 
marked improvement to the NIR correction as it is applied 
operationally to SeaWiFS data. 
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Changes to the Atmospheric Correction Algorithm 
and Retrieval of Oceanic Optical Properties 
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Greenbelt, Maryland 

Menghua Wang 

University of Maryland, Baltimore County 
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Abstract 

In preparation for the fourth SeaWiFS reprocessing, a series of algorithm changes were implemented to enhance 
the performance of the atmospheric correction process and to improve the quality and consistency of oceanic 
optical property retrievals. Included in these changes was the introduction of a filtering process to reduce the 
relative noise between the two NIR channels. In addition, several modifications were made to improve the 
handling of the SeaWiFS out-of-band response, and the normalization of water-leaving radiances was extended 
to account for Fresnel transmittance effects through the air-sea interface. These and other algorithm updates 
are described within this chapter. 


5.1 INTRODUCTION 

While the general approach to SeaWiFS atmospheric 
correction over oceans did not change for the fourth repro- 
cessing, a number of refinements were implemented and 
evaluated. Several of these modifications were found to 
yield significant improvement in the quality and consis- 
tency of oceanic optical property retrievals; the changes 
were included in the final reprocessing software and are 
discussed in detail in this chapter: 

a) A filtering scheme for reducing aerosol model selec- 
tion noise, 

b) A modification to improve algorithm performance 
in very clear atmospheres, 

c) Updates to the corrections for out-of-band response, 

d) An extension of the water-leaving radiance normal- 
ization to account for Fresnel transmittance through 
the air-sea interface, and 

e) A fix for aerosol model ambiguity problems. 

5.2 RELATIVE NOISE REDUCTION 

The SeaWiFS atmospheric correction algorithm (Gor- 
don and Wang 1994) relies on the single-scattering aerosol 


reflectance ratio (e) between the two NIR bands at 765 and 
865 nm to select the aerosol type. The atmospheric correc- 
tion is, therefore, highly sensitive to any relative noise be- 
tween these two NIR channels. A filtering technique was 
developed to reduce the relative noise in the NIR band 
ratio, which thereby reduces the small scale variability in 
aerosol model selection. The smoothing filter adjusts the 
radiance in the 765 nm channel to minimize local variabil- 
ity in the observed NIR aerosol ratio (i.e., the multiscat- 
tering equivalent of atmospheric e, i.e., € ms ), while leaving 
the 865 nm radiance (which governs aerosol concentration) 
unchanged. The effect of this smoothing is to reduce pixel- 
to-pixel variability in the retrieved aerosol type, which ul- 
timately reduces atmospheric correction noise in the re- 
trieved water-leaving radiances. 

The effect of this filtering can easily be seen in level -2 
images of e, the Angstrom coefficient, and, to a lesser ex- 
tent, aerosol optical thickness at 865 nm. The value of 
this smoothing will diminish with increasing spatial and 
temporal averaging, and is more readily seen as reduced 
speckling in level -2 oceanic and atmospheric optical prop- 
erty retrievals. The smoothing was found to induce no 
bias-change in either the aerosol optical thickness or the 
water-leaving radiances. The smoothing algorithm follows: 
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L Define NIR aerosol radiance at pixel i for wavelength A 
as 




1 

^ox (A, i) 


^oz(A } i) 


Z/ r (A, i) 


• (8) 


2. Given a scan-pixel window centered on pixel i.e., i C) 
containing a total of n unmasked pixels, compute mean 
L a ( A) at i c as 


(L«(A,i c )) = iVL a (A,i), (9) 

n z ' 

for i — 1, n and A = 765 or 865 nm. 

3. Compute the mean multiscattering e value at pixel i c as 


(1,(765,1,)) 
(L a (865, i c )) ' 


(10) 


4. Then compute a new L tt (765,i c ) which would yield the 
mean e value when combined with the original L a (865^i c ) 
as 

K(765J c ) = e ms T a (865,i c ). (11) 

5. Reconstruct the TOA radiance at 765 nm as 


L t (765, i c ) = t 0 z(765, i c ) \L' a (765, i c )W(76 5, i c ) 
+ L r (765, z c )l + ££/(765, z c ), 


( 12 ) 


where 


■ L t (A,i) is the observed TOA radiance for wave- 
length A at pixel i; 

■ tLf(\,i) is the white cap radiance, transmitted to 
the TOA; 

■ T r (A, i) is the Rayleigh path radiance; 

■ L 0 (A,i) and L' a (A,i) are the aerosol path radiance 
(including Rayleigh aerosol interaction) and the new 
computed aerosol path radiance, respectively; 

■ *ox(A,i) is the oxygen transmittance; 

■ ^oz(A,i) is the ozone transmittance; and 

■ Mean values are within angle brackets, ( ). 

Using the filter- adjusted TOA radiance at 765 nm, the Sea- 
WiFS atmospheric correction algorithm is then operated in 
the standard manner. 

It is desirable to keep the filter window size as small as 
possible, to limit the reduction of real changes in aerosol 
type. The window size, however, needs to be large enough 
to allow sufficient sample size for the averaging to be ef- 
fective. In addition to varying the size, it is also possible 
to change the shape. This can be achieved by introduc- 
ing the concept of a filter window kernel, which indicates 
which pixels within the window will be considered in com- 
puting the filtered value. Consider these two examples of 
a 5 x 5 filtering window, where the value of 1 indicates that 
the pixel at that location will contribute. 


Square 5x5 
11111 
11111 
11111 
11111 
11111 


Diamond 5x5 
0 0 10 0 
0 1110 
11111 
0 1110 
0 0 10 0 


For the same window size, the diamond filter kernel 
reduces the number of contributing samples by approxi- 
mately 50% over the square kernel, and the radius of influ- 
ence is never greater than two pixels. While this reduces 
the number of samples contributing to the mean, the di- 
amond shape is better designed to minimize line-by-line 
digitization problems such as those associated with Sea- 
WiFS mirror side differences. This is because the diamond 
kernel gives nearly equal weight to the odd and even lines, 
while the square kernel yields a 3-to-2 over weighting of 
opposing lines. 

For SeaWiFS LAC resolution data, it was found that 
the NIR relative noise reduction filter with a 5x5 diamond 
kernel gave the best compromise between noise reduction 
and aerosol smoothing. Figure 23 shows a LAC subscene 
of Lwn values at 443 nm, before and after smoothing. 

Unfortunately, it was found that the filtering approach 
did not always reduce noise in SeaWiFS GAC resolution 
scenes. The problem appears to be that the GAC data 
set, being subsampled at the sensor, does not contain a 
complete record of the bright sources observed by the in- 
strument. This is a fundamental limitation of SeaWiFS, 
as it is simply not possible to identify and correct for all 
stray light contamination in the GAC data set. Any al- 
gorithm that combines neighboring pixels will, therefore, 
increase the probability of stray light contamination in a 
given pixel. The problem is most significant in the vicin- 
ity of scattered clouds. Considering this limitation, and 
the fact that GAC data are primarily used for generating 
spatial and temporal composites (where noise will be sig- 
nificantly reduced through averaging) , it was decided that 
the NIR relative noise reduction filter would not be applied 
to the GAC products in this reprocessing. 

5*3 CLEAR CONDITIONS 

Under very clear atmospheric conditions, the Rayleigh- 
subtracted radiance in the NIR approaches zero. When 
other uncertainties are included, the retrieved aerosol path 
radiances in the NIR may even go slightly negative. The 
aerosol model selection, therefore, becomes highly uncer- 
tain in clear atmospheres. As a result, the SeaWiFS at- 
mospheric correction algorithm often fails to obtain ocean 
color retrievals in the best of atmospheric conditions. A 
simple solution to this problem is to fix the aerosol type 
when the aerosol path radiance in one or both of the NIR 
bands approaches zero, and limit the aerosol radiance at 
865 nm to be greater than, or equal to, zero. With these 
two changes, it is possible for the atmospheric correction 
algorithm to procede when the retrieved aerosol concen- 
tration is effectively zero. 

The aerosol type to which low aerosol pixels will be 
fixed is a simple white aerosol, i.e., p a (A) = p a (865), where 
p a (A) is aerosol path reflectance at wavelength A. The 
threshold below which the aerosol model will be fixed has 
been set at a very conservative value of p a (NIR) = 0.0001. 
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Fig. 23. Sample LAC image of normalized water-leaving radiance at 443 nm, before and after application of 
the NIR relative noise reduction filter: a) the scene without filtering, and b) the same scene with filtering 
applied. 


For the 765 and 865 nm channels, this reflectance value cor- 
responds to slightly greater than 1 digital count. Raising 
this value will force a larger percentage of the pixels to 
assume white aerosols, thus bypassing the aerosol model 
selection process of Gordon and Wang (1994). 

In conjunction with the above enhancement, the occur- 
rence of [( L t — tLf)/t oz — L r \ < 0 in one or more of bands 
2-8 will no longer be considered an atmospheric correction 
failure condition. 

5.4 OUT-OF-BAND CORRECTION 

The SeaWiFS spectral bands cover the range from 380- 
1,150 nm, with nominal band centers at 412, 443, 490, 510, 
555, 670, 765, and 865 nm. The spectral bandwidth, which 
is defined as the full width at half the maximum response, 
is 20 nm for the first six bands and 40 nm for the two NIR 
bands. The SeaWiFS bands, however, are known to ex- 
hibit significant response well beyond the quoted spectral 
range of the bandpasses. Throughout the SeaWiFS at- 
mospheric correction process, adjustments are made to ac- 
count for this out-of-band response (Gordon 1995). In the 
third SeaWiFS reprocessing, additional corrections were 
added to adjust the retrieved Lwn values to correspond to 
the nominal band center wavelengths (Wang et al. 2001). 
For the fourth SeaWiFS reprocessing, several modifications 


were made to the out-of-band corrections for the water- 
leaving radiances and derived reflectances. These changes 
are discussed in the sections that follow. 

5.4.1 Remote Sensing Reflectance 

The solar irradiance ( Fo ) values used in the SeaWiFS 
atmospheric correction processing are band-averaged quan- 
tities. This means that the solar spectrum has been con- 
volved with the relative spectral response (RSR) function, 
where the RSR may include significant out-of-band re- 
sponse. In the third SeaWiFS reprocessing, an algorithm 
was introduced to correct the Lwn retrievals from band- 
averaged quantities to a nominal wavelength. Unfortu- 
nately, the out-of-band corrected Lwn vlaues were still 
normalized by the band-averaged Fo when computing i? ra . 
The R ts ratios between the visible bands were then used 
as input to the OC4 chlorophyll algorithm (O’Reilly et 
al. 2000), so the resulting chlorophyll retrievals may have 
been slightly biased. In the fourth reprocessing, the nom- 
inal band Lwn values were normalized by nominal band 
Fo values when computing R rs . 

5.4.2 Water-Leaving Radiance 

The Lwn retrievals are computed as band-averaged 
values, with an out-of-band correction applied prior to use 
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Fig. 24. Comparison of Morel and Maritorena (2001) bio-optical model with Gordon et al. 1988 model. 



Fig. 25. Mean along-scan normalized water-leaving radiance retrievals, before and after application of the 
Fresnel transmittance correction. The solid line is with correction. 
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in downstream computations such as chlorophyll retrieval. 
The correction applied in the third SeaWiFS reprocessing 
was computed from the Lwn (A 3 / 5 ) ratio, and was based 
on a chlorophyll-dominated Lwn spectrum which uses the 
Gordon et al. (1988) model. A revised set of correction fac- 
tors were generated for this reprocessing using the recently 
published clear- water reflectance model of Morel and Mar- 
itorena (2001). Figure 24 shows the comparison between 
the correction factors from Gordon et al. (1988) and Morel 
and Maritorena (2001), with the Morel and Maritorena 
(2001) model results indicated by the solid lines. 

This change to the out-of-band correction factors re- 
sults in a lowering of chlorophyll values in waters with a 
A3/5 Lwn ratio °f greater than approximately 2. For lower 
ratios, chlorophyll increases slightly. Differencing tests on 
global binned products indicate that the revised out-of- 
band correction results in a net decrease of 2. 5-4.0% in 
global averaged chlorophyll retrievals. 

5.5 FRESNEL TRANSMITTANCE 

The normalization of water-leaving radiance was modi- 
fied to include a correction for Fresnel transmittance 
through the water-atmosphere interface. The Gordon and 
Wang (1994) atmospheric correction algorithm assumes 
that the water-leaving radiance just beneath the ocean 
surface, Lwn( 0"), is uniform. For a flat ocean surface, 
the normalized water-leaving radiance just above the sur- 
face, Lwn{ 0 + ), can be related to the value just beneath 
the surface as 

Lwn{Q) = -4^W(0'), (13) 

n w 

where 0 , tf(0 ), and n w are the sensor zenith angle, the 
Fresnel transmittance of the air-sea interface, and the re- 
fractive index of the water, respectively. It is assumed 
in the above equation that the normalized water-leaving 
radiance just beneath the surface, Lwn( 0“), is uniform 
(independent of the sensor zenith angle 6). Without cor- 
rection, therefore, the SeaWiFS-derived normalized water- 
leaving radiance, Lwn , depends on the sensor zenith angle 
according to tf(0). In fact, the Fresnel transmittance ef- 
fect is part of the ocean bidirectional reflectance factors 
(f/Q correction) discussed by Morel and Mueller (2002). 


A simple correction was implemented to remove the air -sea 
transmittance effect on the SeaWiFS-derived normalized 
water-leaving radiances. The corrected values, L' WN , are 
computed as 

Lwn = ^Lwn(Q+), (14) 

where U WN and Lwn are normalized water-leaving radi- 
ances with and without surface transmittance correction, 
respectively. The U WN values (for all six visible bands) are 
now the SeaWiFS-derived normalized water-leaving radi- 
ances. Note that the correction does not affect the Sea- 
WiFS chlorophyll a concentration values because the trans- 
mittance effects are cancelled in the band ratio. 

The general effect of this change is to increase the 
normalized water-leaving radiance in all bands, with the 
largest increase occurring at the highest view zenith angle, 
reaching approximately 3% at the GAC limit of 56°. Fig- 
ure 25 shows the effect of this Fresnel correction across the 
full SeaWiFS scan. The scan trends were derived by simple 
averaging of water-leaving radiance retrievals within each 
scan pixel over a one year period in the relatively homoge- 
neous waters near Hawaii. The solid line is the corrected 
data, and the dashed line is the uncorrected data. Note 
that the roll off in radiance near the edge of the GAC 
swath and beyond is reduced. 

5.6 AEROSOL MODEL AMBIGUITY 

Under certain geometric conditions, the T99 (tropo- 
spheric, 99% relative humidity) and C50 (coastal, 50% 
relative humidity) aerosol models cross over in e space. 
These discontinuities appear along lines of constant scat- 
tering angle, and they are sometimes visible in images of 
water-leaving radiance and even chlorophyll. Application 
of the aforementioned NIR relative noise reduction filter- 
ing makes these effects even more apparent, as the aerosol 
model selection noise is reduced across the scattering angle 
isolines. A fix was developed which identifies these model 
cross-over conditions and revises the model selection re- 
sult accordingly. The details of these effects are discussed 
in Wang (2003). This is a relatively rare problem which 
will not significantly effect global results. 
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Abstract 

The flags and masks used for the SeaWiFS level-2 and leveT3 processing were updated for the recent fourth 
reprocessing. This chapter discusses the changes and why they were made. In many cases, underlying algorithms 
were changed. Some flags changed their states to either flagging (noting a condition), or masking (denoting 
data excluded from the product) to allow more data to be kept or to improve its quality. New flags were either 
introduced as a part of new algorithms or to denote the status of the data more clearly. The flag and mask 
changes significantly contributed to the improvement in the data quality and increased the amount of data 
retrieved. 


6.1 INTRODUCTION 

As part of the algorithm improvements for the fourth 
reprocessing, several improvements were made in the flag- 
ging and masking done in the level -2 and level -3 SeaWiFS 
data processing. Table 6 shows a list of the current flags 
used in the fourth reprocessing, their status in operational 
processing as a flag (to note a condition) or a mask (to 
note that the data are excluded from the product because 
of this condition), a brief description of the flag, and a 
general indication of how this flag changed from the third 
to the fourth reprocessing. Three general classes of change 
are noted. 

For many flags, the underlying algorithm that gener- 
ates the flag was changed. In all these cases, the new algo- 
rithm is an improvement over the old algorithm in that it 
characterizes the flag condition better, and in some cases, 
is able to allow for a greater number of good retrievals to 
be made. In Table 6, this is denoted as a change. 

Some flags had their status as a mask or a flag changed 
for the fourth reprocessing. Many times, this change was 
a direct result of the algorithm changes referred to above. 
Other changes were prompted by increased experience with 
the operation of the flag, which revealed the data quality 
or amount of retrieved data would be increased with the 


new masking. Table 6 notes cases in which data previously 
masked are now flagged, or vice versa. The changes in 
masking and flagging are noted for the level -2 and the 
level-3 product. For the level-3 product, there is no actual 
flag for the data. In the case where the change in flag status 
is denoted as flagged, it really means that the data which 
were previously excluded from the level -3 product are now 
included as a part of the product. 

Finally, four new flag conditions were added to the list 
of flags. These denote either new general conditions, or are 
part of new algorithms, or are used to more clearly separate 
conditions that occur during the level-2 processing. In 
Table 6, they are denoted as new flags. 

The following sections describe the flag and mask chan- 
ges in more detail. Section 6.2 describes an analysis to 
determine the effect of individual level -3 masks on the 
number of retrievals. Section 6.3 describes the changes 
in the flagging and masking for the level -2 product, Sect. 
6.4 looks at the changes in the masking for the leveT3 
product, and Sect. 6.5 summarizes the changes. 

6.2 FLAG EFFECT ANALYSIS 

Most of the flag and mask changes came naturally from 
the increased understanding of the SeaWiFS data. For 
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Table 6. Flags for the fourth reprocessing. The Mask In columns indicate that no geophysical data is created in 
the level -2 (L2) or level -3 (L3) data set if the flag conditions marked Y exist for that observation. The changes 
from the third reprocessing fall under one or more of the following categories: a) New (a new flag for the fourth 
reprocessing); b) Change (a direct change in the algorithm or threshold was made); and c) Mask in L2 or Mask 
in L3 indicates that in the fourth reprocessing, data are now excluded from the level-2 or level -3 product and 
Flag in L2 or Flag in L3 indicates that data which were previously exluded from the product are now included, 
(e.g., the HIGLINT flag is no longer a mask in L2 or L3, so the change reads: “Flag in L2 and L3”). 


No. 

Flag 

Name 

Mask In 

Description 

Change From 
Third Reprocessing 

L2 

L3 

1 

ATMFAIL 

Y 

Y 

Atmospheric algorithm failure 

— 

2 

LAND 

Y 

Y 

Land 

— 

3 

BADANC 



Missing ancillary data 

— 

4 

HIGLINT 



Sun glint contamination 

Flag in L2 and L3 

5 

HILT 

Y 

Y 

Total radiance above the knee 

Change 

6 

HISATZEN 


Y 

Satellite zenith angle above the limit 

Change 

7 

COASTZ 



Shallow water 

Change 

8 

NEGLW 



Negative water-leaving radiance in bands 7 and 8 

— 

9 

STRAYLIGHT 

Y 

Y 

Stray light contamination 

Change, mask in L2 

10 

CLDICE 

Y 

Y 

Clouds or ice 

Change 

11 

C0CC0LITH 


Y 

Coccolithophore bloom 

— 

12 

TURBIDW 



Turbid, Case-2 water 

Change 

13 

HISOLZEN 


Y 

Solar zenith angle above the limit 

— 

14 

HITAU 



High aerosol concentration 

— 

15 

LOWLW 


Y 

Low water-leaving radiance at 555 nm 

— 

16 

CHLFAIL 

t 

Y 

Chlorophyll not calculable 

Change 

17 

NAVWARN 


Y 

Questionable navigation (tilt change) 

— 

18 

ABSAER 


Y 

Absorbing aerosol index above the threshold 

Change 

19 

TRICHO 



Trichodesmium bloom condition 

— 

20 

MAXAERITER 


Y 

Maximum number of iterations in the NIR algorithm 

— 

21 

MODGLINT 



Glint corrected measurement 

— 

22 

CHLWARN 


Y 

Chlorophyll is out of range 

Mask in L3 

23 

ATMWARN 

! 

Y 

The e value is outside of reasonable range or Lw 

— 





at 490, 510, or 555 nm is less than zero 


24 

DARKPIXEL 



Rayleigh corrected radiance is less than zero for 

— 





any band 


25 

SEAICE 



Sea ice present based on climatology 

New 

26 

NAVFAIL 

Y 

Y 

Navigation of the line is bad 

New 

27 

FILTER 


Y 

Insufficient surrounding pixels for aerosol model filter^ 

New 

28-31 




Spare for future use 

N/A 

32 

OCEAN 



Ocean data 

New 


f The chlorophyll value is not computed, but first guess Lwn values are computed. 
$ The filter algorithm is only applied to LAC and HRPT data. 


the fourth reprocessing, another goal was to increase the 
number of retrievals without degrading the data quality by 
making the masking criteria less strict (known as relaxing 
the criteria). For flags that have a large role in reducing 
the number of retrievals, the algorithms could then be re- 
examined to see if the flagging could be relaxed to increase 
the number of retrievals. 

Early in the preparation for the fourth reprocessing, 
a test was performed to determine which flags, if relaxed, 
would result in a large increase in the number of retrievals. 
A level -3 8-day bin file was made for the period start- 


ing 2 June 2001 with minimal masking criteria (ATMFAIL, 
LAND, FILTER, CLDICE, CHLFAIL, NAVFAIL) enabled to form 
a control case for testing the effect of other flags in mask- 
ing data. The 11 flags usually used as masks in the level -3 
processing were individually turned on and the effect of 
each flag is shown in Table 7. The effect of the flags are 
listed in order from the largest number of bins (Earth area 
coverage) removed to the smallest. 

The stray light flag has the largest effect, decreasing 
coverage by 17% and the total number of retrievals by 
almost 33%. The geographic distribution of the bins lost to 
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Table 7. Effect of the level -3 binning flags in removing data. A baseline 8-day timebin level -3 file was produced 
using the third reprocessing and a minimum amount of level-3 masks (ATMFAIL, LAND, FILTER, CLDICE, CHLFAIL, 
NAVFAIL). The binning was repeated with one extra flag used as a mask to determine the effect of that flag in 
decreasing the bins (coverage) and samples relative to the baseline. The effects are listed in decreasing order of 
bins lost. 


Flag 

Name 

Decrease in 
Bins [%] 

Decrease in 
Samples [%] 

Description 

STRAYLIGHT 

17.26 

32.13 

Stray light contamination 

TURBIDW 

3.71 

7.93 

Turbid, Case-2 water 

L0WLW 

3.39 

8.32 

Low water-leaving radiance at 555 nm 

HIS0LZEN 

1.08 

0.48 

Solar zenith angle above the limit 

ATMWARN 

1.04 

1.43 

The e value unreasonable, Liv a t 490-555 nm < 0 

HISATZEN 

0.45 

1.64 

Satellite zenith angle above the limit 

HILT 

0.25 

0.35 

Total radiance above the knee 

ABSAER 

0.18 

0.51 

Absorbing aerosol index above threshold 

C0CC0LITH 

0.12 

0.39 

Coccolithophore bloom 

MAXAERITER 

0.08 

0.28 

Maximum iterations in NIR algorithm 

NAVWARN 

0.014 

0.028 

Questionable navigation (tilt change) 


the stray light flag is relatively uniform over the oceans and 
is caused mainly by cloud edges. Aggregation over longer 
time periods reduces the lost coverage because of changing 
cloud distributions. The degree to which coverage is lost 
because of the stray light flag prompted the project to re- 
examine the stray light algorithm and flag to see if the loss 
of coverage could be reduced. The results of this work are 
presented in Sect. 6.3.5. 

The analysis also revealed that the high satellite zenith 
angle flag was being set for some of the SeaWiFS GAC 
data. It was previously assumed that the 56° satellite 
zenith threshold would only flag LAC data, but analysis 
showed that it was happening at the very edges of the GAC 
scans. This prompted the examination of the high satel- 
lite zenith angle flag (Sect. 6.3.3). Other flags on the list 
were affected by the fourth reprocessing changes, but only 
the stray light and high satellite zenith angle flags were 
re-examined as a direct outcome of these tests. 

6.3 LEVEL-2 CHANGES 

The flag algorithms are all performed in the level -2 
processing, so a majority of the changes are discussed in 
this section. The effect of these changes on the data is also 
presented. Changes in the use of flags for accepting data 
for level-3 binning are discussed in Sect. 6.4. 

6.3.1 Glint 

In the third reprocessing, the sun glint flag was used 
as a mask in the level -2 processing. The flag was acti- 
vated when the estimated glint reflectance exceeded 0.005. 
Wang and Bailey (2001) showed that glint radiance could 
be effectively removed from the data by using the glint 
correction up to the point where the detectors saturate. 


Because other masks are active at this point, the glint flag 
no longer needs to be a mask. 

In the fourth reprocessing, the HI GLINT flag is changed 
from a mask to a flag in the level -2 processing to allow 
more good data to be retrieved. The change to the glint 
mask left the global parameter means unchanged. The 
number of retrievals increased about 1% while the coverage 
increased by 0.2% for the monthly binned data and 0.9% 
for the daily binned data. The recovered pixels are in areas 
where the sun glint is high, i.e., around the subsolar point. 

6.3.2 High Total Radiance Flags 

For the third reprocessing, the high total radiance 
(HILT) flag was set for a pixel if the digital count value 
was above the knee in any of the eight SeaWiFS bands. 
This was done because radiances taken above the knee have 
substantially lower radiometric precision than those taken 
below the knee (McClain et al. 1995) and were not be- 
lieved to have sufficient precision for accurate chlorophyll 
and water-leaving radiance retrieval. The flag has been 
used as a mask in the level -2 processing. 

It was determined that although the higher precision is 
required for the NIR bands 7 and 8 at 765 and 865 nm to 
accurately characterize the aerosols, it is not required for 
the other bands. For the fourth reprocessing, the HILT flag 
is defined so that it only applies to the radiances in bands 
7 and 8. 

This change yields substantially more good retrievals in 
some coastal areas. Over the globe, the estimated increase 
in number of retrievals is only 0.06%. When binned over 
a month, the coverage is increased by 0.02%. The new 
retrievals are in regions of high chlorophyll, which increases 
the global chlorophyll mean and Lwat(555), but reduces 
the Lwn{^ 12). 
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6.3.3 High Satellite Zenith Angle Flag 

The satellite zenith threshold was extended from 56° 
(used in the third reprocessing) to 60°. This allows pro- 
cessing to be performed for pixels on the edge of the GAC 
swath which were being excluded because of the high satel- 
lite zenith angle. The satellite zenith angle is affected by 
the spacecraft roll, which caused pixels on the edge of the 
GAC scene to have a satellite zenith angle slight larger 
than the 56° limit. The increase in number of retrievals is 
about 1.5% and is distributed evenly over the globe. 

6.3.4 Shallow Water Flag 

The shallow water mask used for the third reprocessing 
was based on the Earth Topography 5 min of arc grid — 
ET0P05 bathemetry data — from the National Geophysical 
Data Center (Sloss 1988). The mask for depths shallower 
than 30 m had a minimum resolution of about 9 km, which 
was noticable at the GAC resolution of 4 km and very not- 
icable at the LAC and HRPT resolutions of 1 km. 

The shallow water mask used for fourth reprocessing is 
based on the Earth Topography 2 min of arc grid — ET0P02 
bathemetry (Sloss 2001) — which has a resolution of ap- 
proximately 4 km. The improved resolution of the mask 
results in a smoother boundary between shallow and deep 
water. The 15 MB size of the new shallow water mask file 
is slightly greater than the previous file size of 13 MB and 
should not pose a storage space problem. 

6.3.5 Stray Light Flag 

Section 6.2 showed that the effect of the stray light flag 
in reducing the number of retrievals was the greatest of 
all the additional flags used in the level -3 processing. If 
some relaxation of the masking criteria was possible, the 
increase in retrievals would be the greatest. 

The stray light algorithm used through the third re- 
processing was developed using laboratory measurements 
of a bright target (Barnes et al. 1995). The correction 
factors (for pixels that could be corrected) and masking 
recommendations (for pixels that could not be corrected) 
were incorporated into the stray light algorithm (Yeh et al. 
1997). In the initial development of the corrections (Barnes 
et al. 1995), it was stated that it may be possible to relax 
the GAC masking so that two more pixels adjacent to a 
bright target could be corrected instead of flagged. Correc- 
tions for these pixels were also provided. The correcting of 
those extra pixels instead of masking results in a significant 
gain in retrievals. 

This change in the stray light algorithm and masking 
was tested on a level -3 data set and on several level -2 data 
sets in open and coastal waters. It was found that the 
change to the stray light algorithm increased the number 
of retrievals by about 17% and did not change the global 
means of retrieved parameters. It was also found that the 
new set of corrected pixels had chlorophyll and normalized 


water-leaving radiance values very similar to neighboring 
pixels in the scene, even though the uncorrected total radi- 
ances showed significant increases from neighboring pixels 
to the new pixels. 

The modified stray light algorithm is incorporated into 
the fourth reprocessing. In addition, the stray light flag is 
now used as a mask (instead of as a flag) for the level -2 
data sets. The exclusion of pixels that are highly contam- 
inated by, and uncorrected for, stray light in the level -2 
data results in more consistent data quality. The stray 
light algorithm is unchanged for both LAC and HRPT 
data. 

As a final note, the correction for stray light was de- 
rived only under laboratory conditions. These conditions, 
however, may not be representative of what is encountered 
by the satellite on orbit. It would be useful if an on-orbit 
analysis could be performed on the stray light correction 
factors to determine if they could be improved. The anal- 
ysis could be performed at known sharp radiance bound- 
aries, such as the moon’s limb in the lunar calibration data, 
or the boundary between the ocean and either land or ice 
features in HRPT data. 

6.3.6 Cloud and Ice Flag 

In September 2000, the SeaWiFS imaging duty cycle 
was increased to include additional data at higher lati- 
tudes. GAC data is now routinely collected at solar zenith 
angles as high as 83°. It was found that the standard 
SeaWiFS cloud detection algorithm was too restrictive at 
these extreme solar zenith angles. In fact, it was pointed 
out that the predicted Rayleigh path radiance alone is suf- 
ficient to trip the cloud flag when solar zenith angles exceed 
75° (H. Fukushima, pers. comm.). A proposed solution to 
this problem was to base the cloud test on a Rayleigh sub- 
tracted radiance value, such as the combined surface and 
aerosol reflectance at 865 nm. A threshold on this quasi- 
surface reflectance field, p s { 86 5), of 0.027 (2.7%) was found 
to yield very similar results to the third reprocessing algo- 
rithm at moderate solar and viewing geometries, while still 
allowing high solar- and viewing-angle observations of ap- 
parently clear ocean pixels to pass through unmasked. 

Figure 26 illustrates how the two quantities of albedo 
and quasi-surface reflectance, vary as a function of so- 
lar zenith angle. Each point in the plots is the average 
albedo or reflectance, computed from all the pixels in a 
typical GAC scan line. Only pixels for which there was a 
chlorophyll retrieval without stray light correction were in- 
cluded in the averages. The plots illustrate that the cloud 
albedo of relatively clear ocean observations (Fig. 26a) gen- 
erally increases with solar zenith angle, with the rate of in- 
crease rising rapidly above 65°. In contrast, the surface re- 
flectance distribution (Fig. 26b) remains relatively flat up 
to 75°, primarily because the change in the Rayleigh path 
radiance was accounted for. The high reflectance around 
40° is due to a dust plume. The solid horizontal lines show 
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Fig. 26. a) Albedo at 865 nm as a function of the solar zenith angle for a typical SeaWiFS GAC orbit. The 
horizontal line indicates the cloud mask threshold for the third reprocessing, b) Quasi-surface reflectance at 
865 nm for the same data as in panel a). The horizontal line indicates the cloud mask threshold for the fourth 
reprocessing. 


the threshold above which values are considered cloud con- 
taminated. The threshold on p 3 ( 865) was chosen to give 
similar results to the current albedo threshold at moderate 
solar and viewing angles, but the surface reflectance test 
will allow more observations at higher solar and viewing 
zenith angles to pass through unmasked. 

Application of the proposed cloud flag was found to 
increase overall coverage (in terms of number of filled 9 km 
bins) by 1-2%, with smaller gains seen in the monthly 
composite and larger gains seen in the daily composite. As 
expected, the increased coverage is seen primarily at high 
latitudes, where the solar zenith angle and the Rayleigh 
radiance are large. The average number of samples per 
bin also generally increased 1-2%. The most significant 
change to the derived products is a general increase in 
mean aerosol optical thickness, especially at higher solar 
angles. This is not suprising, as the cloud flag for the third 
reprocessing tends to mask even moderate aerosol loads 
when the solar and viewing angle is large, thus biasing the 
average toward lower aerosol optical thicknesses. 

6.3.7 Turbid, Case-2 Water Flag 

For the third reprocessing, the turbid water flag was 
set when the measured reflectance in band 5, at 555 nm, 
exceeded the estimated 555 nm reflectance derived from 
Morel’s clear water model (Morel 1988). Because the Morel 
model requires pigment as input, an invalid chlorophyll re- 
trieval would cause the turbid water flag to be erroneously 
set, or prevent the flag from being set under highly tur- 
bid conditions. This problem was remedied in the fourth 
reprocessing by basing the turbid water flag on a thresh- 
old of the remote sensing reflectance in band 6, at 670 nm 
[i? rs (670)]. The water is flagged as turbid if i? rs (670) is 


greater than 0.0012, which is a value 25% larger than the 
expected value for pure water of 0.000954. 

6.3.8 Chlorophyll Failure Flag 

The SeaWiFS chlorophyll fields, although generally 
smooth, have sporadic isolated high values, or speckles. 
In past reprocessings, the Project has implemented checks 
that would limit these speckles, but some still remain. 

For the fourth reprocessing, analysis of the data re- 
vealed some characteristic patterns in the water-leaving 
radiance spectra of pixels that have speckles. It was found 
that many of the high chlorophyll spikes were associated 
with spectra where the 510 or 490 nm bands went negative, 
while bracketing bands were positive. These conditions are 
considered nonphysical, and probably caused by proxim- 
ity to bright sources that cannot be seen in the subsam- 
pled GAC data and which may cause band-dependent stray 
light or ringing effects. To avoid these pixels, new failure 
conditions were implemented in the operational chlorophyll 
algorithm. The following tests must be satisfied or the 
pixel will be given the chlorophyll failure (CHLFAIL) flag: 

1) R rs (510) >0; 

2) R rs (490) > 0, or i? rs (443)R rs (490) > 0; and 

3) Both R rs (443) and R rs (490) > -0.001. 

The first two tests, taken together, enable the 443 nm and 
490 nm bands to be less than zero, but require that the 
shorter wavelengths become negative before the longer 
wavelengths. This typically happens in the blue part of 
the spectrum when a) there is aerosol absorption, or b) 
if an aerosol model is selected which removes too much 
radiance. The third test limits the amount by which the 
reflectances can go negative. These additional tests flag a 
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significant number of the pixels that would otherwise add 
speckles to the data. 

6.3.9 Absorbing Aerosol Flag 

The absorbing aerosol flag algorithm has remained un- 
changed. Some of the thresholds used in the algorithm 
required modification to account for the radiance changes 
at 412 and 443 nm which came as a result of the calibration 
changes described in Chapt. 3. 

6.3.10 Chlorophyll Warning Flag 

Previously, the chlorophyll was stored in the level -2 
data set as a scaled 2-byte integer which could represent 
any value from 0-64 mg m -3 with a precision of 0.001. For 
the fourth reprocessing, the chlorophyll parameter is stored 
as a 4-byte real value so that small values of chlorophyll 
can be stored with the same relative precision as larger 
values. This change also allows the storage of chlorophyll 
values greater than 64. The threshold for the chlorophyll 
warning flag, which was previously set at 64mgm" 3 , was 
increased to 100 mg m -3 to account for the increased range 
available. The limit of 100 was chosen so that high values 
of chlorophyll could be recorded but unreasonably high 
values would be flagged. 

6.3.11 Sea Ice Flag 

The sea ice flag, SEA ICE, is a new flag for the fourth 
reprocessing. It flags pixels that have a high probability 
of containing sea ice, based on a climatology of sea ice. 
The sea ice flag is used mainly in the creation of the PAR 
product. 

6.3.12 Navigation Failure Flag 

The navigation failure flag, NAVFAIL, is new for the 
fourth reprocessing. Previously, if a navigation failure con- 
dition was indicated in the level -la data, the atmospheric 
algorithm failure flag was set in the level -2 data set. There 
was only one such condition — a failure in the navigation al- 
gorithm. In the fourth reprocessing, it became necessary to 
automatically and manually mark time ranges in the data 
where the navigation was unacceptable because of timing 
errors, as described in Chapt. 7. It was decided that for 
this added condition, and to more clearly denote the navi- 
gation failure condition, a separate flag would be assigned 
to failed navigation conditions. The navigation failure flag 
is used as a mask in the level -2 processing. 

6.3.13 Aerosol Model Filter Flag 

The aerosol model filter flag, FILTER, is new for the 
fourth reprocessing. This flag supports a new algorithm 
which reduces the noise in the NIR band ratio and thereby 
reduces the small scale variability in aerosol model selec- 
tion. As a part of the filtering process, a minimum number 


of nonmasked neighboring pixels must be available for the 
filter to operate. If this minimum number is not present, 
the flag is set. The algorithm was not found to improve 
the GAC data, so it is not used in the fourth reprocess- 
ing for GAC data but it is suggested for LAC and HRPT 
processing. 

6.3.14 Ocean Flag 

The ocean flag, OCEAN, is also new for the fourth repro- 
cessing. It indicates that the given pixel is a clear ocean 
pixel, without cloud contamination or land. This flag is 
used to give an estimate of how much useful ocean data is 
in a level -2 file. 

6.4 LEVEL -3 CHANGES 

The following changes were made in the use of flags 
and masks for the binning process which makes level -3 
products. 

1. The glint flag is no longer a masking condition in 
the fourth reprocessing. Section 6.3.1 notes that the 
glint flag does not need to be a mask because the 
glint correction is valid up to where the high total 
radiance or cloud and ice flags mask the pixel. 

2. The chlorophyll warning flag is now a masking con- 
dition in the fourth reprocessing. This was done 
to keep erroneously high chlorophyll values out of 
the level -3 binned file. In the fourth reprocessing, 
chlorophyll values up to lOOmgm -3 are binned. 
More discussion about the chlorophyll flag can be 
found in Sect. 6.3.1. 

3. The navigation failure flag is new for the fourth re- 
processing. This flag is used to mask the level -2 
product (Sect. 6.3.12), and is carried on as a mask 
for the level -3 binning. 

4. The filter flag is set as a mask for the operational 
binning in the fourth reprocessing, but because the 
filtering is not enabled for GAC data, this mask has 
no effect on the operational level -3 product. 

6.5 SUMMARY 

A number of flag and mask changes were made as a part 
of the fourth SeaWiFS reprocessing. When combined with 
the calibration improvements and algorithm changes, this 
results in significantly better retrievals of water-leaving ra- 
diance, chlorophyll, and atmospheric products. A greater 
number of retrievals are also produced. 

The flagging changes resulted, in part, in an overall 
increase of retrievals by 24% and an increase in global cov- 
erage of 24, 14, and 6% for daily, 8-day, and monthly time 
binned level -3 products, respectively. This was made pos- 
sible by a relaxation in the stray light mask, the elimination 
of the glint flag as a mask, the improved treatment of the 
cloud and ice flag, a relaxation in the high satellite zenith 
flag, and a better treatment of the high total radiance flag. 
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Fig. 27. Global distribution of the gained bins and samples for a monthly level -3 binned product (May 
1999). The white indicates where no data exists, the grey indicates where the gain in samples in a bin was 
less than 30%, and the black indicates where new bin coverage occurred or where the number of samples in 
a bin increased by 30% or more. The continental outlines are drawn in black. 


Figure 27 is an image showing the locations where new 
bins were gained and where the number of samples in a 
bin increased by 30% or more. Throughout the ocean ar- 
eas, there have been relatively uniform increases in sam- 
ples and bins. This is mainly the result of the change in 
the stray light mask permitting more data around cloud 
edges. These increases can also be attributed to the gain 
of samples at the ends of the scan lines that were previously 
beyond the old satellite zenith angle cutoff. In addition to 
the uniform increases, several areas show higher amounts 
of new retrievals. The Southern Ocean has substantial in- 
creases in samples and bins, mainly resulting from the new 
cloud flag, which no longer masks out data having mod- 
erately high aerosol radiances. The coastal regions also 
show increases in samples and bins, resulting from the re- 
duction of the stray light mask near coastlines and the 
improved calibration. Other variations are a result of dif- 
fering amounts of clouds in parts of the world. 

Fields of normalized water-leaving radiance and chloro- 


phyll were also improved. The chlorophyll failure flag was 
improved to reject more speckles in the chlorophyll. 

Finally, the flags were improved and expanded to give 
better information about the conditions at any pixel. The 
shallow water flag is improved, using higher resolution 
data. The sea ice and ocean flags are new, and signal sea 
ice and unclouded ocean regions. The navigation failure 
flag separates the flagging for failed navigation conditions 
from different conditions. 

Possible future improvements in the flags include a more 
thorough examination of the stray light algorithm to verify 
and possibly improve the corrections and masking condi- 
tions. The flag for low water-leaving radiance at 555 nm, 
which is a flag for cloud shadows, may be re-examined 
to see if a better cloud shadow algorithm can be imple- 
mented. Users of the data from the fourth reprocess- 
ing should acquaint themselves with the new flagging and 
masking changes so they can make the best use of the Sea- 
WiFS data. 
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Level -la and Level -3 Processing Changes 

Frederick S. Patt, Bryan A. Franz, and Wayne D. Robinson 
Science Applications International Corporation 
Beltsville, Maryland 

Joel Gales 
Futuretech Corporation 
Greenbelt, Maryland 

Abstract 

Improvements were made in the level-la and level-3 data processing for the fourth reprocessing. Two of these — 
handling of the spacecraft time tag anomalies, and modification of the level -3 space binning program — address 
specific problems in the data and processing logic, which excluded otherwise valid data from processing. The 
third improvement, an update to the navigation algorithms, improves the overall data quality by reducing the 
maximum navigation errors. The net effect of these changes is an increase both in coverage and quality of the 
SeaWiFS data products. 


7.1 INTRODUCTION 

This chapter describes changes that were made in the 
level -la and level -3 data processing software as part of 
the fourth reprocessing. These changes were implemented 
independently of the level -2 algorithm updates which are 
described in the previous chapters. They represent either 
solutions to known problems in the previous versions of the 
processing software, or results of ongoing analysis of data 
quality. 

The first of these changes is a method for handling 
anomalies, or glitches, in the spacecraft time tag during 
level -la processing. These glitches have resulted in sig- 
nificant navigation errors, which previously could only be 
handled by excluding entire scenes from further processing. 
As a result, significant amounts of valid data were also ex- 
cluded to avoid including the misnavigated periods. The 
new scheme both corrects simple time glitches and allows 
specific, uncorrected periods to be flagged within scenes, 
minimizing the loss of data. 

As another part of the level-la processing, improve- 
ments were made to the navigation algorithms to reduce 
the seasonal variations in accuracy. This update was ac- 
tually made to the operational software in April 2001, but 
the fourth reprocessing was the first opportunity to apply 
it to the entire mission. In addition, a problem was iden- 
tified and corrected in the initial level -3 processing (space 
binning) logic, which excluded some level -2 files from bin- 
ning. Each of these changes is described in the following 
sections. 


7.2 TIME TAG GLITCH HANDLING 

The data stream produced by the OrbView-2 space- 
craft contains a time tag for each level -0 data record, or 
minor frame. (A level -0 minor frame contains either one 
LAC scan line or five GAC scan lines, along with asso- 
ciated spacecraft and instrument telemetry.) In addition, 
the telemetry data contain time tags for various space- 
craft subsystems. The spacecraft minor frame time tag is 
used during level -0 -to -1 A processing to compute the scan 
line time, which in turn is used in the determination of 
the spacecraft navigation fields (orbit position and atti- 
tude angles). 

Various types of time tag glitches have been observed 
during the mission. The most serious are those that involve 
the minor frame time tag. A 1 s error in this time tag re- 
sults in a pixel location error of 6.75 km, or 6 LAC pixels. 
Glitches that only affect other time tags (e.g., for teleme- 
try fields) have little or no effect, and are largely handled 
in the processing software by data quality checks. Note 
that the glitches discussed here originate on the spacecraft, 
and affect all data types — GAC, LAC, and HRPT — when 
they occur. Time tag glitches may also be caused by data 
transmission errors, but these glitches are source depen- 
dent (specific ground station, and recorded versus direct 
broadcast) and are handled by data quality checks. 

The types of minor frame time tag glitches are dis- 
cussed in more detail below, followed by a description of 
the procedure implemented for handling the time glitches. 
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7.2.1 Minor Frame Time Tag Glitches 

Two types of minor frame time code glitches have been 
observed in the OrbView-2 data. The first is associated 
with resets of the onboard Global Positioning System 
(GPS) receiver. The second is of unknown origin, and con- 
sists of a shift of exactly 1 s over a period of approximately 
30 s. Each of these is described below. 

The GPS receiver is the source of the spacecraft orbit 
data, and also provides a very accurate time signal. On 
OrbView-2, the GPS time signal is used as the reference 
for all other spacecraft time codes. At regular intervals, the 
time codes in the spacecraft computers are compared with 
the GPS time signal, and adjusted to maintain agreement. 

Resets of the GPS receiver have occurred throughout 
the mission. A reset occurs when the receiver temporarily 
stops tracking enough GPS satellites to get a valid orbit 
solution. These events were observed frequently at the 
start of the mission (September-December 1997). After 
the spacecraft engineers analyzed the problem, software 
patches were uplinked to the spacecraft computer and the 
GPS operating parameters were adjusted, which signifi- 
cantly reduced the number of events. Since January 1998, 
resets have been observed at very irregular intervals, av- 
eraging about one or two per month. Only resets which 
overlap with GAC data collection affect data quality. 

A GPS reset results in errors in the minor frame time 
tag, and other spacecraft time tags are also affected. The 
errors appear at the start of the reset event as a series of 
fractional-second steps, spread over about 1 min. In some 
cases, the steps cancel out each other. More frequently, 
they accumulate to a total error of 1 s; this remains for a 
period of from 15-20 min and is corrected in a single step 
at the end of the event. In these cases, because the event 
duration is a large fraction of a GAC scene length, it is 
common to observe only the start or end of an event. 

The second type of minor frame time tag glitch is called 
the 30s glitch, because of the characteristic duration of 
30-35 s. This glitch starts with a single shift of exactly 1 s, 
which is corrected from 30-35 s later. No other spacecraft 
time tags are affected. 

The occurrence of the 30 s glitches was not actually 
documented until early 2001, when they were observed by 
the quality control (QC) personnel on the CVT. That they 
went unobserved for over three years is most likely because 
of their brevity and (at first) infrequency, which limited the 
effect on data quality. Following their discovery, the mis- 
sion data set was analyzed to locate all occurrences. There 
were 17 events in 1998, 24 in 1999, 59 in 2000, 87 in 2001, 
and 48 as of the end of August 2002. The timing of the 
glitches was very irregular. In some cases, multiple glitches 
occurred in a single GAC scene, while at other times none 
were observed for weeks. Geographic and temporal analy- 
ses have revealed no discernable pattern. 

Prior to the fourth reprocessing, the only means for 
handling scenes affected by time tag glitches was for the 


QC team to manually fail entire scenes, to exclude them 
from both further processing and distribution. The de- 
cision of whether or not to fail a given scene was made 
according to the fraction of the data with navigation af- 
fected by the glitch, which was usually a joint decision by 
the QC and Mission Operations team. In the majority of 
cases (all GAC scenes and many HRPT scenes), each such 
decision caused some valid data to be lost as well. This 
approach was highly subjective and, as the frequency of 
the 30 s glitches increased, it also became a burden on the 
QC team. 

7.2.2 Time Tag Glitch Handling 

The approach to handling the minor frame time tag 
glitches during the level -0 -to -1 A conversion is based on 
two principles: 

1. Only those glitches which can be reliably detected 
and corrected, are corrected. 

2. The navigation flags in the level -la data products 
are used to flag scan lines affected by uncorrected 
glitches (as a navigation failure condition, which in- 
hibits downstream processing) and indicate those 
with corrected glitches (no effect on subsequent pro- 
cessing). 

The conditions in the first principle are met by the 30 s 
glitches in the GAC data. They can be detected reliably 
and automatically, because of their consistent character- 
istics (single step time tag shift of 1 s, corrected after 30- 
35 s), their limited effect (no other spacecraft time tags are 
affected), their brevity (for most glitches, both the start 
and end occur within a scene), and the magnitude of the 
shift compared to the GAC minor frame interval of 3.333 s. 
This last element is important because the frame interval is 
neither an integer multiple of, or smaller than, the time tag 
shifts, which can easily be distinguished from missing or 
repeated minor frames. Once the frames affected by a time 
glitch are identified, the error is easily corrected by adjust- 
ing all of the corresponding time tags by the same amount. 
(As indicated previously, a GAC minor frame contains five 
GAC scan lines, and the minor frame time tag is used as 
the first scan line time. The other four scan line times are 
extrapolated from the first using the scan line period of 
0.666 s. The minor frame time tag correction is, therefore, 
automatically applied to all five scan lines.) 

Unlike the GAC, the specific minor frames affected by 
a 30 s glitch in HRPT or LAC data cannot be reliably de- 
tected. The time interval for these frames is 0.166 s, so 
the time shift of 1 s is an integer number of minor frames. 
A forward time shift, therefore, is indistinguishable from 
a data gap, and a backward shift appears the same as 
repeated frames, which the level-0-to-lA software is de- 
signed to filter out. 

As stated previously, the time tag glitches occur si- 
multaneously in all data types (LAC, GAC, and HRPT); 
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therefore, it should be possible to determine the period 
affected by a 30 s glitch based on the G AC data. The lim- 
iting factor is the GAC frame interval, which is 3.333s; 
this means that the start and end times of the glitch pe- 
riod are uncertain by this amount, corresponding to 20 
LAC or HRPT frames. To address this, these frames at 
the start and end of each glitch period are flagged and not 
corrected. 

The following is the approach implemented for handling 
the 30 s glitches: 

1. A database table was developed to store time tag 
glitch information, consisting of the start and end 
times, the data type (GAC, or LAC and HRPT), 
the action to be taken (correct or flag minor frame 
times within the interval) and the time shift (if cor- 
rection is indicated). 

2. During GAC data processing, any 30 s glitches are 
detected, and the affected minor frame time tags 
are corrected prior to the calculation of the scan line 
time tags. (Note that the original, uncorrected mi- 
nor frame time tags are also stored with the data.) 
The glitch time period is added to the database ta- 
ble. 

3. Database table entries are also made for LAC and 
HRPT data as follows: one entry to correct the 
minor frames within the time range of the corrected 
GAC frames, one to flag the 20 minor frames prior 
to this time range, and one to flag 20 minor frames 
after. 

4. During subsequent LAC or HRPT data processing, 
the table entries are used to correct or flag minor 
frame time tags. 

The GPS reset time tag glitches do not meet the reliable 
detection and correction criteria, for the following reasons: 

• Their start is reflected in multiple time shifts, so 
there is no single start time for a glitch; 

• Their duration causes a significant number of these 
glitches to span the start or end of a GAC scene, 
making unambiguous identification difficult; 

• The shifts affect multiple spacecraft time tags, par- 
ticularly in the attitude control system (ACS) tele- 
metry, and; 

• The shifts do not occur simultaneously, making re- 
liable correction virtually impossible. Detection of 
GPS reset glitches is still essentially a manual pro- 
cess. 

The approach described above can also be used to flag 
specified periods with no associated correction, and this 
approach was chosen for the GPS reset glitches. As GPS 
resets are detected by the QC or Mission Operations teams, 
entries are manually added to the database table to flag the 
affected time periods. Afterward, any affected data sets 
previously generated are renavigated to apply the flags, 


and all subsequently processed data are flagged automati- 
cally. While this process still requires manual intervention, 
it has the advantage of allowing only the specific data af- 
fected by a glitch to be excluded from further processing, 
instead of an entire scene being failed. 

7.3 NAVIGATION UPDATE 

The SeaWiFS navigation processing is performed as 
part of the level -0-to-lA conversion. The results of nav- 
igation processing have been under continued evaluation 
since launch, and numerous improvements were made prior 
to the third reprocessing. Following that reprocessing (May 
2000), the navigation results for the entire mission up to 
that time were analyzed, and the results of that analysis 
were used to develop the modification discussed here. This 
change was actually incorporated into the operational pro- 
cessing in April 2001, but was not applied to the mission 
data prior to that time until the fourth reprocessing. 

The SeaWiFS navigation accuracy requirement is 1 pix- 
el (1.1km at nadir) at twice the standard deviation. The 
accuracy of the navigation processing was evaluated using 
the method of island targets, which was developed prior 
to launch (Patt et al. 1997). During the mission, the re- 
sults of this method have been used to develop modifica- 
tions which have steadily improved the navigation accu- 
racy (Patt 1999). The current set of navigation algorithms 
are described in another volume of this series (Patt 2002). 

The results from the third reprocessing showed good 
overall accuracy, but significant seasonal and geographic 
variations were observed. Specifically, errors in the South- 
ern Hemisphere during February and March frequently 
exceeded 2 pixels, because of excessive yaw angle errors. 
These errors were observed, for example, in the HRPT 
scenes from the Pretoria ground station. With the algo- 
rithms in use at that time, however, attempts to reduce the 
yaw errors for this period (e.g., by adjusting attitude sen- 
sor alignments) would have degraded the accuracy during 
other seasons. 

The results also showed that, while the roll angle er- 
rors were within the requirement, there were consistent, 
seasonal variations. Specifically, the errors versus latitude 
during the months of November-February had a small but 
distinct curvature, with positive errors near the equator 
decreasing to negative toward the poles. Figure 28 shows 
the average roll and yaw errors versus latitude for HRPT 
data from the last 10 days of February 2001. This illus- 
trates the typical behavior for this period: both the roll 
error peak near the equator, and the large yaw errors in 
the Southern Hemisphere. 

Because the roll angle is essentially determined from 
the spacecraft horizon scanner data, the curvature observed 
during these months suggested a seasonal variation in the 
shape of the CO 2 layer of the atmosphere as seen by the 
scanner (Patt and Bilanow 2001). The CO 2 height model 
used for the third reprocessing already had a seasonal cor- 
rection for the height and the north-south shift, but the 
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Fig. 28. HRPT data average for the last 10 days of February 2001: a) roll errors, and b) yaw errors. 


shape (as represented by the ellipsoidal flattening factor) 
was a constant which had been determined for all seasons. 
An additional correction was added to the model to vary 
the flattening factor seasonally, with the effect of mak- 
ing the atmosphere slightly more ellipsoidal in November- 
February than in May- August. 

The revised model had the desired effect of minimizing 
the seasonal roll error variations. It was also found that by 
incorporating small changes to the attitude sensor align- 
ments, it was possible with the revised model to reduce 
the maximum yaw errors in the February-March period, 
without degrading the accuracy at other times. This com- 
bination of changes was implemented in April 2001 and 
used for the fourth reprocessing. The accuracy was eval- 
uated using the island targets method for the entire GAC 
mission data, and shows the yaw angle errors to be within 
acceptable limits at all seasons. The HRPT mission data 
have not yet been reprocessed, but the results for February 
and March 2002 showed the desired improvement in the 
Southern Hemisphere. The average HRPT roll and yaw 
errors for the last 10 days of February 2002 (Fig. 29) show 
that the roll errors have been flattened out and, more im- 
portantly, the Southern Hemisphere yaw errors have been 
significantly reduced. 

7.4 SPACEBIN MODIFICATIONS 

This section describes a problem in the initial level -3 
processing that was first identified and resolved during the 


preparations for the OCTS GAC data processing, per- 
formed in November 2001. The problem and its solu- 
tion were also determined to be applicable to SeaWiFS 
GAC processing, and this modification was adopted for 
the fourth reprocessing. 

The level -3 ocean products generated by SeaWiFS and 
SIMBIOS follow a convention known as the data day. The 
data day is defined geographically, according to the time 
when the spacecraft orbit crosses the 180° meridian (here- 
after called the Date Line ) closest to the equator (Podesta 
1995). 

To apply the data day convention, the level-2 ~to _ 3 
binning program, spacebin, splits any level -2 scenes that 
cross the Date Line into two level -3 files. The pixels that 
are west of the date line are assigned to a data day one 
later than those east of the line. The determination of the 
appropriate data day involves two separate steps. A pro- 
gram is run to determine the start and end times of the 
current data day, using the Coordinated Universal Time 
(UTC) day of the scene as a starting point. Then the bin- 
ning program logic determines whether each pixel should 
be associated with the current, previous, or following data 
day, based on the observation time and location. 

During the testing of this processing for the OCTS 
GAC data, a problem was found in which scenes near the 
middle of the data day had alternate data days assigned 
incorrectly. At the time of spacebin execution, each scene 
has been assigned a primary data day, as described in the 
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Fig. 29. HRPT data average for the last 10 days of February 2002: a) roll errors, and b) yaw errors. 


previous paragraph. The program then assigns either the 
previous or following data day as an alternate, to be used 
if the scene crosses the Date Line. This assignment is per- 
formed using the scene time and the data day start and 
end time. If the scene time falls within the first half of 
the data day, the previous day is selected as the alternate; 
otherwise the following day is selected. 

During the test runs, it was found that in December 
and January, scenes near the data day midpoint were some- 
times not being binned. Two problems were identified: 

1. The selection of the alternate data day was being 
performed using the scene start time, which biased 
the selection toward the previous day. 

2. The selection was being performed by first checking 
if the scene time was within 12 h of the data day 
start time. The data day start and end times are not 
exactly 24 h apart, because there is not an integral 
number of orbits per day. For data days less than 
24 h long, this logic also biased the selection toward 
the previous day. 

This caused problems in December and January. The 
low solar declination caused the data collection to extend 
well past the South Pole, which significantly increased the 
longitudinal coverage. Each scene spanned nearly the en- 
tire daylight part of the orbit. As a result, many scenes col- 
lected near 1200 UTC extended from the Greenwich merid- 


ian past the Date Line. The incorrect selection of the pre- 
vious data day as the alternate for these scenes resulted in 
no data being binned. The solution is twofold. First, the 
scene center time is used for this selection, rather than the 
start time. Second, the selection is performed using the 
actual midpoint of the data day, rather than 12 h from the 
data day start time. 

An example of a scene that is affected by this change is 
S20020 10 114002 . L2..GAC. The start and end times for data 
day 2002010 are 00:17:59 and 23:38:16 (all times are UTC). 
Using the previous logic of spacebin, the scene start time 
(11:40:02) was within 12 h of the data day start time, and 
so the alternate data day selected was the previous day. 
This scene crossed the date line at its southern end, re- 
sulting in observations that belonged in the following data 
day. The net result is that no data from this scene were 
binned for data day 2002010. Using this modified logic, 
the scene center time (12:01:52) is compared with the data 
day midtime (11:58:07). Because the scene center time is 
later, the following data day is selected as the alternate, 
and the data are binned correctly. 

These modifications were implemented in the opera- 
tional version of spacebin and were used for the fourth re- 
processing. The result is that no scenes have failed space 
binning because of incorrect selection of the alternate data 
day. 
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Abstract 

A new SeaWiFS derived product was developed which provides an estimate for the amount of PAR reaching 
the ocean surface over a 24 hr period. A description of the algorithm is provided in this chapter, followed by 
comparisons with in situ observations. The in situ observations include several years of data covering a wide 
range of solar illumination conditions. The results indicate good algorithm performance, with RMS differences 
between satellite-retrieved and observed daily average PAR within a few einsteins per square meters per day. 


8.1 INTRODUCTION 

The SeaWiFS PAR product is an estimate of daily (i.e., 
24 hr averaged) PAR reaching the ocean surface. PAR is 
defined as the quantum energy flux from the sun in the 
spectral range of 400-700 nm. It is expressed in einsteins 
per square meter per day. In the sections that follow, the 
algorithm is described in detail and comparisons are pro- 
vided with in situ observations from two moored buoys, 
one at the relatively high latitude of British Columbia, 
and the other in the equatorial Pacific. 


trial solar irradiance at the time of the observation, t d is 
the clear sky total (direct plus diffuse) transmittance, t g is 
the gaseous transmittance, S a is the spherical albedo, and 
A is the cloud-surface system albedo. As the irradiance, 
E 0 cos 0 o t d t g / (l — S a A), passes through the cloud-surface 
system, it is further reduced by (1 — A)(l — A s )~ x . The 
solar irradiance reaching the ocean surface is then given by 


E c (l-A) 

(1 — A s )(l — S a A) J 


(15) 


8.2 ALGORITHM DESCRIPTION 

The PAR algorithm uses plane-parallel theory and as- 
sumes that the effects of clouds and clear atmosphere can 
be decoupled. The planetary atmosphere is, therefore, 
modeled as a clear sky atmosphere positioned above a 
cloud layer. This approach was shown to be valid by 
Dedieu et al. (1987) and Frouin and Chertock (1992). The 
great strength of such a decoupled model resides in its sim- 
plicity. It is unnecessary to distinguish between clear and 
cloudy regions within a pixel, and this dismisses the need 
for assumptions about cloud coverage distribution. 

For a solar zenith angle 0 O > the incoming solar irradi- 
ance at the top of the atmosphere, E 0 cos6q is diminished 
by a factor t d t g /(l — S a A) by the time it enters the cloud- 
surface system. In these expressions, E 0 is the extraterres- 


where A 3 is the albedo of the ocean surface and E c = 
E 0 cos 6 0 tdt g is the solar irradiance that would reach the 
surface if the cloud-surface system were nonreflecting and 
nonabsorbing. In clear sky conditions, A reduces to A s . 

To compute E 3i A is expressed as a function of the 
radiance measured by SeaWiFS in the PAR spectral range 
(i.e., in bands 1-6, nominal center wavelengths from 412— 
670 nm). The algorithm works pixel by pixel and proceeds 
as follows. 

First, for each pixel not contaminated by sun glint, the 
SeaWiFS observed radiance in band i at the top of the 
atmosphere, L t (\i) y is transformed into reflectance, />t(Ai), 
as 


PtPd 


nLtjXi) 

£ o (Ai)cos0 o ’ 


(16) 
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where Eo(Ai) is the extraterrestrial solar irradiance in band where A* and As are nominal wavelengths for SeaWiFS 
z, and 9q is the solar zenith angle at the SeaWiFS obser- bands i and 8, respectively. A monthly climatology based 
vation time. The glint areas are not selected because they on three years of SeaWiFS data (1997-2000) is used for 
would be interpreted as cloudy in the PAR algorithm. r a (A 8 ) and a, because aerosol properties cannot be deter- 
Second, p t ( A*) is corrected for gaseous absorption due mined when the pixel is cloudy. This procedure is also 

justified because, in general, aerosol effects on E s are sec- 

( 17 ) ondaxy compared to cloud or 8o effects. 

To estimate u a and P a , the two closest SeaWiFS aerosol 

models, k and /, that verify a(l) < a < a(k) are selected, 

( 18 ) and a distance d a = [ a(l)—a ] [a(l) -a(fc)] is computed. 
Using this distance, u a and P a are obtained as follows: 


to ozone 


with 


= 

toz(\) = expj — 


Pt( Ai) 


^oz(Ai) 

— Q oz (Xj)C 0 

COS 9 0 


where a oz (A*) is the ozone absorption coefficient in band i 
and C oz is the ozone concentration. The ozone concentra- 
tion is obtained from time and space coincident ancillary 
data sets distributed by the SeaWiFS Project (Ainsworth 
and Patt 2000). 

Third, the reflectance of the cloud-surface layer, p s (Xi), 
is obtained from p[(Xi) following Tanre et al. (1979) as- 
suming isotropy of the cloud-surface layer system. The 
reflectance is calculated using: 


Ps{Xi) — [p*(Ai) Patm(Ai)] ^d(Ai, 0o)td(\iy 6 V ) 


. (19) 


where 6 V is the viewing zenith angle and p a tm(Ai) is the 
intrinsic atmospheric reflectance in band i (photons that 
have not interacted with the cloud-surface layer). The 
assumption of isotropy is made because no information on 
pixel composition is available. 

In (19), Patm(Ai) is modeled using the quasi-single scat- 
tering approximation: 


7m P m -h (V a T a P a 
4 COS #o COS 9 V ’ 


(20) 


u a = d a u> a (k) + [1 — d a \cj a (l), (24) 


and 

Pa = d a P a {k) + [\-d a \P a {l), (25) 

where u> a {l) and uj a (k) are the single scattering albedos of 
aerosol models l and &, and P a (Z) and P 0 (fc) their respec- 
tive phase functions. 

Next, an estimate of daily PAR, E s , is obtained by 
integrating (15) over the length of the day: 


with 



cosQQtgt<i[l - A] A 
[1-4] [l — S a A\ 


(26) 


to — toz tv 


(27) 


_ Si [UMEoM] 

d Ei^o(Ai) 

JJ [gq(AQ^(A < )] 

° " EiEo(Xi) 


(28) 

(29) 


where r m and r 0 are the optical thicknesses of molecules 
and aerosols, respectively, P m and P a are their respective 
phase functions, and u a is the single scattering albedo of 
aerosols. Wavelength dependence has been dropped for 
clarity. The total transmittance t^, and spherical albedo 
S a , are computed using analytical formulas developed by 
Tanre et al. (1979): 


td{8) = exp^ 


-{r m + r a ) | f 0.52r m + 0.83r a 


cos 6 


exp 


cos 0 


, (2: 


and 

S a = (0.92r m + 0.33r a )e- (Tm+T “ ) , (22) 

where 6 is either 9 0 or 9 V . 

The optical thickness of aerosols in band z, t 0 (A*), is 
obtained from the optical thickness in band 8, r a (A 8 ), and 
the Angstrom coefficient, a, between bands 4 and 8: 


A 3 - 


0.05 

1.1 cos 1 *^ $o H - 0.15 

+ 0.0 8 ( 1 -g), 


fd 





fd 

Td( Ai) = exp 


LMXQEojXi) 
ZiEo(Xi) ’ 


T m ( Ai ) + T 0 (Ai) 
COS 


A = Fp s , 


EiPsi Ai) 

E i [p i (A i )^ 0 (A,)]’ 


(30) 

(31) 

(32) 

(33) 

(34) 


T a (Xi) 


r«( A 8 ) 


As 

Ai 


(23) w h ere horizontal line above the variable symbolizes the 
average value over the PAR range of 400-700 nm, and Eo is 
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Fig. 30. Scatterplot of daily averaged SeaWiFS PAR against in situ PAR from moored buoys near a) British 
Columbia (Halibut Bank) and b) the equatorial Pacific (epl). 
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Fig. 31. Scatterplot of 8-day averaged SeaWiFS PAR against in situ PAR from moored buoys near a) British 
Columbia (Halibut Bank) and b) the equatorial Pacific (epl). 
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Fig. 32. Scatterplot of monthly averaged SeaWiFS PAR against in situ PAR from moored buoys near a) 
British Columbia (Halibut Bank) and b) the equatorial Pacific (epl). 
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Table 8. Performance of the SeaWiFS PAR algorithm at the Halibut Bank and epl mooring sites. 


Averaging Period 

Daily 

8-Day 

Monthly 


Halibut Bank 

Correlation coefficient, r 2 
Biasf 

RMS differencef 
Meanf 

Number of points 

0.904 

0.932 (3.3%) 
6.2 (21.7%) 
28.4 
505 

0.984 

0.863 (3.1%) 
2.3 (8.2%) 
28.2 
54 

0.994 

1.10 (4.1%) 
1.8 (6.5%) 
27.2 
24 


epl 

Correlation coefficient, r 2 
Biasf 

RMS differencef 
Meanf 

Number of points 

0.613 
2.9 (6.0%) 
6.2 (12.8%) 
48.7 
882 

0.68 

2.8 (5.8%) 
4.3 (8.9%) 
48.3 
103 

0.673 

2.8 (5.8%) 

3.9 (8.0%) 
49 

38 


Halibut Bank and epl Combined 

Correlation coefficient, r 2 
Biasf 

RMS differencef 
Meanf 

Number of points 

0.883 

2.2 (5.3%) 

6.2 (15.0%) 
41.3 
1387 

0.957 
2.1 (5.2%) 
3.7 (9.1%) 
41.4 
157 

0.978 

2.2 (5.4%) 

3.3 (8.0%) 
40.6 
62 


f In units of einstein m 2 day 1 . 


the mean extraterrestrial solar irradiance within that spec- 
tral range. The direct transmittance through the atmos- 
phere is denoted as T^. The summation over wavelength 
indicated in (28), (29), (31), and (34) is performed over 
SeaWiFS bands 1-6. 

The integration over time expressed in (26) is performed 
from sunrise to sunset of the day and location of the Sea- 
WiFS observation. It should be noted that, with the ex- 
ception of Eq and 5 a , all terms within the integral will 
vary with time of day in accordance with the change in 
solar zenith angle. 

In the expression of t g (27), the effect of both ozone 
and water (£ oz and i y/V) respectively) is modeled according 
to Frouin et al. (1989). 

Following Briegleb and Ramanathan (1982), surface 
albedo (30) is parameterized as a function of solar zenith 
angle and fractions of direct and diffuse incoming sunlight, 
with the diffuse component fixed at 0.08. The direct and 
diffuse components are weighted by the respective fractions 
of direct and diffuse transmittance. This parameterization, 
which takes into account Fresnel reflection and diffuse up- 
welling radiation, is sufficient because the influence of A s 
on surface PAR is small. In some cases, however, the re- 
trieved A might be less than A s . When this happens, A is 
fixed to A s . 

The dependence of A on solar zenith angle (33) is taken 
into account via the angular factor, F . Rather than using 
angular models determined statistically (e.g., Young et al. 
1998), F is computed from analytical formulas proposed by 
Zege et al. (1991) for nonabsorbing, optically thick scat- 


tering layers. The available angular models are fairly sim- 
ilar for partly cloudy, mostly cloudy, and overcast condi- 
tions, and they compare reasonably well with the Zege et 
al. (1991) formulas. 

The cloud-surface system, however, is assumed to be 
stable during the day and to correspond to the SeaWiFS 
observation. This assumption is crude, and PAR accuracy 
will be degraded in regions where clouds exhibit strong di- 
urnal variability. Still, useful daily PAR estimates would 
be obtained by averaging in space and time. Note that the 
algorithm yields a daily PAR estimate for each instanta- 
neous SeaWiFS pixel. 

Finally, the individual daily PAR estimates, obtained 
in units of mW cm“ 2 pm -1 , are converted into units of ein- 
steins per square meter per day (einstein m" 2 day -1 ) and 
averaged into 9 km resolution, daily, weekly, and monthly 
products. The factor required to convert units of mW cm" 2 
pm -1 to units of einstein m" 2 day -1 is equal to 1.193 to 
an inaccuracy of a few percent regardless of meteorological 
conditions (Kirk 1994). In middle and high latitudes, sev- 
eral daily estimates may be obtained over the same target 
during the same day, increasing product accuracy. 

8.3 IN SITU MATCH-UP COMPARISON 

An evaluation of the SeaWiFS PAR estimates was per- 
formed using several years of in situ PAR measurements 
from moored buoys off the west coast of Canada (Hal- 
ibut Bank data set, 49.34°N-123.73° W) and in the central 
equatorial Pacific (epl data set, 0.00°N-155.00°W). The 
total number of days used in the evaluation is 1387 (882 for 
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epl and 505 for Halibut Bank). Scatterplots of SeaWiFS 
versus in situ values are displayed in Figs. 30, 31, and 32 for 
daily, weekly, and monthly averages, respectively, and com- 
parison statistics are summarized in Table 8. Agreement 
with in situ measurements is good, with RMS differences 
of 6.2(15.0%), 3. 7(9.1%), and 3.3(8. 1%) einsteinm -2 day -1 
on daily, weekly, and monthly time scales, respectively, 
when the epl and Halibut Bank data sets are combined. 

The SeaWiFS estimates are higher by approximately 
1 einsteinm -2 day” 1 at Halibut Bank and by approximate- 
ly 3 einsteinm” 2 day -1 at the epl location. Overestima- 
tion at the epl location is due to less cloudiness at local 
noon (about the time of satellite overpass) than during 
the afternoon. A further verification was made using 16 
days of data collected at the Bermuda Bio-optical Program 
(BBOP) site. Similar statistics were obtained for daily val- 


ues, i.e., an RMS difference of 5.6(16%) einsteinm” 2 day -1 
and a negligible bias. 

The results presented above indicate good algorithm 
performance. One should be aware of the limitations of the 
algorithm, which ignores the diurnal variability of clouds. 
This variability will be introduced statistically, as a func- 
tion of geographic location and month of year, in a future, 
improved version of the algorithm. 
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Abstract 

The global atmospheric correction algorithm for SeaWiFS tends to over correct for the atmosphere in coastal 
waters because of water-leaving radiance ( Lw ) in the NIR part of the spectrum, A r . This Lw(K) phenomenon 
occurs particularly in water with high inorganic particulate levels. An iterative solution is used to solve this 
problem. A bio-optical model is used to determine the NIR backscatter from the backscatter at 670 nm, and 
specifically addresses inorganic particulates. This solution requires compensation for absorption by chlorophyll, 
detrital pigments, and gelbstoff (colored dissolved organic matter). The Lw(K) is found and removed from 
the total radiance so that the standard atmospheric models can be applied. Chlorophyll a concentrations, C a , 
in coastal and Case-2 waters are reduced to appropriate levels. The algorithm cannot yet correct areas where 
negative Lw occurs at 670 nm. 


9.1 INTRODUCTION 

The use of satellites to monitor the color of the ocean 
requires effective removal of the atmospheric signal. The 
methods for treating the atmosphere have depended on the 
high absorption of red and NIR light by water. In open 
ocean water where only C a and related pigments deter- 
mine the optical properties, water can be considered to 
absorb all light, so that the signal observed by the satellite 
should result entirely from the atmospheric path radiance. 
While the scattering due to Rayleigh, and absorption due 
to ozone and other gases can be treated through compu- 
tation with appropriate lookup tables to address seasonal 
and latitudinal effects (Gordon et al. 1983 and Gordon and 
Wang 1994), the aerosol optical depth must be computed 
for each pixel. 

The aerosol correction has required determination of 
two major factors: 

1) The amount of aerosol, characterized by the optical 
depth; and 


2) The type of aerosol, which determines the size dis- 
tribution and apparent color, and is characterized 
by either the Angstrom exponent, or e. 

The atmospheric correction for the Coastal Zone Color 
Scanner (CZCS) used a band at 665 nm (CZCS band 4) 
to provide a correction for the aerosol optical depth (Gor- 
don et al. 1983), with additional assumptions about the 
spectra for water with negligible C a (Gordon and Clark 
1981), and calculates the e value for each scene (but not 
each pixel) . The 665 nm band is positioned where the ab- 
sorption of water becomes significant (Table 9), thus, the 
entire radiance in the CZCS correction is presumed to orig- 
inate in the atmosphere. 

It was recognized immediately that water with a de- 
tectable scattering component, in particular inorganic sed- 
iment, had detectable radiance at this wavelength. Smith 
and Wilson (1981) proposed an iterative solution based on 
assumptions of the spectral relationship between 443, 550, 
and 670 nm. Later, Mueller (1984) proposed another so- 


51 





Algorithm Updates for the Fourth SeaWiFS Data Reprocessing 


lution using C a to estimate the radiance at 670 nm. In 
both solutions, the estimated water-leaving radiance (Lw) 
at 670 nm is removed from the signal, then the normal- 
ized water-leaving radiances, Lwn , are recalculated and 
the Lw values are re-estimated. When the change in Lwn 
values becomes negligible, the iteration ends. Gould and 
Arnone (1994) altered the approach by using estimates 
of the diffuse attenuation to obtain the estimated Lw , 
and significantly improved the usefulness of CZCS data 
in clearer Case-2 waters. While these are the best ap- 
proaches for CZCS, these iterations are based on in- water 
optical properties which are limited by their dependence on 
consistent spectral relationships between the bands. These 
spectral relationships will vary with changes in the optical 
constituents. In addition, use of the 670 nm band for at- 
mospheric correction can produce unpredictable results as 
Lw(6 70) can also become greater than the aerosol radiance 
[L a (670)], making the result unreliable. The iteration fur- 
ther depends on consistent spectral relationships between 
the bands. These relationships can vary owing to changes 
in the optical constituents. 


Table 9. Water absorption at the SeaWiFS bands 
(Curcio and Perry 1951, Palmer and Williams 1974, 
Smith and Baker 1981, and Pope and Fry 1997). 


Band 

No. 

Wavelength 

[nm] 

Bandwidth 

[nm] 

Absorption 

[m->] 

1 

412 

20 

0.00450 

2 

443 

20 

0.00700 

3 

490 

20 

0.01500 

4 

510 

20 

0.03250 

5 

555 

20 

0.06000 

6 

670 

20 

0.43000 

7 

765 

40 

2.50000 

8 

865 

40 

4.30000 


With SeaWiFS, bands for atmospheric correction were 
introduced in the NIR at 765 and 865 nm (Table 9). The 
absorption of water at these bands is sixfold and tenfold, 
respectively, of the absorption at 665 nm, producing neg- 
ligible water-leaving radiance in most Case-1 waters. In 
addition, the two bands offer a means of deriving an atmos- 
pheric model that adjusts for aerosol type by determining 
e(Ai,Aj) at each pixel, where e relates the aerosol radi- 
ance at band A i to a reference band A j (Gordon and Wang 
1994). In coastal waters having high concentrations of 
scattering material, water-leaving radiance can still occur 
in the NIR, an effect noted previously by researchers using 
other sensors that had NIR bands [Moore (1980), Stumpf 
and Tyler (1988), and Stumpf and Pennock (1989)]. 

The presence of water-leaving radiance in the NIR in- 
troduces two sources of error into the removal of the aerosol. 
First, the total aerosol is overestimated as some of the total 
radiance ( L t ) at 865 nm derives from the water. Second, 
as the absorption of water changes from 765-865 nm, the 


selection of the appropriate atmospheric model is affected, 
causing an error in the extrapolation of the aerosol radi- 
ance to the shorter wavelengths. As a result, the atmos- 
pheric radiances will be overestimated at all bands with 
increasing severity for shorter wavelengths, even leading 
to negative radiances in the blue bands in coastal water. 
This results in severe errors, if not complete failure, of var- 
ious algorithms for C a and optical properties. 

To solve this problem, iterative solutions were proposed. 
Land and Haigh (1996) developed a solution for SeaWiFS 
that involved modeling the water reflectance and atmos- 
pheric aerosol at all wavelengths to convergence. They 
attempted to solve simultaneously for both the ocean and 
atmosphere at all wavelengths. While a promising solu- 
tion, this attempt puts severe demands on the accuracy of 
the bio-optical model. 

For the Medium Resolution Imaging Spectrometer 
(MERIS), Moore et al. (1999) developed an iteration based 
on three NIR bands to estimate the NIR water radiance 
based on the Gordon and Wang (1994) solution. This so- 
lution is promising, but does not apply to SeaWiFS or 
MODIS, both of which have fewer NIR bands. 

Gould et al. (1998) proposed an iteration for SeaWiFS 
that determined NIR scattering from 670 nm scattering, 
which is the basis for the scattering component presented 
here. The solution addressed only the absorption by water. 

Ruddick et al. (2000) developed a method that solves 
for the aerosol radiance and Lw simultaneously in the NIR 
to good results in highly scattering systems in the North 
Sea. This method uses a single aerosol type (constant e 
value) determined manually, which poses problems for au- 
tomated processing. 

Hu et al. (2000b) developed a technique to transfer the 
e value from the nearest clear water. While potentially 
effective, it depends on the aerosol type remaining spatially 
constant in the coastal zone, potentially over hundreds of 
kilometers. 

Siegel et al. (2000) developed an iterative technique 
that presumed the backscatter to covary with C a and used 
the Gordon and Wang model to solve for the NIR backscat- 
ter. This model has been implemented by NASA, but is 
most effective in Case-1 water or in water where sediment 
covaries with C a . 

The information presented in this chapter is a partially 
coupled solution for SeaWiFS (and ultimately MODIS), 
where the scattering problems of atmosphere and water 
scattering in the NIR are coupled. A variation on this 
model was implemented into the atmospheric correction 
program within the processing software for the fourth re- 
processing of the SeaWiFS data set (Chapt. 4). 

9.2 METHODS 

The development of an algorithm that couples to the 
Gordon and Wang (1994) atmospheric correction involves 
several components. The first involves identifying the nec- 
essary change in the atmospheric correction. The second 
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involves the theory and bio-optical models used to deter- 
mine the water-leaving radiance Lw( V), and third is the 
iterative process. 


to the linear form 


^rs(A) — y 


T w 65 (A) 
0(A) a** (A)’ 


(38) 


9.2.1 NIR Concept 

Gordon and Wang (1994) (henceforth referred to as 
GW94) describe the solution for L t {\): 

Lt( A) - L r ( A) 4- [La( A) + Lra(A)] + t(A)Lwr(A) # (35) 

where L t ( A) is the radiance at the top of the atmosphere, 
L r ( A) is the Rayleigh scattering radiance, L a ( A) is the 
aerosol scattering radiance, L ra ( A) is the interaction be- 
tween molecular and aerosol scattering radiance, t{ A) is the 
diffuse transmittance through the atmosphere, and Lw(X) 
is the water-leaving radiance. A basic assumption of the 
GW94 atmospheric correction approach is that Lw(A r ) is 
negligible. This assumption allows for the selection of an 
aerosol model using 

La( At) + L ra (A r ) = L t ( A r ) - L r (A r ). (36) 

If Lw(T6 5) and Lvv(86 5) are not negligible, the right 
side of (36) is increased, introducing two different errors 
in the determined aerosols. First, if Lw(865) is not neg- 
ligible, then L a (865) + L ra (865) is overestimated. This 
means that more aerosol is determined than is present, 
resulting in overcorrection for aerosols at all bands, pro- 
ducing Lw(X r ) to be lower than is the case. Second, with 
Lvv(76 5) not negligible, the aerosol type will be in error. 
The absorption by water at the SeaWiFS bands varies (Ta- 
ble 9), and at 765 nm it is 57% of that at 865 nm. As a 
result, Lvv(765) > Lw(86 5), and water-leaving radiance 
that is interpreted as aerosol will have a high e value. 
For the extreme case of an atmosphere with no aerosol, 
water- leaving radiance at 765 and 865 nm would result in 
6(765,865) « 2.2, almost double the highest value for a 
true aerosol reported in GW94, and much higher than 
€(765,865) « 1 observed for marine aerosols. The resul- 
tant atmospheric model will lead to an overcorrection er- 
ror where the overcorrection increases with shorter wave- 
length. Where the water already has low reflectance in the 
blue bands, such as in Case-2 water, the overcorrection 
would produce negative Lw in blue bands. Accordingly, a 
solution for t(A)Lw(A r ) is needed. 


9.2.2 Bio-Optical Models for NIR Iteration 

Because of the strong absorption at the wavelengths 
of interest the relationship of remote sensing reflectance 
(Rr 8 ) to the inherent optical properties, backscatter {h) 
and total absorption (a to t) is reduced from the general form 
(adapted from (4) in Gordon et al. 1988): 


-Rrs(A) — y 


T w h{ A) 

Q( A) MA) + a tot (A)’ 


(37) 


where in (37) and (38), T w is the transmission and refrac- 
tion loss at the air-water interface; and Q is the factor 
E u /L u , where E u and L u are the upwelling irradiance and 
radiance, respectively. The Q-factor is often assumed to 
be 7T, although the value is somewhat larger and variable. 
The variable y is the value from Gordon et al. (1988) 
times Q. 

The linear form of (38) assures a stable iteration at 
high reflectance, although (37) shows that i? rs (A) should 
approach {y[T w /Q( A)]) asymptotically without exceeding 
it. Ruddick et al. (2000) also showed that the linear solu- 
tion is an accurate estimator of the Gordon et al. (1988) 
solution in the NIR. 

The backscattering term ( h ), is the sum of h from pure 
water (bbw) and h from particles (hp)- The particle con- 
centration and the scattering efficiency (size characteristics 
and the index of refraction) influence spectral h • Gould et 
al. (1999) determined that the spectral shape of h is linear 
in coastal waters. The backscatter is sufficient to produce 
measurable reflectance in the NIR part of the spectrum 
(Ruddick et al. 2000 and Siegel et al. 2000). 

For the iterative solutions, R rs is estimated at the crit- 
ical NIR wavelengths, A r , from R ra from a reference band, 
Aj, using a solution of (38): 

i?rs(Ar) = /?rs(Ao)^^r bb (A r) A j ), (39) 
G'tOti'V ) 

where the backscatter relationship (rbb) is 


7 'bb(Ar? A j ) 


h(Xr) 

h(Xj 


(40) 


and 77 is a constant. As 2>bp » &bw, it is only necessary to 
determine rbb rather than the actual backscatter. For the 
two NIR bands, A r becomes A*, and A j is 670 nm. 

Gould et al. (1999) concluded that the Petzold volume 
scattering function, 6, varied linearly with wavelength, and 
h & 0.02 b. Their result gives a spectral relationship for 
backscatter in coastal waters: 


b b ( A) = — 0.00113A + 1.62517, (41) 

with 7] — 1 in (40). If no spectral dependence existed, then 
either r) = 0 or B 0 = 0. 

The total absorption is 


^tot(A) — flu>(A) + Uph(A) + &dg(A), (4 2) 

where a w , a p h, and a, i g are the coefficients of absorption 
due to water, phytoplankton, and dissolved or detritus 
matter, respectively. For the 670 nm band, all three terms 
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are determined; for the 765 and 865 nm bands, a p h and a dg 
are presumed to be negligible. 

The a w term is determined from tabulated values from 
Palmer and Williams (1974) and Pope and Fry (1997), with 
additional data from Curcio and Petty (1951) and Smith 
and Baker (1981). The a ph term is found using data from 
Morel and Gentili (1991) with tuning of a ph (440) to the 
SeaBAM data set (O’Reilly et al. 1998, and Maritorena 
and O’Reilly 2000): 


The iteration first uses L t ( A r ) as input to GW94 in order 
to solve for Lw{§ 70), as would be done in any nonitera- 
tive method. Then, Lw(6 70) is used through the model 
described in (39) to determine Lw(A r ) [i.e., Lw( 765) and 
Lw(865)]. These are propagated to the top of the atmos- 
phere (correcting Lw for the direct transmittance T(A) de- 
termined from the GW94 atmospheric model). This TOA 
Lw( Ar) is removed from L t (K) which is then input into 
the GW94 computation. 


<z ph (440) = 0.08(C o )°' 65 . (43) 

Using Bricaud et al. (1998), and Lee et al. (1998) a p h is 
obtained: 


0-iteration: 

L t (X r ) is input into GW94, whose output goes into R rs (A v ), 
where (A v ) are the visible wavelengths; i? rs (670) is input 
into (39), whose output R v s (A r ) is 


a ph (A) - a ph (440) [A 0 ( A) + A x { A) lna ph (440)]. (44) 

While a p h is determined everywhere; for C a < l|xgL -1 , it 
becomes negligible as a p \ l /a w < 0.05. 

In coastal waters such as river plumes, detrital and 
gelbstoff absorption can be significant. Absorption at 
400 nm of 1-20 m -1 , corresponds to an absorption at 
670 nm of 0.04-0.7 m -1 , which is greater than a dg (400) in 
some Case-1 waters. An analysis of field data has shown 
that the a d g value may be approximated by using the 555 
and 670 nm bands. For coastal waters, which includes wa- 
ters having extremely high a dg values, the extreme absorp- 
tion in the blue eliminates the information content of the 
412 and 443 nm bands for determining a dg (412) (Carder 
et al. 1999). As a result, a semi- analytical solution is em- 
ployed using the 555 and 670 nm bands, using the ratio of 
R rs (555) and R rs (670), where 

a dg (670) = 0.147 - 0.18X, (45) 


when 


i? rs (555) - R rs (670) 
i? rs (555) 


(46) 


The coefficients in (44) were determined by fitting X to 
derived a dg (Fig. 33). The effectiveness of this method 
is shown from field data where it is dominated by a dg 
(Fig. 34). For areas offshore, a dg (670) is negligible. An ad- 
ditional analysis compared estimated R rs (412) [using (41)- 
(46) in (38)] to observed R rs (412), including water with 
extremely high dissolved absorption. The solution uses 
a dg (412) = a dg (670) exp[0. 013(670 - 412)]. The compar- 
ison shows the ability to estimate a dg (670) is meaningful 
even when using it to extrapolate to 412 nm (Fig. 35). 


9.2. 2.1 NIR Iteration Application 

In the SeaWiFS processing, (39) was implemented us- 
ing an iterative computation of L^(765) and Lw{ 865). 
The procedure is based on the original SeaWiFS process- 
ing code using the GW94 atmospheric model. The goal is 
to remove the Lw(A r ) component from L t (A r ), so that only 
the atmospheric component of L t (A r ) is input into GW94. 


L t ( \ r )i - Lt(X r ) - T(A r )R rs (A r )oF 0 (A r ). (47) 

1-iteration: 

L t ( A r )i is input into GW94, the output of which goes 
into J?rs(A v )i; R r s(670)i is input into (39), whose output 
l?rs(Ar)l Is 


L t ( A r ) 2 = L t (A r )i - T(A r )R r8 (A r )!Fo(A r ). (48) 

If AR rs (765) = R r s(A ir )i — i?rs(Ai r )o > 10 -5 sr -1 then con- 
tinue the iteration; else stop the iteration. 

The iteration is not performed if the initial Rr S (765) < 
5 x 10“ 5 sr _1 or if Lw(670) < 0, and the iteration stops 
when AR rs (765) (the change between iterations) is less 
than 10~ 5 sr _1 . 

The minimum threshold to initiate iteration constrains 
the solution to water having a significant scattering sig- 
nal. Even in oligotrophic Case-1 water, such as the Loop 
Current in the Gulf of Mexico, SeaWiFS has sufficient sen- 
sitivity that two iterations could be performed if some con- 
straint were not applied. The process is run for up to eight 
iterations; although convergence is generally achieved in 
two iterations, and rarely takes more than four. 

The atmospheric error tends to cause a greater error 
at shorter wavelengths, which has a potential impact on 
the iterative solution. The underestimate results in the C a 
and a p h being overestimated. To prevent overestimation 
of a to t(670) and the resultant overcorrection, the C a con- 
centrations in the computation are limited on the first two 
iterations. The C a value is not allowed to be greater than 
10 jxgL -1 as an input to the computation of the first iter- 
ation, and not greater than 20 jxgL -1 in the computation 
of the second iteration. (The C a product is not limited in 
any way.) If Lw(5 55) < 0 during any iteration, then only 
0^(670) is used in the computation of a tot (670), as nei- 
ther a ph (670) nor a dg (670) is determined. In most areas, 
C a and detrital-gelbstoff absorption are not as critical as 
water absorption in determining a to t(b70). In some estu- 
aries and river plumes, detrital and dissolved absorption 
produce a significant effect. 
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Fig. 33. Comparison of X = [R r s (555) - ^*(670)1/7^(555) to a dg (555) derived from spectral model. The 
labels indicate locations of field data: P=Pamlico, T=Tampa Bay, A= Alabama shelf, F=Florida Bay, S=South 
Atlantic Bight, G= Gulf of Mexico, H=North Carolina shelf; the best fit to this solution is (45). The derived 
a dg came from the total absorption, (38), less the phytoplankton absorption, (44). 



Fig. 34. The derived total absorption (a to t) at 412 and 555 nm compared with AC -9. The coastal a dg 
computation applies at these sites and a dg produces over 90% of the total absorption at these bands. The 
derived total absorption (a to t) used R rs with (44) and (45). 
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Fig. 35. Estimated i? rs (412) versus observed i? rs (412) using the Gulf of Mexico-South Atlantic Bight data 
set developed as part of this study. The SeaBAM tuning was applied here. The boxes are samples that did 
not have a 670 nm value, so the coastal a & g could not be used. The letters denoting locations are the same as 
in Fig. 33. The estimate accounts for 83% of the variance in the measured i? rs (412). 


The iterative solutions were implemented in the Sea- 
WiFS Data Analysis System (SeaDAS) version 4 and were 
applied to pixels containing the station taken within one 
day of the sample. The standard calibration was applied. 
The C a value was found using version 2 of the Ocean 
Chlorophyll (OC2) algorithm (Maritorena and O’Reilly 
2000). 

9.3 RESULTS 

SeaWiFS imagery was processed using the iterative 
technique for the northern Gulf of Mexico, South Atlantic 
Bight, and the East Coast. In the coastal zone of the Gulf 
of Mexico, scattering due to suspended particulate mat- 
ter is significant, as well as absorption from dissolved and 
detrital materials. 

The NIR iteration produces a substantial decrease in 
the correlation of L a (865) with Lvv(0 70) (Figs. 36 and 
37). Without the NIR iteration, a strong correlation ex- 
ists between aerosol radiance and water-leaving radiance 
[the smallest values of L a (865) for each L o (670)j. For pix- 
els with the lowest aerosol, L a (865) increased at about 
0.07 Lw( 670) with a noniterative process. After the it- 
eration, most of the correlation has been removed, with 
L a (865) showing a negligible change for Lw{6 70) < 

1.5mW cm~ 2 pm~ 1 sr“ 1 (or R rs (670) < 0.0096 Sr - * 1 ). The 
L a (865) value shows a weak correlation for greater 
Lw(6 70), with L a (865) changing at approximately 0.015 x 
Lw(6 70). These areas are among the highest scattering 
found, including the core Mississippi River plume and Flor- 
ida Bay during major resuspension events (of calcium car- 
bonate sediments). The linearization of the relationship of 


R T8 to the optical properties (65 and a) in (39) produces 
a negligible effect for low R rs (670). Ruddick et al. (2000) 
noted that for the ratio of NIR reflectances, the linearized 
form produced only small errors. 

A comparison of the iterative and noniterative method 
shows that the iterative techniques increase reflectance at 
shorter wavelengths (Fig. 38). The noniterative and Siegel 
process produces a strong spectrally dependent bias against 
field measurements (Fig. 39), while the NIR iteration re- 
duces both the bias and the RMS error against the mear 
sured data (Figs. 39 and 40). With the reduction in the 
spectral component of the bias, the calculated C a value 
from SeaWiFS becomes more consistent with that calcu- 
lated from field radiometry (Fig. 41). 

In order to prevent iteration from occurring unneces- 
sarily in extremely clear waters, iteration does not take 
place if R rs (765) > 5 x 10~ 5 . This suggests that some NIR 
water-leaving radiance may be present in most Case-1 wa- 
ter (Siegel et al. 2000). The iteration is based on the spec- 
tral characteristic of backscatter by inorganic sediments. 
In Case-1 water (which does not have inorganic scatter- 
ing), the spectral backscatter relationship, rbb(A r ,670), 
will be influenced by phytoplankton and may differ from 
that used here. Phytoplankton tends to produce much 
less backscatter than inorganic particles, so the errors in 
using a sediment- based rbb(A r ,670) may be small, even 
though phytoplankton does produce some NIR scattering 
(Siegel et al. 2000). The potential for error, however, in 
high chlorophyll Case-1 waters should be examined. 

The NIR iteration does not remove all negative ra- 
diances, indicating that other factors, probably absorb- 
ing aerosols, are also involved in this error (Figs. 39 and 
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Fig. 36. Comparison of L 0 (865) with a Luw(670) with a standard noniterative process indicating correlation 
of presumed aerosol with water-leaving radiance (in units of mW cm -2 pm -1 sr -1 ). 
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Fig. 37. Comparison of L a (865) with the NIR iteration, compare with Fig. 38 to show the decrease in 
correlation between the presumed aerosol and water-leaving radiance (in units of mW cm -2 pm -1 sr -1 ). 
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Fig. 38. Mean spectra of noniterative (standard), Siegel iteration, and NIR-iteration methods of processing 
SeaWiFS data are shown in the station data from the Gulf of Mexico and South Atlantic Bight. The 64 
stations were taken over a two-year period. 



Wavelength [nm] 


Fig. 39. Bias of SeaWiFS-derived reflectances against measured reflectances taken within one day of over- 
passes for the noniterative (standard), Siegel iteration, and NIR-iteration methods. 



Fig. 40. RMS error between SeaWiFS and field reflectances taken within one day of overpasses for the 
noniterative (standard), Siegel iteration, and NIR-iteration methods. 
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Fig. 41. Comparison of chlorophyll estimated by the Ocean Chlorophyll 4 version 4 algorithm (OC4v4) from 
satellite- and field matchups for the noniterative (standard), Siegel iteration, and NIR-iteration methods. 


40). By improving the aerosol model selection, however, 
the algorithm results in less error, particularly at shorter 
wavelengths. As expected, once the scattering effect was 
removed, the calculated C a from the OC2 algorithm de- 
creased. The NIR iteration produces a substantial de- 
crease in C a value in coastal waters. The C a algorithm 
is based on the ratio, i2 rs (490)/i? r8 (555), and the errors 
reduce i? rs (490) relative to i? rs (555). Removing the scat- 
tering effect removes this error, bringing the C a value to 
reasonable levels, 

9.4 CONCLUSIONS 

This chapter has described an atmospheric correction 
procedure that will extend SeaWiFS products into the 
coastal regions. This procedure provides an improved 
method of estimating the water-leaving radiance and re- 


mote sensing reflectance in turbid coastal waters. This 
procedure is partially coupled with the Gordon and Wang 
solution for scattering and molecular absorption and is 
based on semi-analytic solutions to the spectral behavior 
of the remote sensing reflectance of the water. This iter- 
ation presumed a fixed spectral shape of the backscatter, 
n>b(A r ,670). In water that is dominated by chlorophyll, 
this shape may change somewhat. 

The full correction described here improves the accu- 
racy of retrieval of remote sensing reflectance, which in- 
herently reduces the number of pixels with unacceptably 
low reflectances. It will further permit extension of image 
products into bays and estuaries, however, it will not elim- 
inate all negative radiance. Areas of highly turbid coastal 
water that produce negative radiances at 555 nm will r€>- 
quire additional algorithms to achieve a valid atmospheric 
correction. 
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Chapter 10 


A Comparison of SeaWiFS LAC 
Products from the Third and Fourth 
Reprocessing: Northeast US Ecosystem 

John E. O’Reilly 

NOAA National Marine Fisheries Service 
N arragansett, Rhode Island 

James A. Yoder 
University of Rhode Island 
Narragansett, Rhode Island 

Abstract 

A comprehensive set of 2,524 LAC scenes of the northeast US continental shelf and adjacent waters scanned 
by SeaWiFS from 4 September 1997 to 11 November 2002 was processed using NASA standard algorithms and 
methods employed in the fourth reprocessing. Estimates of chlorophyll a concentration ( C a ) and normalized 
water-leaving radiance ( Lwn ) from the fourth reprocessing are statistically compared with results obtained 
previously from the third reprocessing. Chlorophyll a concentration from the fourth reprocessing is lower than 
those from the third reprocessing in Gulf Stream and Sargasso Sea waters, along the outer continental shelf 
and slope water, and over the deep basins in the Gulf of Maine. In nearshore shelf waters approximately less 
than 50 m, C a from the third and fourth reprocessings are comparable, except over the shoal water on Georges 
Bank, Nantucket Shoals, and the northern nearshore Gulf of Maine, where C a from the fourth reprocessing is 
approximately 1.1-1. 2 times greater than C a from the third reprocessing. The median Lwat( 412) and median 
Lwn{ 443) from the fourth reprocessing are substantially greater than values from the third reprocessing, and 
the frequency of negative water-leaving radiances for the 412 nm and 443 nm bands is significantly lower with 
the fourth reprocessing. Statistical match-up comparisons between SeaWiFS C a and in situ C a indicate that 
the fourth reprocessing improved the accuracy of C a estimates for this region. 


10-1 INTRODUCTION 

Previous chapters in this volume describe a number of 
improvements to the SeaWiFS algorithms and processing 
methods that were incorporated into the fourth reprocess- 
ing, with a focus on results from the global ocean. This 
chapter focuses on high resolution SeaWiFS LAC scenes 
of the northeast US coastal ecosystem, comparing results 
from the third and fourth reprocessing. 

The northeast US coastal ecosystem (NEC) is a good 
representative area for studying changes in SeaWiFS LAC 
products resulting from changes in algorithms and data 
processing methods. This study area (Fig. 42) encom- 
passes several large eutrophic estuaries, highly productive 
continental shelf water (O’Reilly and Zetlin 1998), and 
oligotrophic blue waters of the Gulf Stream and Sargasso 
Sea. SeaWiFS estimates of surface chlorophyll concentra- 
tion (C Q ) in this area span over three orders of magni- 
tude, with most (99.8%) from 0.02-30 mg m“ 3 (Table 10). 


This region encompasses a wide range of bio-optical condi- 
tions (Case-1 and -2 waters), capturing all but the clearest, 
most oligotrophic waters (0.001-0.02 mg m -3 C a ) scanned 
by SeaWiFS. Moreover, this and other regions adjacent to 
densely populated and industrialized areas present formi- 
dable challenges to the atmospheric correction of ocean 
color data and the derivation of accurate water-leaving ra- 
diances that are required for reliable C a estimates. 

10.2 METHODS 

These analyses are based on 2,524 SeaWiFS LAC scenes 
covering the NEC region and the period from 4 Septem- 
ber 1997 through 11 November 2002. Level-1 SeaWiFS 
data were obtained from the NASA DAAC and processed 
to level-2 using SeaDAS (Fu et al. 1998). Level-2 process- 
ing employed the NASA standard algorithms and standard 
(default) SeaDAS processing parameters used in the third 
and fourth reprocessings (McClain 2000, and Chapt. 1 of 
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Table 10. Statistical characteristics of SeaWiFS products from the third and fourth reprocessing. Values in the 
“0.1%” and “99.9%” columns represent the 0.1 and 99.9 percentiles. The percent of data below zero is indicated by 
P neg . The statistics are based on clear pixels (third reprocessing: 3.07 x 10 8 ; fourth reprocessing: 3.52 x 10 s ) present 
in 2,524 SeaWiFS NEC scenes. 



-76 -72 -68 -64 



Longitude 

Fig. 42. The study area and its major oceanographic regions: Middle Atlantic Bight, Georges Bank, Gulf 
of Maine, Scotian Shelf, Bay of Fundy, Slope Water, Gulf Stream, and the Sargasso Sea. In addition, Cape 
Hatteras (C.H.), Nantucket Shoals (N.S.), and the locations of the 50m, 500m, and 4,000m isobaths are 
shown. The mean position of the north wall of the Gulf Stream is illustrated as a dashed line intersecting the 
4,000 m isobath. 
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this volume, respectively). Scenes collected before 15:30 
GMT and after 18:40 GMT were discarded, leaving from 
2-3 scenes each day with some coverage of the NEC region. 
HRPT stations represented in the final data set included 
2,496 from the NASA Goddard Space Flight Center, 22 
from the Bedford Institute of Oceanography, and 6 from 
the Bermuda Biological Station. SeaDAS (version 4.3) was 
also used to map level -2 data to the standard Lambert 
conic projection (1024x1024 pixels with a resolution of 
1.25km per pixel). 

A number of SeaWiFS products were generated with 
SeaDAS: 

a) OC4v4 chlorophyll a concentration; 

b) Normalized water-leaving radiance for six bands, 
i.e., Lwn( 412), Lwn( 443), Lwn( 490), Lwn(P 10), 
L\vn( 555), and Lvw(670); 

c) Aerosol optical depth at 510 nm, i.e., r a (510), and 
865 nm, i.e., r a (865); 

d) Aerosol Angstrom coefficient at 510 nm, a(510); 

e) PAR; and 

f) Level -2 processing flags, L2 flags. 

The L2 flags used to mask (i.e., eliminate) data from sta- 
tistical analyses of all products except PAR were land, 
sunglint, highlt, straylight, cldice, solzen, and 
lowlw. 

SeaWiFS C a estimates were compared (i.e., matchups) 
with in situ C a measured in surface water collected dur- 
ing surveys of the continental shelf region conducted by the 
National Oceanic and Atmospheric Administration 
(NOAA) National Marine Fisheries Service (NMFS). In 
these matchups, SeaWiFS C a was represented as the me- 
dian of retrievals from a 7 x 7 pixel box centered over the in 
situ sampling coordinates, with at least 11 clear pixels, and 
within ±12 hr of the in situ observation. The in situ C a 
measurements were made at sea on 90% acetone extracts 
of pigments from surface (2 m) water samples using the 
Welschmeyer (1994) fluorometric method and a Turner De- 
signs model 10 AU digital fluorometer fitted with nonacid- 
ification optical filters and a lamp kit. The fluorometer 
was calibrated with Turner Designs certified fluorometric 
chlorophyll standard (19.6 pgL -1 C a ). Measurements of 
fluorescence of a secondary standard — Turner Designs’ red 
solid secondary standard for chlorophyll a — were also made 
with each sample reading to correct for minor variations 
in fluorometer sensitivity during surveys. Matchups were 
statistically evaluated by using a Type II linear regression 
(Press and Teukolsky 1992) of log 10 -transformed data. 

10.3 RESULTS 

The relative and cumulative frequency distributions for 
C a and Lwn values from the third and fourth reprocessings 
are illustrated in Fig. 43. Table 10 summarizes the statis- 
tical characteristics of these frequency distributions. The 


fourth reprocessing yielded more clear pixels (3.52 xlO 8 ) 
than the third reprocessing (3.07xl0 8 ), with the same 
number of scenes and L2 flags applied to eliminate sus- 
pect data. This increase is related to the changes of the 
flags and masks that increased the number of valid pixels 
with the fourth reprocessing (Chapt. 6). 

10.3.1 C a 

When considering the entire region, the median C a 
from the fourth reprocessing (0.441 mgm" 3 ) is lower than 
that from the third reprocessing (0.593 mgm" 3 ). A simi- 
lar decrease in the mode and the 0.1 and 99.9 percentiles 
is also evident (Table 10). The relative frequency distri- 
bution for C a from the fourth reprocessing (Fig. 43b) is 
smoother and has a more log-normal shape than the dis- 
tribution from the third reprocessing (Fig. 43a). 

The spatial distribution of the ratio of C a (from the 
fourth to the third reprocessings) is portrayed in Fig. 44. 
In continental shelf waters shallower than 50 m, estimates 
of C a from the third and fourth reprocessings are com- 
parable, except along the Maine coast, Nantucket Shoals, 
and the tidally mixed shoal water on Georges Bank, where 
C a from the fourth reprocessing is approximately 1.1-1. 2 
times greater than C a from the third reprocessing. 

In areas of usually low C ai such as the slope water, Gulf 
Stream, and Sargasso Sea (Fig. 45), the mean C a from the 
fourth reprocessing is lower (0.8) than the mean C a from 
the third reprocessing. In the turbid, chlorophyll-rich areas 
within the Middle Atlantic Bight estuaries and the Bay of 
Fundy, C a from the fourth reprocessing tends to be lower 
than C a from the third reprocessing. 

10.3.2 Normalized Water-Leaving Radiances 

The median Lwn{ 412) from the fourth reprocessing 
was substantially greater than that from the third repro- 
cessing (Table 10; Fig. 43). The median values for both 
Lwtv(443) and Lwn( 490) also increased, while medians 
for L\yn( 510), Lvw(h 55), and ±vwv(6 70) decreased. The 
percent frequency of pixels with negative water-leaving 
radiances (P neg ) for the 412 nm band was substantially 
lower in data from the fourth reprocessing (14%) than 
from the third reprocessing (25%). In a similar manner, 
7 ? ne g (443), P ne g(490), P neg (510), and P neg (5 55) decreased 
from the third to the fourth reprocessing whereas P ne g (670) 
increased from 7.6-10.7% (Table 10). 

The spatial distributions of P neg (412) and P neg (443) 
are shown in Fig. 46. The general pattern of P neg is corre- 
lated with the pattern for the mean C a , with the highest 
Pneg in the nearshore and estuarine areas and the lowest 
Pneg in offshore waters (compare Fig. 45 with Fig. 46). 
The P neg (412) above 50% from the fourth reprocessing is 
restricted to a narrow band adjacent to the coast, whereas 
P neg (412) from the third reprocessing exceeded 50% over 
most of the inner half of the continental shelf area. 
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Fig. 43. The frequency distributions for SeaWiFS data from the third and fourth reprocessing. The relative 
frequency distribution (thick curve) and cumulative frequency distribution (thin curve) are shown in each 
panel: a) C a , third reprocessing; b) C a , fourth reprocessing; c) Lwn{^ 12), third reprocessing; d) Lwn{4 12), 
fourth reprocessing; e) Lwn{ 443), third reprocessing; and f) Lwn(& 43), fourth reprocessing. 
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Fig. 43. (cont.) The frequency distributions for SeaWiFS data from the third and fourth reprocessing. The 
relative frequency distribution (thick curve) and cumulative frequency distribution (thin curve) are shown 
in each panel: g) Lwn(* 90), third reprocessing; h) Lwn( 490), fourth reprocessing; i) Lwtv(510), third 
reprocessing; j) Lwn($ 10), fourth reprocessing; k) Luw(555), third reprocessing; and 1) Lvw(555), fourth 
reprocessing. 
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Fig. 43. (cont.) The frequency distributions for SeaWiFS data from the third and fourth reprocessing. The 
relative frequency distribution (thick curve) and cumulative frequency distribution (thin curve) are shown in 
each panel: m) Lwn( 670), third reprocessing; and n) Lww(670), fourth reprocessing. 



Fig. 44. The spatial distribution of the ratio of mean C a from the fourth reprocessing to mean C a from the 
third reprocessing. The 500 m isobath is represented as a white line. 
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Fig. 45. The spatial distribution of the mean SeaWiFS C a from the fourth reprocessing. The mean C a is 
a geometric mean based on 2,524 SeaWiFS LAC scenes from 4 September 1997 to 11 November 2002. The 
500 m isobath is represented as a white line. 
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Fig. 46. The spatial distribution of the frequency of negative water-leaving radiance retrievals (P ne g ): a) 
P neg (412), third reprocessing; b) P neg (412), fourth reprocessing; c) P neg (443), third reprocessing; and d) 
P neg (443), fourth reprocessing. Ocean areas where no negative water- leaving radiances were detected are 
designated as white (e.g., portions of the Gulf Stream and Sargasso Sea regions). The P neg value is based on 
the clear pixel retrievals in 2,524 SeaWiFS LAC scenes. The 500 m isobath is represented as a white line. 
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Fig. 47. Reduced Major Axis (RMA) regression of SeaWiFS C a with in situ C a \ a) matchup results from the 
third reprocessing; and b) matchup results for data from the fourth reprocessing. The number of matchups 
(N), regression intercept (INT), regression slope (SLOPE), R 2 coefficient of determination (RSQ), and the 
regression line are shown in each plot. The dashed line represents a 1:1 relationship. 


10.3.3 Matchups 

The statistical results for the space-time matchups be- 
tween SeaWiFS and in situ C a (Fig. 47) indicate that the 
accuracy of SeaWiFS C a from the fourth reprocessing is 
superior to that from the third reprocessing (the regres- 
sion slope is closer to 1.0, and the intercept is closer to 
0.0). The coefficient of determination ( R 2 ) also improved 
from 0.566 to 0.599 (Fig. 47). Note that the fourth re- 
processing yielded more matchups than the third, using 
identical match-up criteria and SeaDAS masking flags, be- 
cause more clear pixels were retrieved (Table 10). 

10.4 CONCLUSIONS 

The suite of improvements to algorithms and process- 
ing methods incorporated into the fourth reprocessing have 
markedly improved the accuracy of SeaWiFS LAC esti- 
mates of C a and water-leaving radiances for this region, 
relative to results from the third reprocessing. The use of 


stray light corrected MOBY data in the vicarious calibra- 
tion (Chapt. 3) and other modifications, have resulted in 
substantial decreases in the frequency of negative water- 
leaving radiances for the 412 nm and 443 nm bands. De- 
spite these improvements, the fourth reprocessing still pro- 
duced relatively high P neg (412) (>20%) and P n eg(443) 
(>10%) throughout a large portion of the continental shelf 
area. The occurrence of negative Lwn values is an obvi- 
ous symptom of inadequate atmospheric correction. Fur- 
ther improvements in atmospheric correction schemes are 
required, particularly in the nearshore areas, and such im- 
provements will likely further improve the accuracy of Sea- 
WiFS estimates of C a for this region. 
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Glossary 

ACS Attitude Control System 
AOT Aerosol Optical Thickness 

BBOP Bermuda BioOptics Project 

CVT Calibration and Validation Team 
CZCS Coastal Zone Color Scanner 

DAAC Distributed Active Archive Center 

ET0P02 Earth Topography 2 min grid 
ET0P05 Earth Topography 5 min grid 

GAC Global Area Coverage (SeaWiFS 1 km resolu- 
tion, subsampled to 4 km) 

GPS Global Positioning System 
GSFC Goddard Space Flight Center 

HRPT High Resolution Picture Transmission 

LAC Local Area Coverage (SeaWiFS 1km resolu- 
tion) 

MERIS Medium Resolution Imaging Spectrometer 
MOBY Marine Optical Buoy 

MODIS Moderate Resolution Imaging Spectroradio- 
meter 

NASA National Aeronautics and Space Administra- 
tion 

NDVI Normalized Difference Vegetation Index 
NEC Northeast US Coastal Ecosystem 
NIR Near-Infrared 

NIST National Institute of Standards and Technol- 
ogy 

OCTS Ocean Color and Temperature Scanner 
Orb View- 2 Not an acronym, but the current name for the 
SeaStar satellite. 

PAR Photosynthetically Available Radiation 

QC Quality Control 

RMA Reduced Major Axis 
RMS Root Mean Squared 
ROLO Robotic Lunar Observatory 
RSR Relative Spectral Response 

SBRS Santa Barbara Research Systems 
SeaDAS SeaWiFS Data Analysis System 
SeaWiFS Sea-viewing Wide Field-of-view Sensor 
SIMBIOS Sensor Intercomparison and Merger for Bio- 
logical and Interdisciplinary Oceanic Studies 
SIRCUS Spectral Irradiance and Radiance Responsi- 
tivity Calibrations Using Uniform Standards 

TOA Top of the Atmosphere 

USGS United States Geological Survey 
UTC Coordinated Universal Time 

Symbols 

do Coefficient of the lunar calibration function of time. 
ai Coefficient of the lunar calibration function of time. 
d 2 Coefficient of the lunar calibration function of time. 
dd g Absorption coefficient for detritus and gelbstoff. 
a oz Ozone absorption coefficient. 
a p Absorption coefficient for particles. 


d p h Absorption coefficient for phytoplankton. 

dtot Total absorption coefficient, the sum of da g , a p , and 

0>xu ■ 

a w Absorption coefficient for water. 

A Cloud-surface system albedo. 

A s Albedo of the ocean surface. 

b Scattering function. 
bb Backscatter coefficient. 

&b P Particulate backscatter coefficient. 

&wv Pure water backscatter coefficient. 

C a Chlorophyll a concentration. 

C 0 z Ozone concentration. 

d a Distance of aerosol model. 

Eo Solar flux at the top of the atmosphere. 

E c Solar flux that would reach the surface if the cloud- 
surface system were nonreflecting and nonabsorb- 
ing. 

E 3 Solar flux reaching the ocean surface. 

Es Estimate of daily PAR. 

E u Upwelling radiance. 

Ee Mean extraterrestrial solar irradiance. 

/ Lunar calibration function of time. 
fr Detector temperature correction. 
f /Q Bidirectional reflectance at the water’s surface. 

F Dependence of the cloud-surface system albedo on 
solar zenith angle. 

Fo Solar irradiance. 

Gi Gain for a given band. 

i An index variable for either a given pixel, band, or 
number of iterations. 
i c Centered pixel i. 

j An index variable for a given band. 

k SeaWiFS aerosol model indicator. 

K 2 Detector temperature correction factor. 

Kd (490) Diffuse attenuation coefficient at 490 nm. 

1 1 SeaWiFS aerosol model indicator. 

L(\m) The Lwn values for MOBY. 

L(\s) The Lwn values for SeaWiFS. 

L(\s/m) L(Xs) averaged over the 5x5 pixel subscene and this 
average value is divided by 

La(Xyi) Aerosol path radiance, including Rayleigh- aerosol 
interaction for wavelength A at location i. 

La(Xyi) The computed aerosol path radiance, including 
Rayleigh-aerosol interaction for wavelength A at lo- 
cation i. 

Lr(Xyi) Rayleigh path radiance for wavelength A at location 
i. 

L ra Interaction between molecular and aerosal scatter- 
ing radiance. 

Lt Radiance at the top of the atmosphere. 

Lt(Xyi) Observed TOA radiance for wavelength A at loca- 
tion i. 

L u Upwelling radiance. 

Lw Water-leaving radiance. 

Lwn(X) Normalized water-leaving radiance at wavelength A. 

Lwn{ A r ) Normalized water-leaving radiance in the NIR part 
of the spectrum. 

L'wn Corrected normalized water- leaving radiance. 

Lwn(0~) Normalized water-leaving radiance, just below the 
sea surface. 

Lwn{ 0 + ) Normalized water-leaving radiance, just above the 
sea surface. 
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n w Refractive index of the water. 

N Number of matchups. 

P a Aerosol phase function. 

Pm Molecular phase function. 

Pneg(A) Percent frequency of pixels with negative water- 
leaving radiances. 

Q The factor E u /L u . 

rhb Backscatter relationship. 

R 2 Coefficient of determination. 

Rrs(A) Remote- sensing reflectance at wavelength A. 

S a Spherical albedo. 

t Time. 

£(A) Diffuse transmittance. 
td Clear sky total (diffuse + direct) transmittance. 
tf(6) Fresnel transmittance of the air-sea interface. 
t g Gaseous transmittance. 

tLf(X,i) White-cap radiance, transmitted to the TOA for 
wavelength A at location i . 

£ox(A, i) Oxygen transmittance for wavelength A at location 
i. 

t ox Gaseous transmittance due to ozone. 

£oz(A, i) Ozone transmittance for wavelength A at location i. 
£ wv Gaseous transmittance due to water vapor. 
to Reference time for lunar calibration time series. 

T Detector temperature. 

T(A) Direct transmittance. 

Td Clear sky direct transmittance. 

T re f Detector reference temperature. 

T w Transmission and refraction loss at the air-water 
interface. 

X See (46). 

y The l\ value from Gordon et al. (1988) times Q. 

a Angstrom coefficient. 

e Single-scattering aerosol reflectance ratio. 
e m8 Multiscattering equivalent of e. 

r i Constant. 

9 Sensor zenith angle. 

Go Solar zenith angle. 

9 V Viewing zenith angle. 

A Wavelength. 

A i Nominal wavelength for band i. 

A i/j Radiance ratio of band i to band j. 

A j Reference band wavelength. 

Am The Lwn values for MOBY. 

A r The NIR part of the spectrum. 

A s The Lwn values for SeaWiFS. 

A s/m Mean of the SeaWiFS and MOBY Lwn values. 

At, Wavelength in the visible part of the spectrum. 

p a Aerosol path reflectance at wavelength A. 

Patm Reflectance due to scattering by molecules and 
aerosols in the atmosphere. 
p s Reflectance of the cloud-surface layer. 
pt Reflectance at the top of the atmosphere. 
p[ Corrected p t for gaseous absorption due to ozone. 

T a Aerosol optical thickness. 
r m Optical thickness of molecules. 

Ua Single scattering albedo of aerosols. 
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